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The Proposed ‘“‘Atomic Energy Act of 1945” 


6 LENARY supervision and control, so far 
as the jurisdiction of the United States 
extends, over all sources of atomic energy and over 
all matters connected with research on the trans- 
mutation of atomic species, the production of nu- 
clear fission, and the release of atomic energy.” 


That, in short, is the scope of the May-Johnson 
Bill introduced into both Houses of Congress 
on October 3 (S 1463 and HR 4280). Under this 
Bill the President would appoint a nine-man 
Atomic Energy Commission which would set 
up its own Administrator, Deputy Adminis- 
trator, and necessary staff and working facilities. 
The powers granted the Commission and the 
Administrator are totally unprecedented in a 
democratic form of society during peacetime. 
Under it there would be transferred into the 
custody and control of the Commission all stocks 
of thorium and uranium and their ores, all 
plants, facilities, equipment, and materials for 
the refining, production, or utilization of thorium, 
uranium, or elements higher than uranium in the 
periodical table; and all processes, technical in- 
formation, contracts, agreements, leases, patents, 
patent applications, and discoveries relating to 
thorium or uranium production or utilization. 
Then, at times set by the Commission or the 
Administrator, all persons having any right, title, 
interest, or claim in such property or processes 
or any ores from which such materials may be 
produced must declare these in the form and 
to the extent prescribed by the Commission or 
the Administrator. 

Under this Bill the Administrator ‘‘is author- 
ized to conduct research and experimentation in 
the field of nuclear fission, the transmutation of 
atomic species, and closely related phenomena, 
and to proceed with the development of any 
and all processes or methods for the release of 
atomic energy, and for the exploitation and use 
thereof for military, industrial, scientific, or 
medical purposes: provided, however, that it shall 
be the policy of the Commission and of the 


Administrator, in accord with the objectives of 
this Act, to utilize, encourage, and aid colleges, 
universities, scientific laboratories, hospitals, and 
other governmental, nonprofit, or private insti- 
tutions equipped and staffed to conduct research 
and experimentation in this field.”” Further, 
however, the Bill provides that “it shall be un- 
lawful for any person to conduct research or 
experimentation involving the release of atomic 
energy in amounts deemed and to be prescribed 
by the Commission as constituting a national 
hazard or being of military or industrial value, 
without the consent of the Administrator and 
upon such conditions as the Administrator, with 
the approval of the Commission, may prescribe.” 
To soften this provision there is added, “In 
administering this provision, the Commission 
and the Administrator shall interfere to the least 
possible extent with small-scale experimentation 
in research laboratories of nonprofit institutions.” 

Under ‘‘Security Regulations” the Bill states, 
“The Commission is authorized and directed to 
establish, and to provide for the administration 
of, security regulations governing the collection, 
classification, dissemination, publication, trans- 
mission, handling, and communication” of infor- 
mation and material in this field. The penalties 
for security violations are extraordinary, amount- 
ing to as much as a fine of $300,000 as well as 
imprisonment for 30 years. 

It is inconceivable that a bill which would | 
throttle scientific research in the manner here 
proposed would receive serious consideration by 
the Congress. As Dr. K. T. Compton has said, 
supposing that when the automotive engine was 
in its infancy, similar secrecy restrictions had 
been imposed because of its possible military 
applications; in such a case, one of our greatest 
assets—our overwhelming superiority in me- 
chanical equipment—would have been com- 
pletely lacking in the present war. It behooves 
every scientist to read this Bill carefully and to 
use his influence in keeping open all channels that 
lead to progress in science.—E. H. (Nov. 1, 1945) 
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Gas Flow in the Mass Spectrometer 


RICHARD E. HONIG 


Physics Department, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received May 21, 1945) 


If a mass spectrometer is to be used successfully for gas analyses, the gas flow through the 
system must meet a number of requirements. The three modes of gas flow—molecular, inter- 
mediate, and viscous—are reviewed, and their working ranges are discussed in terms of the 
pressure of the system. The upper limit for molecular flow, through capillaries or through 
small openings, is found to lie at a pressure where the mean free path of the molecules is 
between 15 and 40 times the diameter of such a gas leak. Expressions are developed to show 
that, for molecular flow, the ion beam intensity is independent of molecular weight but is a 
linear function of the sample pressure, and that the principle of linear superposition applies 
to mixture peaks. Since the composition of a gas mixture in a reservoir is found to Fag 
with time (a numerical example is worked out), the necessity for analyses of short duration 
becomes evident. The practical aspects of designing the gas system for molecular flow and for 
efficient ion production are discussed in detail. The various flow requirements are best met 
by a gas leak consisting of a small hole in a diaphragm, or several holes in parallel. A relatively 


high gas pressure within the electron beam boundaries and resulting large ion currents are 
produced by introducing the gas as a molecular beam into a gas-tight ionizing region. 





INTRODUCTION 
HILE the design, construction, and opera- 


tion of mass spectrometers have received 


considerable attention,’ studies of the gas flow 
in these instruments have not been published. 
At first, the spectrometer was used mainly to 
measure isotopic abundance ratios? and to study 
the products of ionization and dissociation by 
electron impact.’ Recently its use has been ex- 
tended to gas analyses,** and it is in this appli- 
cation that a knowledge of gas flow is of para- 
mount importance. An instrument of the Nier 
type® used for gas analyses is shown in Fig. 1. 


1For a review article, including developments until 
1942, see: E. B. Jordan and L. B. Young, “A Short History 
of Isotopes and the Measurement of Their Abundances,” 
J. App. Phys. 13, 526-538 (1942). 

. Che for example: A. O. Nier, “A Mass-Spectrographic 
Study of the Isotopes of Hg, Xe, Kr, Be, I, As, and Cs,” 
Phys. Rev. 52, 933-937 (1938). 

A: O. Nier, “The Isotopic Constitution of Strontium, 
Barium, Bismuth, Thallium, and Mercury,” Phys. Rev. 
54, 275-278 (1938). 

Review Article: D. Rittenberg, ‘Some Applications of 
Mass Spectrometric Analysis to Chemistry,” J. App. 
Phys. 13, 561-569 (1942). 

*J. A. Hipple, “Gas Analysis with the Mass Spec- 
trometer,” J. Tes. Phys. 13, 551-560 (1942). 

4H. Washburn, D. D. Taylor, E. E. Hoskins, and R. V. 
Langmuir, ‘Application of Mass Spectrometry to Analysis 
of Hydrocarbon Gas Mixtures,” Transactions of the Cali- 
fornia Natural Gas Association, Fall Meeting (1940). 

Herbert Hoover, Jr. and Harold Washburn, “Analysis 
of Hydrocarbon Gas Mixtures by Mass Spectrometry,” 
California Oil World (Nov., 1941). 

H. W. Washburn, H. F. Wiley, and S. M. Rock, “‘The 
Mass Spectrometer as an Analytical Tool,”’ Ind. Eng. 
Chem. Anal. Ed. 15, 541-547 (1943). 

H. W. Washburn, H. F. Wiley, S: M. Rock, and C. E. 
Berty, “‘Mass Spectrometry,”’ Ind. Eng. Chem. Anal. Ed. 
17, 74-81 (1945). ; 

§ A. O. Nier, “A Mass Spectrometer for Routine Isotope 
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To perform satisfactory gas analyses with the 
mass spectrometer, the following primary re- 
quirements should be fulfilled. (1) The ion beam 
intensity should be directly proportional to the 
pressure of the gas in the sample, but should be 
independent of its molecular weight. (2) In a gas 
mixture, the presence of one component should 
not affect the peak intensities due to the other 
components, i.e., all mixture peaks are formed by 
the linear superposition of individual intensities. 





Fic. 1. Photograph of M.I.T. spectrometer* used for gas 
analyses on American Petroleum Institute Project 43c 
(“Studies of the Effect of Radioactivity on the Trans- 
formation of Marine Organic Materials into Petroleum 
Hydrocarbons’’). 


rr Measurements,” Rev. Sci. Inst. 11, 212-216 
1940). 
* Constructed by author in partial fulfillment of the re- 
uirements for the degree of Doctor of Philosophy from 
the Massachusetts Institute of Technology. 


JOURNAL OF APPLIED PHYSICS 


4 






















ar-> a @ «=lew 
ey es Smm Hg  __ __* SSA? 
2 

~6 a 
p.b-—— ied. AM a 











IONIZING 
REGION 


RESERVOIR 


Fic. 2. Essential components of a spectrometer gas system. Approximate pressure levels of various 
parts are indicated by dotted lines. 


(3) The gas flow, and as a result the ion beam 
intensity, should remain essentially constant 
during an analysis. 

All of these requirements can essentially be 
fulfilled if molecular flow through a properly 
designed gas system is established by the correct 
choice of pressures. Figure 2 is a schematic 
diagram of the essential components of a spec- 
trometer gas system. Gas contained in reservoir 
1 at pressure p; flows through a suitable leak L 
into the ionizing region 2 of the spectrometer 
proper. In this ionizing region, which is main- 
tained at pressure pe, excitation, ionization, and 
dissociation of the gas molecules are produced 
by electron impact. A certain fraction of the 
ionized molecules and fragments, after being 
accelerated by adjustable electric fields, is col- 
limated by slits to form an ion beam and 
emerges into the tube proper. The ion beam is 
analyzed by magnetic field B, and the com- 
ponents of the resulting mass spectrum are 
measured electrically at the collector end. The 
neutral molecules and fragments are removed, 
through cold trap T, by diffusion pump 3 which 
is at pressure ps3. If p; be chosen so low that the 
corresponding mean free path of the molecules 
is much larger than the effective diameter of the 
leak opening, then the probability of inter- 
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molecular collisions within the leak is negligibly 
small compared to the probability of molecules 
colliding with the walls. The condition that each 
molecule proceeds through the leak individually 
is known as molecular flow. In Fig. 2 the orders 
of magnitude shown for pressures 1, pe, ps, and 
ps represent average values which will permit 
molecular flow as well as proper operating con- 
ditions for the spectrometer tube. 


LAWS OF GAS FLOW 


The laws of molecular flow were first derived 
by Knudsen in his classical work on the flow of 
rarefied gases through tubes and small openings.®” 
Defining Q,, as the time rate of change of the 
pressure-volume product which, at unit pressure, 
becomes equal to the volume V of gas discharged 
per second by a capillary tube or by a small 
opening, the Knudsen equation may be written 
for cylindrical tubes as: 


Qn =d(pV)/dt = 3800(Dr/L)(T/M)*(p: 
=mr7(pi-— ; (1a) 


®M. Knudsen, ‘“‘Molekularstrémung und innere Rei- 
bun eer Ann. d. Physik [IV] 28, 75-131 (1909). 

TM. Kn Knudsen, “Die Molekularstrémun der Gase durch 
Offoun ee vad & Effusion,” Ann. d. Physik [IV] 28, 
999-1 


— po) 


p2) bar-cm*/sec., 
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Fic. 3. Plot of quantity G vs. mean pressure p for the 
flow of O, through a capillary of diameter 0.066 cm. 
Knudsen’s experimental results are shown as circles, while 
the solid curve represents values computed by Eq. (3). 


pl i To 


and for circular openings as: 
Qn = 2860Do(T/M)}(p1— p2) 


=mo(pi1—p2) bar-cm?/sec., (1b) 


where 
Dr =diameter of tube in cm, 
Do=diameter of opening in cm, 
L=length of tube in cm, 
T =absolute temperature in °K, 
M=molecular weight (scale: O% = 16), 
pi=pressure at high side of tube or opening in bars, 
p:= pressure at low side of tube or opening in bars, 
P= (p1+p2)/2 =mean pressure in bars, 
mr = coefficient of molecular flow through tubes, defined 
by Eq. (la), and 
mo= coefficient of molecular flow through openings, de- 
fined by Eq. (ib). 


In deriving Eqs. (1a) and (1b) above, it is 
necessary to assume that the mean free path Ly 
of the molecules is large compared to the tube 
or opening diameter. Knudsen found experi- 
mentally that Eq. (1a) is applicable in the low 
pressure range provided the mean free path, 
corresponding to the mean pressure of the 
system, is greater than about 30 Dr. At high 
pressures (Ly<Dr) the gas flow Q, is viscous 
and follows Poiseuille’s law: 


Q.=(xD1/128nL)p(p1— ps) 
=vp(pi—p2) bar-cm?/sec., (2) 


where »=coefficient of viscosity in poises, and 
v=coefficient of viscous flow, defined by Eq. (2). 

No simple relation exists for the intermediate 
pressure range (Ly=Dr). However, Knudsen 
was able to establish the following semi-empirical 
relation, which, by the proper choice of the two 
constants C,; and C2, could-be made to fit satis- 
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factorily his experimental data for Hz, Oz, and 
CO, in this intermediate region as well as over 
the entire pressure range: 


1+Cip 
Q Pan. cm*/sec., (3) 
Pi-p) 14+Cop 





where 


G=quantity of gas (=pV) flowing through the tube 
per sec. per unit pressure differential, 

Q=Q,+Q,=total flow in bar-cm*/sec., 

C,=(1.00Dr/n)(M/RT)!, 

C2=(1.23Dr/n)(M/RT)}, and 

R=gas constant = 8.31 X10? erg/°K mole. 


It is seen that Eq. (3) reduces to the coefficient 
of molecular flow for p-—>0 since the fraction in 
the first term approaches unity and the second 
term vanishes. For large pressures, the viscous 
flow term becomes predominant, while the 
molecular flow term is reduced to mrC,/C2 
0.8 mr. Klose* and Adzumi’ repeated the 
experiments on flow through capillaries and 
generally confirmed Knudsen’s results, although 
Adzumi’s value of the molecular flow term is 
closer to 0.9 mr for pure gases and is as low as 
0.6 mr for certain gas mixtures. 

An interesting graph is obtained by plotting 
the quantity G as a function of the mean pres- 
sure p. Figure 3 shows Knudsen’s results for Oz 
flowing through a capillary. The region of 
molecular flow is represented by the horizontal 
line, intermediate flow by the minimum region 
and viscous flow by the steep slope on the right. 
If it be arbitrarily assumed that the region of 
molecular flow extends to the point where the 
quantity G has dropped 1 percent below its con- 
stant value (=mr) (see arrow in Fig. 3), this 
upper limit of pressure fmax is such that the 
minimum mean free path (Ly) min=30 Dr. This 
value of Ly lies between 15 Dr and 40 Dr when 
Knudsen’s results for Hz and CO, and Azdumi’s 
for He, CeHe,.C3Hs, and their mixtures are 
plotted. Thus the inequality Ly>Dr, yielding 
molecular flow as mentioned above, is numeri- 
cally quite well defined for pure gases as well as 
gas mixtures. Knudsen’s results’ for gas flow 
through openings do not give rise to a minimum 
region of G for intermediate values of , but here 
again the region of molecular flow is clearly 


8 W. Klose, “Uber die Strémung verdiinnter Gase durch 
Kapillaren,” Ann. d. Physik [V] 11, 73-93 (1931). 

*H. Adzumi, ‘“‘Studies on the Flow of Gaseous Mixtures 
through Capillaries. III. The Flow of Gaseous Mixtures at 
Medium}Pressures.” Bull. Chem. Soc. Japan 12, 292-303 
(1937). 
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defined and extends to a pmax such that 
(La) min=20p. 


GAS FLOW IN THE SPECTROMETER 


The assertion that all the requirements stated 
in the introduction are satisfied by molecular 
flow must now be proved. To show that for 
molecular flow the ion beam intensity is directly 
proportional to the reservoir pressure ; and is 
independent of molecular weight (requirement 
(1)), it is first assumed that the intensity depends 
linearly on the number of ions formed per second. 
This number, in turn, is proportional to the 
density of gas molecules in the ionizing region, 
and therefore to pressure p:. Since, for equi- 
librium, the gas flow into the ionizing region 
(Qiz) must equal the flow out (Qo3), one may 
write, neglecting the presence of the cold trap T: 


Qio=(T/M)*(p1— po) /ri2= Qos 
=(T/M)*(pe—ps)/res, 


_— Yi2, %23=coefficients including the line 
resistances between 1 and 2, and 2 and 3. But 
pi(10 mm)>p.(=10- mm)>p;(—10-* mm). 
Therefore, to a good approximation : 


Po= Pires/Tie. (4) 
Equation (4) indicates that the gas pressure in 
the ionizing region, and, therefore, the ion beam 
intensity, is directly proportional to the reservoir 
pressure and depends only on the geometry of 
the lines, but is independent of molecular weight. 
By generalizing the basic condition for mo- 
lecular flow stated above, it can be shown that 
in a gas mixture the different components do not 
influence each other’s flow (requirement (2)). If 
it be true that, within the leak boundaries, the 
progress of any one molecule is not affected by 
other molecules present in the reservoir, it fol- 
lows that the nature of these other molecules is 
immaterial. Therefore the molecules of one 
species present in a gas mixture will proceed 
through a leak in a manner independent of any 
other components, and their discharge rate is a 
function only of their own partial pressure dif- 
ferential and molecular weight, and of the leak 
geometry—provided, of course, that the sum of 
all partial mean pressures lies below the upper 
limit for molecular flow. Consider a gas mixture 
consisting of components a, --:7, ---m with 
partial reservoir pressures Pia, -**Pii, ***Pins 
whose total amount in the reservoir 


M1N=VN Dd pu=M1Nnid fi 
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where f; represents the mole fraction. Then 


la =f.p1, etc., 
and similarly, because of Eq. (4) 


P2a =fape, etc., 
which yields: 


oo, puede (S) 


Therefore the total discharge rate may be written 
as 


0=5 0:=(TY/ns)E (p1—bs)s/ MA, 


or 


G=Q/(pi— 2) =(T*/ris) De fi/MH. (6) 
Experiments have shown" that for a fixed elec- 
tron energy and a given gas a the intensity J,+ 
of any one of the peaks in its mass spectrum is 
proportional to the electron beam intensity J- 
and the concentration of molecules in the 
ionizing region, and depends on the kind of gas: 


I+ =kI-poaSa, (7) 


where k=a proportionality factor and S,=spe- 
cific ionization, characteristic of gas a. Since in 
a gas mixture the components contribute 
linearly to the pressure in the ionizing region 
(Eq. (5)), it appears that the total intensity of 
any mixed peak 


I*=kI- > poaiSi=kI-pi(re3/riz) > FSi. (8) 


If the specific ionization numbers Sa, ---S;, 

--S, are determined experimentally from Eq. 
(7) by running each component gas separately, 
a mixture of m components can be completely 
analyzed by applying Eq. (8) times. Analysis 
by this procedure has been described in detail 
by Washburn and his collaborators.‘ 

Constancy of gas flow is the third main re- 
quirement for performing proper gas analyses. 
By making the reservoir volume V, sufficiently 
large, the change in reservoir pressure Ap; 
occurring during the course of an analysis may 
readily be kept down to the order of 1 percent 
or less. Thus the changes in ion beam intensity 
due to changes in gas pressure will be of the same 


10 See for example: D. P. Stevenson ant 4.4 A. et? Norma Jr., 
“Tonization and Dissociation by Electron 7 sg 

Butane, Isobutane, and Ethane,” J. Am. C 

1588-1594 (1942). 
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order as the accuracy with which the peaks are 
measured. However, a far more serious consider- 
ation arises from the fact that the rate of molecu- 
lar flow of a gas depends inversely on the square 
root of its molecular weight (Eq. (1)). While it 
is true that the ion beam intensity is independent 
of molecular weight (see Eqs. (4) and (7)), the 
fact remains that the composition of a gas 
mixture in the reservoir will change slowly due 
to the different flow rates of its components. To 
express this change mathematically, consider a 
reservoir of volume V;, containing initially an 
amount of gas p» V; which after a time is reduced 
to ~iV;. Since the rate of gas lost from the 
reservoir equals the rate of flow through the 
system, one may write at time ?: 


—d(p Vi) ‘dt= = Vidpy /dt=Q, 
=(T/M)*(pi— p2)/rie. 


But since poe= pires/Pi2, 


— Vidp,/dt=(T/M)¥(r12—123)p1/r12= C’p1/ M) 
and 
dpi/pi= —C'dt/ Vi:M*= — Cdt/ M}, 
which yields on integration 


In (pi Po) =— Ct/M' 


or 


Pi=po exp (—Ct/M'). (9) 
650 





Applying Eq. (9) to a binary/mixture of com- 
ponents a and 3, one obtains: 


In (p1/Pa)e_ 
In (p:1/Po)s \Ma 


(p1/Po)a= (P1/Po)s(Mo/ Ma, 


To take a rather extreme case, suppose the 
mixture consists initially of equal parts of 
butane (Mg=58) and hydrogen (My=2), and 
it is desired to determine its composition after 
a certain interval. The reservoir pressure of 
hydrogen corresponding to, -say, a 1 percent 
decrease in butane pressure is computed from 
Eq. (10) as 


(P1/Po) x = 0.9908/2)! — 0.947, 


from which the new composition is obtained as 


or 
(10) 





0.990 
—————100=51.1 percent butane 
0.990 +0.947 
and 
0.947 
100 = 48.9 percent hydrogen. 
0.990 +-0.947 


These changes show up very plainly in Fig. 4 
where, over the region considered, the mole 
percent of the components changes linearly with 
the butane pressure drop. To perform an analysis 
accurate to +0.5 percent, it is therefore neces- 
sary to limit (Api/po)z in this case to about 0.5 
percent. While this example represents one of the 
most unfavorable cases encountered, it points 
out clearly the necessity of large reservoirs and 
analyses of short duration. 

It is interesting to note here that the enrich- 
ment in heavier components occurring in the 
reservoir for molecular flow is not encountered 
with pure viscoys flow. Equation (2) shows that 
such flow is independent of molecular weight, 
but is inversely proportional to the viscosity 
coefficient 7, a quantity which can only be deter- 
mined by experiment. It turns out that in the 
case of a gas mixture the components flow out 
together at a rate determined by a new coef- 
ficient. Not only is this coefficient characteristic 
of the mixture, but its value depends in an 
entirely unpredictable manner on the relative 
abundances of the components.*" Since the 


1 W. Schudel, ‘‘Kinetische Gastheorie und Viskositit 
binaérer Gasmischungen,” Monatsbulletin Schweiz. Gas- 
und Wasserfach 22, No. 2, 21-31 (1942). 
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principle of linear superposition does not apply 
in this case, it is not possible to perform analyses 
in the fashion outlined above. Furthermore, it 


must be realized that even though pure viscous 


flow could easily be maintained between reservoir 
and ionizing region, the flow between ionizing 
region and pump would of necessity be molecular. 
Therefore, enrichment in heavier components 
would still occur in the spectrometer. 


DESIGN OF GAS SYSTEM 


Besides the major requirements for proper gas 
flow discussed in detail above, there are a number 
of practical considerations that are encountered 
when the actual design of a gas system is under- 
taken. For a complete gas analysis, provisions 
must be made to admit in turn to the reservoir 
the gas mixture and all the necessary calibrating 
gases at a pressure that is known and can be 
readily controlled. Figure 5 shows in diagram 
form a satisfactory gas-handling system, adapted 
from a design published by P. T. Smith e¢ al.” 
Gas from any one of the bottles is admitted into 
the manifold to the right of stopcock S; at a 
pressure measured on the manometer. The known 
amount of gas contained in the gas pipet between 
Se and S3; is then expanded into reservoir 1. 
After the measurement has been completed, the 


2 P. T. Smith ef al., ““A High Sensitivity Mass Spec- 
trograph with Automatic Recorder,’’ Rev. Sci. Inst. 8, 
51-55 (1937). 
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reservoir,and gas-handling system are evacuated 
through S,. 

The physical dimensions of the gas system are 
determined by the following considerations. The 
cross section of the manifold should be as small 
as possible to minimize dead space, yet not so 
small that appreciable composition changes in 
the unknown can occur as it flows from the 
sample bottle into the pipet. An interior diameter 
of 3 to 4 mm, corresponding to about 1 cm? per 
10 cm of line, appears satisfactory. Stopcock S; 
serves to cut off the manifold when unknown 
samples are run. A reservoir volume, Vi, of 10 
liters is sufficiently large to keep the gas pressure 
drop at about 2 percent per hour when used with 
an appropriate leak. If the gas pipet volume 
V, is made 4 cm? and filled to a pressure of 10 cm 
Hg, measurable on the manometer to better than 
one part in 100, the initial reservoir pressure is 
given by Boyle’s law as pp>=4X10-* cm Hg or 
40 microns. Then the mean pressure p=20 
microns, and it follows that the corresponding 
mean free path Ly is about 0.2 cm for CO: and 
somewhat higher for most other gases. Thus, to 
insure molecular flow, the leak diameter should 
be smaller than about 0.01 cm. This limitation 
suggests the use of an opening, since capillaries 
are not usually available in such small sizes. 
With the upper limit fixed, the actual diameter is 
determined by the gas flow desired. This flow 
should be large enough to yield ample peak inten- 


CALIBRATION GASES 











| To LEAK 


Fic. 5. DiagramYof gas-handling system. Sto 
amounts of unknown samples are run. 
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UNKNOWN SAMPLE 


ock 5S; serves to minimize dead space when small 
eservoir 1 is made up of two five-liter bottles. 
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Fic. 6. Enlarged cross-sectional view (not drawn to 
scale) of ordinary 7-mm Pyrex tube serving as diaphragm 
for opening of 1-mil diameter. Note how wall thickness is 
tapering off toward hole. 


sities above background, yet not be too large 
in order to keep the reservoir pressure drop down 
to a minimum. It turns out that, together with 
the reservoir pressure and volume specified above, 
an opening about 1 mil (= 0.0025 cm) in diameter 
satisfies all requirements. If in certain cases the 
diameter limited by molecular flow should not 
satisfy the discharge requirements, it is an easy 
matter to arrange a number of small-diameter 
leaks in parallel—a device first employed by 
Knudsen.*® 

Although it would be difficult to make an 
opening in an ideal diaphragm, whose wall 
thickness is negligible compared to the hole 
diameter, small openings can be obtained very 
easily in the following manner: ordinary thin- 
walled glass tubing is heated at one end in a 
flame and allowed to shrink while at the same 
time air is blown through. By judicious heating 
and blowing, openings of any desired diameter 
are obtained, of a shape as shown in Fig. 6. While 
the diaphragm is neither plane nor uniform in 
wall thickness, it serves very well-in actual prac- 
tice. In order to keep the hole from plugging, the 
leak is oriented in such a way that the gas always 
flows from the convex to the concave side. 

Of the many considerations involved in the 
proper design and construction of the ionizing 
region, only the gas flow aspect is to be dis- 
cussed here. Since the production of ions depends 
on the gas pressure in the ionizing region, it is 
desirable to have a relatively high pressure po, 
preferably well localized within the electron 
beam boundaries where the actual formation of 
ions takes place. The requirement that the 
number of collisions between molecules (neutral 
as well as ionized) should be negligibly small 
within the ionizing region limits p2 to about 10~* 
mm Hg. On the other hand, the average pressure 
in the spectrometer tube proper, p,, should be 
képt as low as possible in order to minimize the 
number of encounters of the ions, either with one 
another or with stray molecules, over their 
extended path which amounts to about 60 cm 
in the case of the 60° Nier type instrument. In 
practice, p, may be about 10-* mm Hg. To 
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obtain this pressure differential (p2.—),), two 
different methods are utilized. First, the average 
pressure in the ionizing region is raised by in- 
serting a high flow resistance between the region 
and the rest of the tube. Secondly, the local 
pressure within the electron beam boundaries 
can be boosted further by introducing the gas 
as a well-defined intense molecular beam. 

A relatively gas-tight ionizing region is ob- 
tained by the simple construction indicated in 
Fig. 7 (a), where the narrow gap between side 
walls and bottom is needed for electrical insula- 
tion. Gas-tightness can be improved by insulat- 
ing spacers as shown in Fig. 7 (b), in which case 
the gas molecules can only escape through the 
electron and ion slits. However, such insulating 
spacers must be mounted so they cannot ‘‘see”’ 
the ion beam; otherwise the field distortion 
resulting from stray charges collecting on them 
may completely suppress the beam. If the 
spectrometer is to be used for hydrocarbon 
research, it is important that the effects of 
thermal cracking of gases at the hot filament 
should be minimized. To prevent the products 
of thermal cracking from entering the ionizing 
region, the filament should be located in a region 
where the mean pressure is considerably lower 
than po, i.e., close to a pumping lead. Another 
pumping lead should be connected directly to 
the tube proper to keep p, to a minimum value. 

The formation of molecular beams has been 
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Fic. 7. Geometry of ionizing region. The narrow gap 
required for electrical insulation in (a) has been eliminated 
in (b) by the use of insulating spacers. 
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discussed in detail by Fraser!* and reviewed by 
Loeb." It is found that the net rate of molecular 
flow from region a to } through a tube (Fig. 8) 
is given as 

N=}(n.—m)¢PA, (11) 
where ; 


N =total number of molecules per sec., 

na=N,/Va=average molecular density in region @ in 
per cm#, 

n y= N,/V,=average molecular density in region } in 
per cm, 

t=average molecular velocity in cm/sec., 

A =area of tube, and 

P=probability that any molecule arriving at tube will 
pass through it. 


P becomes unity in the limiting case where the 
tube is reduced to a hole in a thin diaphragm, 
and it approaches 4D/3L for long tubes of 
length Z and diameter D. Equation (11) assumes 
that the condition for molecular flow is fulfilled. 
The spacial distribution of the emerging mole- 


se) 














ray 




















Fic. 8. Molecular beam geometry. 


3 R.G. J. Fraser, Molecular Rays (The Macmillan Com- 
pany, New York, 1931), pp. 10-17. - , 

4 L. B. Loeb, The Kinetic Theory of Gases (McGraw-Hill 
_— Company, Inc., New York, 1934), Sec. 79, pp. 301- 
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Fic. 9. Plot of P(@) cos @ as a function of angle 6. The 
solid curve applies to the special case of a short tube 
(L=D), the dotted circle to a thin diaphragm. 


cules may be written as 


dNo=1(na—m)tAP(0) cos @dw, (12) 


where 


dN»=number of molecules emerging per sec. through 
solid cone dw (see Fig. 8), 
@=angle made by dw with normal to A, and 
P(@)=probability for molecules to emerge at angle 0; 
equal to unity for hole in diaphragm. 


Clausing'® has computed and plotted the quan- 
tity P(@) cos @ in polar coordinates for the special 
case of a short tube (L=D) (see Fig. 9) and 
compared it to the cosine law valid for dia- 
phragms. While the total gas flow N through the 
short tube, represented by the area under the 
solid curve, reduces to about one-half of that 
through the diaphragm, the number of molecules 
emitted in a forward direction (dNo) is undi- 
minished. Thus it becomes clear that well- 
collimated molecular beams can be readily 
formed with the help of relatively short tubes. 

Equation (12) is easily modified to yield an 
expression for the molecular concentration 
n(r, 6) existing at any point (r, 6) of region b 
due to the molecular beam. Since dw=ds/r* (see 
Fig. 8), one may write 


n(r, 0) =dN/dV=dN(r, 6)/eds 
= (1/4xr*)(n.—m,)AP(@) cos 8. (13) 


If region b is totally enclosed, the molecules in 
the beam will be scattered by the walls according 
to the cosine law, as first shown by Knudsen.'® 


% P, Clausing, ‘“‘Uber die Strahlformung bei der Moleku- 
larstrémung,”’ Zeits. f. Physik 66, 471-476 (1930). 

16 M. Knudsen, ‘“‘Das Cosinusgesetz in der kinetischen 
Gastheorie,”” Ann. d. Physik [IV] 48, 1113-1121 (1915). 
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After a number of these reflections, the scattered 
molecules will form a diffuse background, of in- 
tensity mo, which is fairly uniform over the entire 
space. Thus the total molecular concentration 
n(r, @)7 at any point (r, @) is due to the original 
beam as well as a scattered background and 
given by . 


n(r, 0)r=(1/4xr’)(n.—m)AP(0) cos 0+. (14) 


It follows that if this total concentration is 
integrated over the entire space and divided by 
V,, the average density m, will result. 

In the mass spectrometer, the ionizing region 
2 corresponds to the region 6, mentioned above 
in the discussion on molecular beams. The rela- 
tive importance of the first and second terms of 
Eq. (14) will depend on the individual spec- 
trometer design. Since part of the molecular con- 
centration, and therefore of the local pressure pe, 
depends inversely on r’, it is obviously desirable 
to bring the molecular beam as close to the elec- 





tron beam boundaries as is possible without 
introducing series distortion of the electric field. 
This is most readily done by introducing the 
beam through the top plate, as already indicated 
in Fig. 7. Since the ion slits in most instruments 
are quite long, it seems advantageous to replace 
the cylindrical inlet tube by a tube whose cross 
section resembles the ion slit. With this geom- 
etry, the density of ions produced will be quite 
uniform over the entire slit length. 

It is a pleasure to thank Professor Clark 
Goodman for reading the manuscript. and 
making many constructive suggestions. Stimu- 
lating discussions with Dr.,C. W. Sheppard 
helped greatly in clarifying several issues. The 
author also wishes to take this opportunity to 
thank Dr. J. A. Hipple, Jr., Westinghouse Re- 
search Laboratories, the late Dr. O. L. Roberts, 
Atlantic Refining Company, Philadelphia, and 
Dr. H. W. Washburn, Consolidated Engineering 
Corporation, for their cooperation during the 
construction of the mass spectrometer. 





66 JF the voice of the people is heard we shall win this peace, but it will be a tougher, longer job 
than winning the war. For war, after all, is a simple business. Its objectives are clear and 
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heartily supported by the overwhelming majority of citizens. We are united by an impressive com- 
mon purpose. : 

“Not so ts the struggle to maintain a stable, happy society. Here the objective, while clear enough, 
has not the dramatic intensity, the immediacy, the urgency to specific individual and collective 
action. ... 

“Science will march on against disease, ignorance and want, giving us an ever higher material 
standard of living, but along with scientific thinking we need ‘an upsurge of sound, constructive, 
humanistic thought.’ * 

“The supreme danger lies in the possibility that a group of men, despairing of the competence 
of their own thinking to meet a particular set of problems and becoming frustrated, will seek cowardly 
escape of individual responsibility by delegating the authority to think for them to a leader or a 
central planning authority whom they permit, without safeguards, to control their action and con- 
duct. This voluntary slavery of the mind is more degrading, more inherently vicious and potentially 
destructive than mere slavery of the person. Its beginnings are insidious. I sense some of them in the 
organization of various pressure groups, in some of which the leaders are permitted to use methods 
no self-respecting member of the group would countenance as an individual.” 

Joun T. TATE—Commencement Address 
Case School of Applied Science—A pril 22, 1945 
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Here and There 








New Appointments 


Dr. Lauchlin M. Currie has been elected vice president 
in charge of research of National Carbon Company, Inc. 
Dr. Currie has been Acting Director of Research since 
1942, except for 15 months during which he was Associate 
Director of the Division of War Research of Columbia 
University. 


Dr. Earl K. Fischer has begun research work as head of 
the Division of Physical Chemistry at the Institute of 
Textile Technology, Charlottesville, Virginia. 


Louis J. Chatten, who resigned July 31 as Director of the 
Radio and Radar Division of the War Production Board 
in Washington, has been appointed vice president and 
— commercial manager of North American Philips 

ompany, Inc. 


Dr. M. L. Redman has been placed in charge of the 
Electron Microscope field engineering group of Radio 
Corporation of America, with headquarters at Camden, 
New Jersey. For the past two years Dr. Redman has been 
engaged in supervising RCA Electron Microscope instal- 
lations in various hospitals and industrial research labora- 
tories throughout the country. 


W. James Lyons, secretary of the Division of High 
Polymer Physics of the American Physica! Society, has 
joined the staff of the Research Laboratory, Firestore Tire 
and,Rubber Company, Akron 17, Ohio. He was formerly 
with the Radio Corporation of America at Camden, New 
Jersey. 


Gerald Pickett, formerly research physicist with the 
Portland Cement Association, has taken a position with 
Kansas State College, Manhattan, Kansas. 


James G. Potter, recently a physicist with the Bell 
Telephone Laboratories, has been appointed head of the 
department of physics of the Agricultural and Mechanical 
College of Texas at College Station, effective September 1. 


Albert Ray Olpin, executive director of the research 
foundation of the Ohio State University, has been elected 
president of the University of Utah. He will take office on 
July 1, 1946. 


W. S. Gillam, professor of agricultural chemistry at 
Purdue University, Lafayette, dine has joined the 
staff of the Midwest Research Institute, Kansas City, 
Missouri, where he will specialize in research in soils, plant 
nutrition, and analytical services. 


John F. G. Hicks, who had been on leave of absence from 
Corning Glass Works to assist in the development of the 
atomic bomb, left in August for Sao Paulo, Brazil, where 
he will be technical and scientific consultant for Cia. 
Vidraria Santa Marina, affiliate of Corning Glass Works of 
South America. 


Cyril S. Smith has resigned from his position as research 
metallurgist for the American Brass Company, Waterbury, 
Connecticut, to accept appointment as professor of 
metallurgy and head of the newly formed Institute for 
the Study of Metals at the University of Chicago. He has 
been on leave of absence from American Brass for war 
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work since early in 1942, at first as Research Supervisor 
for the War Metallurgy Committee of the National Re- 
search Council, and since April, 1943, with the Los Alamos 
atomic bomb project near Sante Fe, New Mexico. 


Paul Cramer, assistant professor of physics at Denison 
University, Ohio, has been appointed chairman of the 
department of mathematics and physics of Huron College, 
South Dakota. 


Howard F. MacMillin has become associated with 
Arthur D. Little, Inc., Cambridge, Massachusetts. For- 
merly president and general manager of the Hydraulic 
Press Manufacturing Company, Mt. Gilead, Ohio, Mr. 
MacMillin will direct the application to industry of 
developments in the pean: press men in and applied 
physics fields by the laboratories of Arthur D. Little, Inc. 


Institute of Textile Technology 


The Institute of Textile Technology, a graduate school 
chartered by the State of Virginia, has announced a building 
program. The Institute will be enlarged and ——— on 
its present site 14 miles west of Charlottesville. It is a 
non-profit corporation, owned by its members and governed 
by a board of trustees selected he the members, consisting 
of executives prominent in the textile industry. Its field of 
service now extends from Maine to Alabama. 

The Institute will accept as fellows only applicants who 
have their Bachelor’s degrees with excellent standing from 
recognized colleges and universities. It contemplates a 
course of two calendar years leading to a Master’s degree 
and four years leading to a Doctor’s degree in the various 
fields of science. Its prime functions are to train men for 
the textile industry and to utilize its staff in research in 
subjects of interest to the textile industry. 


New Institutes at University of Chicago 


Research which led to invention of the atomic bomb will 
be continued at the University of Chicago through two 
new institutes, devoted to the study of nuclear physics and 
metals. The Institute of Nuclear Studies will be headed by 
Samuel K. Allison, professor of physics at the University 
of Chicago, and his staff will include Enrico Fermi, pro- 
fessor of physics, and Harold C. Urey, professor of chem- 
istry. The staff of the Institute of Metals will consist of 
Cyril Smith, director, and Clarence Zener, professor of 
metallurgy. 

In announcing the establishment of these institutes, 
Chancellor Robert Hutchins emphasized that the purpose 
of them is to advance knowledge and not primarily to 
develop the military or industrial applications of nuclear 
research. The University will accept contributions from 
the government and industry for the support of work, but 
plans to finance the central core of the institutes from its 
own funds. 

As a further outgrowth of its work on the atomic bomb, 
the University has established an institute for applying 
the results of nuclear physics research to such problems as 
cancer, heredity, and the aging process. To be known as 
the Institute of Radiobiology and Biophysics, this organ- 
ization will seek also to perfect techniques for the pro- 
tection of workers using radioactive materials. The Insti- 
tute will be headed by Diatone Raymond E. Zirkle. 


Honorary Degree 


The degree of Doctor of Science was conferred on Dr. 
John R. Dunning, professor of physics at Columbia 
University, at the 56th commencement exercises of the 
Nebraska Wesleyan University. Dr. Dunning, an alumnus 
of the university, gave the commencement address. 
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Administrative Changes at Caltech 


Dr. Robert A. Millikan has retired as administrative 
head of the California Institute of Technology, at the age 
of 77. He will continue as vice president of the board of 
trustees and will assist the president of the board in public 
relations work and institutional development. He plans 
to devote much time to his own research and writing. 


Science reports that Dr. Max Mason has retired from the 
Executive Council at California Institute of Technology. 
He will continue as chairman of the Observatory Council 
and expects to resume work on the 200-inch telescope for 
the Palomer Observatory. 


E. C. Watson, professor of physics, California Institute 
of Technology, has been appointed dean of the faculty. 
Other appointments include William V. Houston, chair- 
man, division of physics, mathematics, and electrical 
engineering; and Frederick C. Lindvall, chairman, divi- 
sion of civil and mechanical engineering. 


Death of Eric Loewenstein 


Dr. Erich Loewenstein, 58, biological chemist at Bausch 
& Lomb Optical Company, died at Rochester, New York, 
following six weeks illness. He was engaged in major 
research on colorimetric methods of determining vitamin 
concentration, and had published a number of papers on 
the subject. 


Physical Society of Pittsburgh 


Elected officers of the Physical Society of Pittsburgh for 
the year beginning July 1, 1945, were: Dr. Mary E. Warga, 
University of Pittsburgh, President; Dr. I. Estermann, 
Carnegie Institute of Technology, Vice President; Dr. 
R Mason, Westinghouse Research Laboratories, Sec- 
retary- Treasurer. 


Physics Club of Chicago 


The following are the new officers of The Physics Club 
of Chicago: President, L. I. Bockstahler; Vice President, 
W. W. Sloane; Secretary, H. A. Leedy; Treasurer, Louis 
Robertson. 


Technical Translations 


The R. T. P.’s listed below are available from the 
Durand Reprinting Committee in care of California Insti- 
tute of Technology, Pasadena 4, California. 


2428 The Effect of Rounding-Off on Originally Sharp 
Wolff, E. Leading Edge on the Resistance of Wings. (In- 
nieur-Arch. Vol. IV, No. 6, December, 1933, Part 

» p. 521, Part Il, pp. 533-537.) (10 pages.) 

2429 Detonation and Peroxides in the Internal Com- 

Sokolik bustion Engine. (Acta Physicochemica, USSR Vol. 9, 
Nos. 3-4, 1938, pp. 593-619.) (15 pages.) 

2432 New Curves of MHeat-Energy Transformation, 

Colombi, Ch. (Bull.-tech. Suisse Romande Vol. 66, No. 14, 13.7.40, 
pp. 149-154.) (8 pages.) 

2433 On an Oscillation Phenomenon in Fluids with Stable 

Gortler, H. Density Stratification. (ZAMM Vol. 23, No. 2. 

i April, 1943, pp. 65-71.) (10 pages.) 

2439 The Interaction of Nitrogen with Metals at High 

Gorvisy, A. Temperatures. (Soc. Chim. de Paris 4th Series, Vol. 
X, pp. 159-160.) (1 page.) 

2441 The Electrolytic Welding of Metals. (Gas and Elek- 

CZ. trowarme 1944, Pt. 4, pp. 74-75.) (2 pages.) 


Photographic Striation Methods Applied to the 
Supersonic Wind Tunnel of the E.T.H. Zurich. 
(E.T.H. Proceedings No. 8, 1943, pp. 44-49.) 
(4 pages.) 

2444 Optical Multi-Curve Recorder of the Finnish State 
Aircraft Factory. (Official File.) (7 pages.) 

New Curves of Heat-Energy Transformation. (Bull- 
tech. Suisse Romande Vol. 66, No. 16, 10.8.40, pp. 
169-174.) (8 pages.) 

Practical Applications of the Diffraction of Light on 


Supersonic Waves. (Physikalische Zeitschrift Vol. 36, 
1935, pp. 142-144.) (4 pages.) 


2443 
de Haller, P. 


2451 
Colombi, Ch. 


2453 
Sokoloff, S. J. 
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2454 Regulations for Altitude Flying and Instructions on 
the Use of Respirators. (D. Luft 1205.) (24 pages.) 

2455 The “S.O.30N"’ and its Design Features. (Aviation 
Francaise No. 7, 2ist March, 1945.) (3 pages.) 

2456 Fundamental Problems of Aeroengine Design. 

Gimm, A. ae Vol. 9, No. 12, 1942, pp. 356-360.) 
10 pages.) 

2459 Structure Investigation in Aluminum Metallurgy. 

Schrader, A. (Aluminum Vol. 25, No. 3, March, 1943, pp. 116- 
126.) (16 pages.) 

2465 On the Development of Play in Bolted Joints under 

Dirksen, B. Fatigue. (L.F.F. Vol. 19, No. 4, 6.5.42, pp. 153- 
156.) (7 pages.) 

2466 The S.0.30N. A very Special Prototype. (Aviation 

St. Arnaud, A. Francaise No. 5, 7th March, 1945.) (2 pages.) 

2467 Friction and Leakage Losses in Piston Rings. 

Eweis, M. (Forschung Vol. 6, No. 2, p. 104.) (1 page.) 

2468 The Mechanics of the Plastic Deformation of Mild 


Hohenemser, K. Steel (ZAMM Vol. 12, No. 1, February, 1932, pp. 


» We 1-14.) (17 pages.) 
2470 The Recent Development of Two-Stroke Engines. 
Zeman (Z.V.D. 1. Vol. 87, No. 1/2, 9.1.43, pp. 7-14.) (12 
pages.) 
2471 The Principles of the Magnetophone Method. 
Lubek, H. (AEC Mitteilungen, No. 9, Sept., 1938, pp. 453- 
459.) (11 pages.) 
2472 Buckling Stresses on Rectangular Fixed Plates. 
Schleicher, F. (Mitt. Forschungshutte Outcheftungshutte Kon- 


zerns, Vol. 1, 1931, pp. 186-193.) (14 pages.) 





New Booklets 


The Electronic Corporation of America, 45 West 18th 
Street, New York 11, New York, has published an illus- 
trated booklet, The Amazing Electron, which is available 
to educational leaders, teachers, and students. Illustrated 
in cartoon technique and simply written, the booklet 
makes clear the composition of the atom, the historical 
development of electronic science, and explains its appli- 
cation in radio (including f-m), television, medicine, and 
industry. The booklet can be obtained without charge to 
individuals or for a 5-cent handling charge in quantities of 
20 copies or more. 








The Taber Instrument Corporation of North Tonawanda, 
New York, in their Bulletin 4506, announce their new 
triple cut specimen shear designed for use with the Taber 
V-5 stiffness gauge to assure accurate cutting of test speci- 
mens for uniform and comparable test results. 


The Applied Research Laboratories of Glendale, Cali- 
fornia, and Harry W. Dietert Company, 9330 Roselawn, 
Detroit 4, Michigan, have published a 4-page folder 
describing a new two-meter grating spectrograph. This 
new instrument can be supplied with 36,600 lines per inch 
or 91,500 total grating lines. Copies available free from 
Harry W. Dietert Company. 


The General Radio Experimenter for July and August, 
1945, contains articles on “‘A Heterodyne Frequency Meter 
for 10 to 3000 Megacycles” and ‘‘How Humidity Affects 
Insulation.”’ 12 pages. General Radio Company, 275 Mas- 
sachusetts Avenue, Cambridge 39, Massachusetts. 


North American Philips Company, Inc., 100 East 42nd 
Street, New York 17, New York, has published a 12-page 
booklet entitled Engineering-Design Development of X-Ray 
Spectrometer. The text covers in detail the basic design 
principles involved in the Geiger-Counter x-ray spec- 
trometer, recently developed industrial control tool. Illus- 
trated with drawings and photos. 
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The Use of Direct Current in Induction Electron Accelerators 


W. F. WESTENDORP 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received September 14, 1945) 


The application of a direct-current bias in the magnetic circuit of induction accelerators 
is described. In its simplest form this bias reduces power consumption to approximately 30 
percent; in especially designed machines a reduction of the weight to 50 percent or less is an 
additional advantage. A machine with a closed central core is discussed and the use of orbit 
contraction by saturation over a range of x-ray energies is described. 





INTRODUCTION 


HE induction electron accelerators or beta- 
trons so far described in the literature’~? 

have been machines operating entirely on alter- 
nating current. During one-quarter cycle of the 
current wave, electrons in a toroidal vacuum 
tube are accelerated by the electrical gradient 
associated with the change of magnetic flux and 
are retained in a circular orbit by the magnetic 
field between two pole pieces. The machine is 
utilized only one-quarter of the time even though 
the electrical gradient is present half the time in 
the proper direction. Much can be gained in 
this respect by the use of direct current. Various 
systems will be described utilizing direct current 
with a saving of electrical power of more than 50 
percent or a saving of materials of well over 50 
percent or both. They are described in the order 
of their electrical and mechanical complication. 


I. DIRECT CURRENT ADDED TO THE CON- 
VENTIONAL MACHINE 


The superposition of direct current upon the 
alternating current used in the magnetizing coils 
of the conventional machine allows the reduction 
of the alternating current to 54 percent of its 
original value. It was found experimentally that 
the iron losses are very nearly proportional to 
the square of the flux density and that the d.c. 


1M. Steenbeck, U.S. Patent 2,103,303, Dec. 28, 1937. 

2D. W. Kerst, Phys. Rev. 60, 47-53 (1941). 
(1941), W. Kerst and R. Serber, Phys. Rev. 60, 53- 58 

1 

4D. W. Kerst, Rev. Sci. Inst. 13, 387-394 (1942). 

5 W. F. Westendorp and E. E. Charlton, J. App. Phys. 
16, 581 (1945). 

* D. Iwanenko and I. Pomeranchuk, Phys. Rev. 65, 343 
(1944). 

7 John P. Blewett, paper to appear in Phys. Rev. 
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component of the magnetic flux did not produce 
an increase in the hysteresis losses for a given 
a.c. excitation. This is of importance because the 
hysteresis losses form the major portion of the 
losses in the silicon steel used for the magnetic 
circuit. The a.c. power consumption can therefore 
be approximately reduced to 30 percent and even 
with the additional d.c. power requirement the 
total power consumption will remain below 50 
percent. 

Figure 1 shows a cross section of the conven- 
tional machine. The magnetizing coils form the 
principal electrical circuit elements. The so-called 
2:1 relation* necessary for a circular orbit is 
established by properly proportioning the gaps 
between the pole pieces and centerpieces with 
respect to the main gap in which the vacuum 
tube is located. The flux density near the orbit 
follows the equation B= Bop in which p= R/Ro 
and 1>n>0 for stability of the electron orbit or 
the magnetic focusing of the electrons toward 
the stable circular orbit.* 

None of these fundamental requirements 
interferes with the superposition of direct current 
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Fic. 1. Cross section of conventional design of betatron 


_ in which direct current can be used with a saving of 


power. 


* 2xRoBo=A®o in which Ro is the radius of the electron 
orbit, Bo the flux density at the circular orbit, and A» 
the change in the total flux through the circular orbit, 
from the instant at which Byo=0. 
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Fic. 2. Magnetizing current as a function of time. 


and the corresponding reduction of the alter- 
nating component. However, the d.c. component 
should never exceed the a.c. crest value since 
otherwise By would not go through zero and no 
electron injection and capture could take place. 
Furthermore, it is important to have a con- 
siderable electrical gradient present during 
injection and for this reason it is advisable to 
make the direct current, e.g., 0.866 of the a.c. 
crest value with an electrical gradient during 
injection equal to one-half of the crest gradient 
occurring during the cycle and locating the 
injection point 30° beyond the negative crest of 
the a.c. wave as illustrated in Fig. 2. 

Figure 3 shows the electrical hook-up which 
will allow the superposition of the direct current. 
The additional circuit elements are 


. a reactor; 

. a d.c. generator or any other source of d.c.; 

. a by-pass capacitor; 

. additional capacitors in the capacitor bank to neu- 
tralize the lagging current drawn by the reactor. 


- wn 


Contrary to the first impression the reactor is 
not of the same order of cost and weight as the 
machine itself. It can very economically be built 
as an air-core type and can be either forced air- 
cooled or convection oil-cooled. The weight of 
the copper ‘is a function of the R.M.S. value of 
the total current squared multiplied by the 
inductance, viz. (J?L)*/5, and can be proved to be 
a minimum for such an inductance as will make 
the R.M.S. value of the a.c. component equal to 
the direct current. However, considerably less 
compensating capacity is required in the ca- 
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Fic. 3. Circuit diagram for the application of direct current 
to a conventional machine. 
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pacitor bank for a slightly larger reactor in which 
the alternating current equals half the direct 
current. Under these conditions an economic 
minimum is obtained for the total cost. 

The d.c. source can be a motor generator set 
or a high current copper oxide rectifier such as 
used in the electrochemical industry. 

The by-pass capacitor need not be very large 
since its reactive a.c. drop can be 100 percent 
neutralized by a small fraction of the reactor 
drop obtained through an appropriate tap, or by 
means of a few turns around the outside of the 
non-metallic reactor tank or housing. 


II. ADDITIONAL COILS AND SPECIAL POLE 
PIECES 


An analysis of the system first described shows 
only a saving in power, no saving in materials, 
and also leads to the conclusion that the maxi- 


CENTERPIECE 
MAIN COILS 


GROOVE COILS 





Fic. 4. Machine with grooved pole pieces, extra 
coils and a large centerpiece. 


mum electron energy attainable in a given 
machine is limited by the saturation of the 
centerpiece laminations which upsets the 2:1 
relation and causes orbit contraction. The maxi- 
mum attainable electron energy has not been 
changed by the use of direct current but the 
saturation of the centerpiece is now caused by 
the combined current a.c.+d.c. If the d.c. com- 
ponent in the central flux could be suppressed 
and the d.c. component at and near the orbit 
be retained still higher flux densities could be 
used at the orbit before causing saturation in 
the center. This purpose can be achieved by 
means of direct-current bucking coils located in 
grooves in the pole pieces as illustrated in Fig. 4. 
Measurements on a solid steel model with 
exploring coils and a flux meter indicated the 
necessity of the use of one large centerpiece and 
two air gaps. Again the 2:1 relationship is main- 
tained by the proper proportioning of these two 
air gaps with respect to the main gap. If one 
gap was used in the center of the model the 
measurements clearly indicated a_ horizontal 
magnetic field component above and below the 
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central plane near the electron orbit caused by 
the strong direct current in the bucking coils in 
the grooves. These radial magnetic fields would 
interfere with the injection and capture of elec- 
trons. The construction with the large center- 
piece and the two gaps eliminated these stray 
fields since outér pole pieces and center-piece are 


at the same magnetic potential when the main’ 


magnetizing coils carry zero current. Con- 
sequently no fringing fluxes would pass through 
the vacuum tube during electron injection. 

Figure 5 shows the electrical connections. The 
groove coils carry slightly fewer direct-current 
ampere turns than the main coils since they just 
buck out the d.c. field in the central core. This 
requires the same number of ampere turns as are 
expended in the main gap, whereas the main 
coils also serve to magnetize the yokes, legs, and 
the small gaps between the various parts of the 
magnetic circuit. 

The additional parts are an a.c. blocking filter 
consisting of a reactor and a capacitor, a d.c. 
generator, and the groove coils. The pole-piece 
construction with grooves offers an additional 
mechanical problem. However, it can be shown 
that the weight of the silicon steel for, e.g., a 50- 
Mv machine can be reduced to approximately 50 
percent of the weight of the a.c. design, while for 
higher energies the reduction is even greater. 

The blocking filter can economically be de- 
signed so that the a.c. current through the 
reactor is 50 percent of the direct current. The 
exciting current of the filter can readily be made 
less than 0.5 percent of the direct current. The 
exciting current will flow through the d.c. 
generator and through the groove coils. It will 
have a negligible effect on the location of the 
orbit. 

Figure 6 shows a different circuit which 
permits all the power capacitors to be placed in 
one bank. In this case the groove coils have been 
given a higher number of turns than the main 
coils so that the induced e.m.f. in the groove coils 
equals the voltage across the main coils and 
consequently no alternating current will flow 
through the groove coils. However, in order to 
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Fic. 5. Circuit diagram for the machine of Fig. 4. 
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Fic. 6. Circuit a —- to the machines 
of Figs. 4 and 7. 


keep the d.c. ampere turns at the proper value, 
the groove coils now will have to carry a lower 
direct current. A control resistance and a by-pass 
reactor allow adjustment of the direct current 
through the groove coils. It is possible to design 
a machine in which the reactor and the resistor 
can be omitted; however the d.c. ampere turns 
in the groove coils would over-compensate the 
d.c. ampere turns in the main coils and a d.c. 
flux through the centerpiece would result, causing 
saturation during the negative peak of the a.c. 
current wave. The magnetic circuit of such a 
machine without reactor would be considerably 
larger than the one for which the complete 
diagram of either Fig. 5 or Fig. 6 is carried out. 

In Fig. 6 the conductor loops shown in series 
with the by-pass capacitors across the d.c. control 
resistance and the d.c. generator represent turns 
coupled with the reactor in order to compensate 
the capacitive voltage drop. 


Ill. ELIMINATION OF THE CENTRAL GAPS 


Examination of Fig. 6 leads to an interesting 
consequence of this diagram. The groove coils 
have a larger number of turns than the main 
coils in order to match voltages and are directly 
connected across the main coils as far as the a.c. 
operation is concerned. The 2:1 relation now 
depends no longer on the ratio of magnetic re- 
luctances but is fixed by the turn ratio of the 
groove and main coils, their voltages being 
identical and their magnetic fluxes being in 
inverse proportion to their numbers of turns. A 
fine adjustment of the electron orbit diameter 
as determined by the correct 2:1 relation may be 
made by adjusting the flux ratio by means of an 
additional tapped autotransformer connected 
between main coils and groove coils. This auto- 
transformer is not shown in Fig. 6. 

Figure 7 shows a cross section of a machine in 
which the central gaps are eliminated and Fig. 8 
shows a model. For ease of manufacture the two 
pole pieces and the centerpiece are still three 
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Fic. 7. Machine with closed central core. 


separate units but the gaps are reduced to the 
lowest practical minimum so that a sheet of kraft 
paper or thin pressboard can be used as a sepa- 
rator. The groove coils in this machine will have 
to carry both d.c. and a.c. bucking ampere turns 
since the reluctance of the closed central core re- 
quires very few ampere turns for magnetization 
and therefore an alternating current will flow 
counteracting the effect of the main coils. More 
space is therefore needed for the groove coils. 
However, the elimination of the central gaps re- 
sults in a reduction of the magnetic energy stored 
in the machine at the peak of the current wave. 
This energy is stored principally in the air gaps 
where the product HB is high. By the elimination 
of the central gaps this stored energy can be re- 








Fic. 8. Solid steel model of a betatron with ring-shaped 
pole piece, groove with coil and large centerpiece. Top yoke 
and top pole piece have been removed to show details. 
This model was used to study the field distribution with 
direct current excitation, by means of a pick-up coil and 
a flux meter, 
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duced by 50 percent or more depending on the 
space taken up by the vacuum tube. The final re- 
sult will be an appreciable reduction in the size 
of the capacitor bank. 

The use of the circuit diagram of Fig. 6 in 
conjunction with the machine shown in Fig. 7 
has the added advantage of adjustment of the 
electron and x-ray energy when orbit con- 
traction by means of saturation of the central 
core is used.? A relatively small reduction in the 
direct current in the groove coil by means of the 
control resistor will cause saturation of the 
central core at the crest of the positive half cycle 
of the a.c. wave since the ampere turns of the 
main coil are then insufficiently compensated by 
the groove coil current. The fact, that the direct 
connection of groove coils to main coils tends to 
maintain the 2:1 relation, will cause a distortion 
in the a.c. wave form of the current drawn by the 
groove coils. This distortion is limited by the 
relatively high leakage reactance of this winding 
and does not offset the effect of saturation. 

It is possible to produce saturation of the 
central core over a wide range of values of the 
main a.c.+d.c. operating current. The ratio of 
a.c. to d.c. crest is maintained throughout this 
range at, e.g., 0.866 and the control resistor in 
the groove coil circuit is used to produce satura- 
tion at the crest of the wave for each setting of 
the main current. A wide range of output x-ray 
energies can thus be covered in a simple manner 
without resort to orbit shift circuits involving 
high current vacuum tubes.*® 


CONCLUSION 


Various methods have been described which 
will make possible the construction of induction 
accelerators with peak flux densities at the orbit 
of 6000 to 8000 gausses and correspondingly 
smaller radii resulting in considerable saving of 
materials and electric power. 

For machines of more than 50-Mv electron 
energy such a reduction in the size of the orbit 
will introduce the new problem of the radiation 
by high speed electrons in a circular path.® 
Experience with the 100-Mv induction ac- 
celerator in our laboratory® has clearly brought 
this problem to the foreground as will be de- 
scribed in a forthcoming paper.’ 
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Crucible Failure in the Induction Melting Process 


ANDREW GEMANT AND JOSEPH STICHER 
The Detroit Edison Company, Detroit, Michigan 


(Received June 4, 1945) 


Experience shows that in the course of the induction 
melting process occasional local overheating of the crucible 
takes place, which in extreme cases leads to perforation of 
the latter. An explanation of this occurrence is presented 
on the basis that it is essentially a phenomenon of thermal 
instability. The theory shows that overheating of the 
crucible is more probable the steeper the electrical re- 
sistivity—temperature characteristic of the crucible ma- 


terial, a condition brought about by impurities present. 
Numerical computations, showing the transition from 
stable to unstable operating conditions, are carried out 
for a small experimental and a large industrial furnace. 
The equations develo show how the various operating 
conditions, particularly the frequency of the current in the 
coil, have to be selected in order to insure operation in the 
stable region. 





1. THE PROBLEM 


NDUCTION heating is rapidly developing 

into a highly important tool in engineering 
processes. The versatility of the method is shown 
in a paper by Sherman,! and is emphasized also 
by Baker.? A particular application, with which 
this paper is concerned, is the high frequency 
induction furnace. Until recently this method 
was chiefly used for steel alloys, hence up to 
temperatures of 1600°C, but nowadays it is 
commonly used for alloys requiring tempera- 
tures appreciably beyond that limit. The ques- 
tion of the crucible material is then of con- 
siderable importance; zirconium oxide crucibles, 
with a melting point of around 3000°C, are 
among those that have been developed for this 
purpose. 

Experience has shown that occasional trouble 
arises causéd by local overheating of the zirconia 
crucible, which causes the latter to melt and 
eventually allows the molten metal charge to 
run out. This local disturbance occurs while the 
bulk of the charge and of the crucible is appre- 
ciably below the melting point of the latter. The 
present paper is concerned with a discussion and 
explanation of such crucible failures in induction 
heating. 

In ordinary heating processes, the heat sup- 
plied is constant and the temperature of the 
material being heated rises until the rate of heat 
dissipation is equal to the rate of heat input or, 
in other words, until equilibrium conditions are 
established. The basis for the explanation of the 
phenomenon discussed in the previous paragraph 
lies in the electrical conductivity versus tem- 
perature characteristic of the crucible material. 
This characteristic may be such that the portion 


1V. W. Sherman, Ae. ras Rev. 2, No. 5 (1943). 
2 R. M. Baker, Trans. A.I.E.E. 63, 273 (1944). 
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of the heat that is generated in the crucible may 
not be constant. It may, indeed, as the crucible 
material becomes more and more conducting with 
rise in temperature, increase more rapidly than 
the heat dissipating characteristic. Thus equi- 
librium conditions cannot be reached and the 
temperature continues to rise until the crucible 
material melts. 

It should be mentioned that an explanation as 
here suggested occurred to the authors as an 
analogy to the so-called thermal breakdown of 
electrical insulators. This latter mechanism was 
first quantitatively developed by K. W. Wagner* 
and extended, among others, by Gemant and 
S. Whitehead.*® In the case of insulation it is an 
electric field due to high potential electrodes that 
causes a current to flow, develops heat, and leads 
to failure; in the case of induction melting it is 
the electromagnetic field that by a similar process 
leads to failure. 

In the process of induction melting, conditions 
are generally such that stable operation is pos- 
sible. However, in localized spots of the crucible 
conditions might be different, either owing to 
inhomogeneity of the material or to secondary 
reactions, for instance, reduction of the zirconia. 
If in this localized region conditions of stability 
are not fulfilled, local failure will be the result, 
even though the bulk temperature of the crucible 
is below its melting point. 

Although in practice the failure is probably 
always a localized process, the mathematical 
treatment given herein will be based on homo- 
geneous conditions extended over the whole cir- 
cumference. The reason for this postulate is that 
the results probably do not differ essentially in 

3K. W. Wagner, J. Am. Inst. Elec. Eng. 41, 1034 (1922). 


*K. W. Wagner and A. Gemant, Preuss. Akad. Ber. 
(1934). 


( . al Gemant and S. Whitehead, Phil. Mag. 27, 582 
1 j 
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Fic. 1. Cross-sectional diagram of induction 
heating furnace. 


the two cases, and that the simplified treatment 
as presented will be sufficient to show the prin- 
cipal features of the mechanism. Besides, the 
treatment remains simple only so long as circular 
symmetry holds, becoming too involved for 
simple mathematical analysis when conditions 
vary along the circumference of the crucible. 

This presentation shows in a semi-quantitative 
manner the conditions that must be fulfilled to 
make stable operation possible. One of the 
essential characteristics is the frequency of the 
applied current, and the results indicate how 
the conditions of stability vary with the fre- 
quency. This circumstance permits the selection 
of a frequency range assuring optimum operating 
conditions. 


2. THE CRITERION FOR THERMAL 
STABILITY 


The treatment of the problem given refers to 
a cylindrical charge arranged coaxially with a 
cylindrical coil (Fig. 1). In the cross-sectional 
figure, Ch indicates the charge, Co the coil, and 
the space Cr between them the crucible. The 
radius of the metal charge is R, that of the inner 
surface of the coil is R’. The space Cr is either 
completely occupied by the crucible, or the 
crucible is only the inner part of the space, as 
indicated by the dotted circle, the space outside 
the crucible being filled with a granular mass 
acting as spacer and consisting of ZrO2 pebbles 
or the like. The wall thickness of the crucible in 
the former case is R’ — R, in the latter b. 

In the subsequent analysis it is assumed that 
the device is sufficiently long for end effects to 
be negligible in a first approximation. This is 
more justified, the better the heat insulation at 
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thetop and bottom of the furnace. Hence, heat 
flow is assumed to take place only in a radial 
direction. All data concerning heat generated 
and dissipated refer, therefore, to unit axial 
length (cm). 

The consideration of transient conditions is 
outside the scope of this paper, in other words, 
the process of the gradual heating of the charge 
is not considered. The question is only whether 
or not conditions permit stable equilibrium. 
=a: | Wo quantities will have to be computed, the 
heat, W,, generated per second per cm length 
and the heat, Wa, dissipated per second per cm 
axial length, in order to determine under what 
conditions W,= Wag, the criterion of stability. 
In making use of the theory of induction heating, 
the reader is referred to various texts, or to 
articles such as that by Greig.* The induction 
melting process in particular is discussed by 
Wall.’ 

The quantity W, is the sum of two terms, the 
heat generated in the charge, Wen, and in the 
crucible material itself, Wce,. The heat generated 
in the coil does not have to be considered in this 
connection. 

For the value of Wc, the following equation 
will be used : 


Wen=(RH¢?/8)(focn)', (1) 


where Hy,=amplitude of the magnetic force in 
the empty coil (in oersteds), f=frequency of 
applied current in coil, pc, =electrical resistivity 
of charge (electromagnetic units). This relation 
is valid only if the condition 


tR(8f/pcr)'>4 (2) 


holds. Expressed in words, formula (2) states 
that the ratio of the charge radius to the depth 
of penetration of the magnetic field should be 
larger than, say, 4. Conditions of operation are 
nearly always such that (2) is satisfied, and, 
accordingly, Eq. (1) can be used. 

The permeability in Eq. (1) as in all following 
equations is taken as unity. This is correct for 
non-magnetic charges, such as non-ferrous metals 
and ferrous metals at temperatures above the 
Curie point (about 800°C). 

Concerning W¢,, it has to be pointed out that 
in general the total temperature drop takes place 
across the crucible, since the coil in practice is 
usually cooled to a relatively low temperature. 
While the electrical resistivity, p, of the crucible 


* James Greig, World Power 25, 246 (1936). 
7 T. F. Wall, The Engineer 173, 209 (1942). 
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is very high at low temperatures, it falls to low 
values at high temperatures (more will be said 
about this relation in Section 3). The current 
flowing in the crucible is, therefore, negligibly 
small except for a narrow region adjacent to the 
charge that has a temperature only slightly 
below that of the metal. A rigorous treatment 
would have to start with the differential equation 
of the heat balance in each annulus of the 
crucible. A good approximation could be carried 
out on the basis of the logarithmic drop of tem- 
perature across the crucible, but the present 
treatment is even further simplified. 

It is assumed that a thin annulus adjoining the 
charge is the only heat source due to the crucible. 
This amount of heat, Wc,, is then added to Wea 
and the total is considered as a cylindrical heat 
source acting at the boundary between the charge 
and the crucible. The bulk of the crucible is con- 
sidered free of heat sources from the standpoint 
of heat conduction. In the case of a thin crucible 
in a granular mass the thin annular source is 
simply identified with the crucible; the thickness 
of the source is then b (Fig. 1). In the case of a 
thick-walled crucible the thickness of the source 
(also designated by 0b) is taken arbitrarily. The 
value of } in this case might also be estimated by 
inspecting the depth of discoloration that can 
be observed at the inner wall of crucibles. 

Since the resistivity of the crucible, p, at any 
point, even at the inner surface, is appreciably 
larger than pcx, it is probably correct to say 
that the magnetic field in the crucible is Ho, the 
same as in a coil unaffected by the charge. This 
statement is utilized in computing We,. 

If i9 denotes the maximum of current density 
in the crucible material, then 72p/2 is the heat 
generated per cc. As to the current density, 7: 


t= —(do/dt)/24Rp (3) 


(with permeability unity), where ¢ is the flux 
linked with the annulus of radius R and thick- 
ness b. 

The flux ¢ is the sum of two terms, the flux in 
the charge, ¢ca, and the flux in the thin layer 
of crucible, ¢c,. For the first it is assumed that, 
because of the large p/pcx ratio, its value is not 
affected by the presence of the crucible; we then 
have (with t=time) : 





ocra= 


HyR PCh ; 
(>) (cos 2xft+sin 2xft). (4) 


This relation is valid if condition (2) holds; 
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otherwise certain Bessel functions have to be 
considered. 


~ Concerning ¢c,, the field intensity in the 
crucible is Hy cos 2xft. Evidently for the inner- 
most layer of the annulus b this flux component 
is zero; for the outermost layer it is 2rRbH, 
Xcos 2xft. Taking the arithmetic average : 


¢cr=mRbH, cos 2zft. (5) 
The sum of Eqs. (4) and (5) is the total flux 
linked with the annulus. From Eq. (3), 
° Ho 
4 es, “ (perf )* cos 2nft 
p 
+((penf)!+2xfb) sin 2xft]. (6) 


Adding the two vector components : 
. Ho 
to= > Lecaf+ (oenf)'+ 2x fb)*}}. (7) 
p 


Computing %9’p/2, the heat generated per cc 
becomes : 


H?k/8p, 
with , 
k=penft ((ecaf)*+2nfb)?. (8) 
Hence, 
Wer=aRbH ek /Ap. (9) 


The total heat generated, W,, is then the sum 
of Eqs. (1) and (9): 


RH, abk or) 
4 p 


v 








(10) 


The next step is the calculation of Wa, the 
heat dissipated radially across the crucible 
toward the coil. We have 


Wa=G(T-T>), (11) 


where 7 =temperature of the charge, 7)>=tem- 
perature of the coil, both in absolute degrees, and 
G=thermal conductance of the crucible. For G 
we have 


G=(g/2m) log (R’/R), (12) 


with .g=thermal conductivity of the crucible 
material. 

In the following analysis pcs, the electrical 
resistivity of the metal, is taken as a constant, 
its moderate increase with temperature being 
disregarded. The temperature dependence of p 
on the other hand is considerable, and will be 
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Fic. 2. Specific resistivity of zirconia as a function of 
temperature. Dotted line: straight extension of experi- 
mentally known portion; (1) and (2): assumed breaks, 
according to information on other inorganic crystals. 


represented by the equation 
p=po exp (8/T), (13) 


with pp and 6 =constants. 

In order to determine the criterion of stability, 
both W, and W, are considered as functions of 
temperature. If plotted in this manner, one single 
curve results for Wz, but a whole family of curves 
for W,, containing Hy as the parameter. Most of 
the W, curves show two intersections with W,: 
the one, corresponding to the lower temperature, 
represents the stable state, i.e., the equilibrium 
temperature at that particular value of Ho; the 
other is of no significance, since it determines an 
unstable equilibrium. The W, curve is being 
shifted toward higher values with increasing Hp. 
If a certain critical value of the latter which may 
be called the limiting field H; is reached, W, and 
Wz will just touch at one point. For still higher 
fields there will be no common points between 
the two curves, i.e., the temperature rises con- 
tinuously without ever reaching an equilibrium. 

‘From this analysis the method of finding H; is 
evident. It must fullfill two conditions, one being 


W,=Wa, 
the other 
dW,/dT =dW./dT. 
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These two equations contain two unknowns: H; 
and 7; (the limiting temperature) which can be 
calculated. 

The two equations are : 


—( arbk (perf)? 





) =G(Ti-— To), (14) 








4 \po exp (8/T:) 2 
and 
aBRH *bk 
=G (15) 
4T ;* po exp (8/T)) 
Eliminating H;, one obtains 
po exp (8/T1) (perf)! B(T:—To) 
teweecerar, fot 
2xbk T? 


which allows computation of 7). The limiting 
field is then from Eq. (15) 


j= 





4T :* po T)\3 
( po exp (8/ ’) a7) 


aBRbkG 


Equation (17) in conjunction with Eqs. (8) 
and (16) contains the final result of our analysis. 


3. THE ELECTRICAL RESISTIVITY OF 
ZIRCONIA 


Of the various quantities occurring in Eq. (17), 
the frequency of the applied current plays im- 
plicitly an essential part, as will become evident 
from subsequent numerical computations. As 
explained, 0 is arbitrary, but this weakness could 
be remedied by a more elaborate treatment. 

The most important constant is perhaps 8, 
introduced by Eq. (13), since it occurs in the 
exponent. Unfortunately, its value is not known 
from experiment, and for this reason this subject 
deserves a brief discussion. 

The electrical conductivity of zirconia has 
been measured only up to 1500°C, as far as the 
authors could ascertain. The few existing data‘ 
are plotted in Fig. 2. As pointed out, however, 
recent applications extend up to 1800°-1900°C, 
and it is rather unlikely that the same 8 (slope) 
that is valid up to 1500° should hold up to higher 
temperatures (dotted line). 

There is information available® on a large 


8 International Critical Tables (McGraw-Hill Book 
Company, Inc., New York, 1927), Vol. II, p. 86. 

® Landolt-Bérnstein, Physikalisch-Chemtsche Tabellen 
(Verlagsbuchhandlung Julius Springer, Berlin, 1936), 
Complem. Vol. IIIc, p. 2017. 
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number of lower melting crystals for which the 
data cover the whole temperature range up to 
the melting point. These crystals are essentially 
the halides of alkali metals, silver, thallium, and 
the like. It appears that at temperatures of 
several hundred degrees C below the melting 
point the slopes of the curves suddenly increase 
to values of about twice the previous values. 
Lacking experimental data on ZrO, for tem- 
peratures above 1500°C, two arbitrary exten- 
sions (lines 1 and 2) were assumed as possible 
bases for numerical calculations. Because of 
these assumptions the following numerical data 
can be considered as only semi-quantitative. A 
research filling this gap in our knowledge would 
be very valuable. 

It is known that admixtures to ZrO, lower the 
resistivity considerably (the values for the Nernst 
glower—essentially ZrO» with a few other oxides 
—for instance, are 1000 times lower than those 
for pure zirconia). The actual p vs. T curve will, 
then, depend on the degree of purity of the 
material, the resistivity being lower the more 
impurities are present. The lines 1 and 2 repre- 
sent two possible curves corresponding to two 
different degrees of contamination. Secondary 
processes, as mentioned in Section 1, result in a 
lowering of the p vs. T curve, and this, in its 
turn, influences the values of po and 8 in Eqs. 
(13) and (17). It will be shown in the following 
section how the value of these constants influ- 
ences the conditions determining stability. 


4. NUMERICAL COMPUTATIONS FOR TWO 
INDUCTION MELTING FURNACES 


In the following, two special cases are dis- 
cussed and the equations deduced above are 
used in determining the limiting field intensity, 
beyond which stability is not possible. The 
ampere-turns (r.m.s. value) per cm axial length, 
AS, are related with Hy by the equation 


AS=H,/0.44v2, (18) 


and the final results will be presented in ampere- 
turns. 

The numerical values for those various con- 
stants that occur commonly in both cases are as 
follows: One of the higher melting metals for 
which sufficient data are available in the litera- 
ture is tungsten, and this metal is used, therefore, 
in these examples. The electrical resistivity of 
tungsten, pcr, at around 2000°C is 60 10-* ohm 
cm or 60X10* electromagnetic units. The mag- 
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Fic. 3. Heat generated and dissipated vs. temperature 
for a crucible with a radius of 2 cm, a frequency“of 4X 10 
c.p.s., and 78 amp.-turns/cm. Curve (1) shows stable, 
curve (2) unstable operating conditions. 


nitude of b is taken 0.3 cm, this being the order 
of magnitude of the wall thickness of small size 
crucibles. The outer surface of the container, 
adjoining the tubes carrying cooling water, has 
a temperature of, say, 100°C or, in a rounded 
number, 7»)=400°K. 

Concerning p, the resistivity of ZrO, slopes 
(1) and (2) of Fig. 2 are used, as pointed out 
above. For curve (1), po is 0.32 electromagnetic 
unit, B=53X10*, for curve (2), po is 0.810 
electromagnetic unit, 8=65 X 10°. 

(a) The first case under consideration is an 
experimental type of furnace, utilizing a small 
zirconia crucible of 2 cm radius, surrounded by 
zirconia pebbles in a thickness of 6 cm (R=2, 
R’=8). The thermal conductivity, g, of the 
zirconia pebbles was taken as 0.8X10° ergs per 
sec. per cm per deg. C or 0.5 B.t.u. per hr. per 
ft. per deg. F, as estimated from data on zirconia 
brick.” From Eq. (12), G=3.7X 10° ergs per sec. 
per cm per deg. C. 

The results of a detailed calculation are shown 
in Fig. 3, referring to a frequency of 4X10° 
c.p.s., a high value suited to the small size of 
the charge. The curve shows watts/cm vs. tem- 


1% Max Jakob and G. A. Hawkins, Elements of Heat 
Transfer and Insulation (John Wiley & Sons, Inc., New 
York, 1942). 
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Fic. 4. Heat generated and dissipated vs. temperature 
for a crucible with a radius of 25 cm, a frequency of 4X 10° 


c.p.s. and 260 amp.-turns/cm. Curve (1) shows stable, 
curve (2) unstable operating conditions. 


perature in centigrades. The line Wz is the heat 
dissipated. For W, branch (1) in Fig. 2 was 
assumed. Calculating the limiting field and the 
limiting ampere-turns, AS), a value for the latter 
of 80 was obtained from Eqs. (17) and (18). 
This means that stable operation is possible only 
below that value. As an example, AS=78 was 
selected—a value slightly below 80—and the 
curve (1) as shown was calculated from Eq. (10). 
The left-hand intersection indicates the condition 
for equilibrium. The steady-state temperature 
under these conditions would be somewhat 
below 2200°C. 

Suppose now that at one localized spot the 
crucible contains impurities that cause the p—T 
characteristic to follow branch (2) of Fig. 2. 
Calculating W, for such a case (with AS=78), 
curve (2) of Fig. 3 results, showing that no stable 
condition is possible. The temperature at that 
localized spot will rise until the crucible melts and 
and the charge runs out. Calculating AS; for this 
case, a value of 77 is obtained. Hence, with our 
ampere-turns of 78, it is obvious that we are 
now in the unstable region. 

(b) The second case refers to a large, industrial 
furnace, having a charge radius of 25 cm and a 
crucible thickness of 8 cm. The thermal con- 
ductivity of the solid crucible is taken as 1.7 K 105 
ergs per sec. per cm per deg. C or 1.0 B.t.u. per 
hr. per ft. per deg. F. From Eq. (12), G=3.8X 10° 
ergs per sec. per cm per deg. C. 

The results for this case are shown in Fig. 4 


666 





which is analogous to Fig. 3. The frequency is 
410° c.p.s. With branch (1) of the p—T 
characteristic (Fig. 2), AS;=270 amp.-turns/cm. 
Thus, we apply only 260 amp.-turns, and obtain 
curve (1), indicating stable operation. Assuming 
branch (2) of Fig. 2 to be valid at a localized 
spot, curve (2) results, having no intersections 
with the Wz, curve. AS; for this latter case turns 
out to be 250 amp.-turns/cm. 

When the failure point is reached in a furnace 
of this type, pitting of the inner surface of the 
crucible, rather than perforation takes place. 
This process, however, might be cumulative and 
eventually leads to a destruction of the crucible. 


5. PRACTICAL CONCLUSIONS 


Practical conclusions from the analysis pre- 
sented in this paper can be discussed from two 
angles. 

(1) If the explanation of crucible failure, as 
given, is correct, it follows that the main reason 
for a failure is the steepness of the resistance— 
temperature characteristic of the crucible ma- 
terial at temperatures above 1600°C. This again 
is influenced by the degree of purity of the 
zirconia. These impurities are either originally 
present in the material or are introduced as a 
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Fic. 5. Limiting amp.-turns/cm and limiting watts/cm vs. 
frequency for crucibles of 2- and 25-cm inner radius. 
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result of secondary chemicai reactions during the 
melting process. 

If troubles of this type are to be avoided, the 
quality of the zirconia has to be improved by 
proper selection of the raw materials and manu- 
facturing processes, and operating conditions 
have to be selected so as to prevent the occur- 
rence of unwanted reactions. 

(2) If an improvement along these lines is 
not feasible or is otherwise unwarranted, the 
process has to be carried out under conditions 
permitting a high value of AS;. The role of the 
various factors is evident from Eqs. (16) and 
(17) and a systematic analysis would show the 
effect of each factor upon the value of AS, 
quantitatively. 

As an example, the role of the frequency, one 
of the most essential factors, is shown in Fig. 5, 
both for the experimental set-up (R=2) and the 
large furnace (R=25). The abscissa gives the 
logarithm of the frequency which is chosen as 
the independent variable. The ordinates show 
(left side) the respective values of AS), the 
limiting ampere-turns/cm, and (right side) W,, 
the corresponding watts/cm generated in the 


charge, calculated from Eq. (1). This latter 
quantity determines the speed of the melting 
process, the time to reach a given temperature 
being shorter the larger the heat input rate. 
The range of possible frequencies used is 
limited toward the lower end of the range by 
conditions imposed by the depth of penetration. 
These conditions are known and need not be 
discussed here. The present analysis shows that 
there is equally a limit toward the high end of 
the frequency range. The higher f, the lower both 
the limiting ampere-turns/cm and the limiting 
heat input W:. The conclusion can be drawn 
from Fig. 5 that the lowest frequency com- 
patible with the depth of penetration criterion 
permits the highest limiting input. If, however, 
other considerations suggest the use of higher 
frequencies, as in practice often might be the 
case, the graph shows how far the ampere-turns 
can be increased without damage to the crucible. 
It is evident from the section dealing with the 
resistivity of zirconia that the graphs presented 
in this paper will have to be recalculated when 
true resistivity—temperature characteristics of 
the crucible material are obtained experimentally. 





A Double-Solenoid Coil System for Magnetic Measurements 


E. W. FRoweE Anp C. J. ARONSON 
Naval Ordnance Laboratory, Navy Yard, Washington, D. C. 


(Received August 2, 1945) 


A double-solenoid coil system is described which can be used as a neutralizing coil having 
large internal and small external magnetic fields. Its application in the measurement of mag- 
netic moments of small magnets is also discussed. gag 


INTRODUCTION 


EXPERIMENTAL techniques occasionally 

e necessitate the use of current-carrying coils 
to neutralize magnetic fields within very limited 
volumes. Since it is often not desirable for the 
effect of the neutralizing coil to extend appreci- 
ably beyond the volume where the field is to be 
neutralized, an ordinary solenoid coil cannot 
always be used as a neutralizing coil. 

2. One of the solutions to the problem of ob- 
taining a neutralizing coil having a limited 
extent of external field employs a toroid-shaped 
coil to provide the neutralizing field. In Fig. 1 
volume C represents a region within the toroid 
in which the component of field parallel to the 
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y axis is to be neutralized. The distances A and 
B are measured from the center of the volume C 
to the top and side of the toroid, as shown. The 
volume C is small enough, so that the curvature 
of the neutralizing field of the toroid may be 
neglected. The principal disadvantage of such a 
system is that to produce a uniform neutralizing 
field in the volume C, it is necessary to have a 
comparatively large toroid, i.e., A>B, pre- 
cluding the use of such a coil within small spaces. 
A second disadvantage is that the toroid neu- 
tralizes only the component of field parallel to 
the y axis, as shown. If this component is 
parallel to a surface, it can be neutralized at a 
distance B from the surface; but if this com- 
ponent is perpendicular to a surface, it cannot be 
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neutralized closer than the distance A from the 
surface. A third disadvantage of the toroid is 





Fic. 1. 


that such a coil is difficult to wind. The double 


solenoid was developed to overcome these disad- 
vantages of the toroidal neutralizing coil. 


DESCRIPTION OF THE DOUBLE SOLENOID 


3. The double-solenoid neutralizing coil which 
is shown schematically in Fig. 2 consists of two 
coaxial, solenoidal windings of the same length, 
L, but with different diameters and different 
numbers of turns. The two solenoids are con- 
nected in series, so that they produce magnetic 
fields of opposite polarity. The design is such 
that the magnetic moment (product of the total 
number of turns times the cross-sectional area) 
of one solenoid is equal to the magnetic moment 
of the other solenoid. Thus, the external field of 
the system decreases rapidly with distance from 
the center and becomes negligible for most pur- 
poses even at distances as small as one solenoid 
length from the center of the coils. The internal 
field, however, is not zero, but is equal to the 
difference in the fields of the two coils. The inner 
coil, having the greater number of turns per 
centimeter and the smaller radius, produces the 
stronger field at the center of the system. 


FIELD AND MOMENT CALCULATIONS 


4. Computations of the internal field of the 
double-solenoid neutralizing coil are based on 
the equation for the field, H, at the center of a 
solenoid of finite length, Z centimeters, 


H=0AnrNI cos a, (4.1) 


where H=magnetic field intensity, in oersteds, 
at the center of the solenoid, N=number of 
turns per centimeter of length, J=current, in 
amperes, and a=the angle determined by the 
axis of the solenoid and a line intersecting the 
axis at the center of the solenoid and passing 
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through the winding at the end of the solenoid 
(tangent a =the radius divided by the half-length 
of the solenoid). The magnetic moment M in 
e.m.u. of a solenoid is given by the relation 


M=0.inIa, (4.2) 


where n=total number of turns= NL, a=cross- 
sectional area of the solenoid in square centi- 
meters (xR? for solenoid of radius R centimeters), 
and J=current, in amperes. Figure 2 shows a 
double solenoid of length Z, having an inner 
winding of radius R; and N, turns per centimeter, 
and an outer winding of radius R, and N; turns 
per centimeter. From relation (4.1) and the ge- 
ometry of Fig. 2 it can be shown that the net field, 
H, of a double solenoid is equal to the difference 
between the fields of the inner solenoid, H,, and 
the outer solenoid, Ho, 


H = H,—H2=0AnIL(Ni/(L?+4R2)! 
—N2/(L?+4R??)']. (4.3) 


From relationship (4.2) and the geometry of 
Fig. 2, it may be shown that the net moment, M, 
of the double solenoid is equal to the difference 
between the moments of the outer solenoid, M2, 
and the inner solenoid, M,, 


M= M,—M.2=0.17ILLNiR2— NR?’ |. (4.4) 


Thus for a net zero external moment M, shall 
equal Mz, i.e., the relation 


N,/N2= (R2/R:)? (4.5) 


should hold. (If the solenoids under discussion 
consist of more than one layer of wire it can be 
shown that the radius, r, of a single-layered 
solenoid which has the same total number of 
turns and the same moment as a multi-layered 
coil is 

r=((r2+roritr?)/3)}}, (4.6) 


where 7, is the outer radius and r; the inner 
radius of the multi-layered coils.) _ 

5. Several double-solenoid neutralizing coils 
have been made and tested. Shown in Table I 
are the dimensions of a typical one of these. 
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The constants of this double-solenoid neutral- 
izing coil substituted into relations (4.3) and 
(4.4) gave: 
H/I=6.5 oersted/ampere, 
M/I=0.02 e.m.u./ampere. 


The above data can be compared with the data 
calculated for an ordinary solenoid of length, 
L=2.54 cm, 2R=1.59 cm, and total number of 
turns n= 16. This gives nearly the same ratio of 
field intensity to current as the double solenoid, 
but it has an external moment about 150 times 
larger, i.e., for a single solenoid: 


H/I=6.73 oersted/ampere, 
M/I=3.17 e.m.u./ampere. 


6. Tests were made to determine the amount 
of external field associated with the double- 
solenoid neutralizing coil. With sufficient current 
in the double solenoid to produce a net internal 
field of 3.5 oersteds, the maximum external field 
along the axis was 0.006 oersteds measured one- 


‘half inch from the end of the coil. Use of a three- 


inch diameter toroidal neutralizing coil like that 
shown in Fig. 1 to produce a similar internal 
neutralizing field resulted in an external field, 
parallel to the internal field, of 0.02 oersteds 
measured at a point one-half inch to the side 
of the toroid. 


MEASUREMENT OF MAGNETIC MOMENT 


7. If a small magnet is placed in the center of 
a long solenoid, the flux linkage with the coil is 
proportional to the component of the magnetic 
moment of the magnet along the axis of the coil, 
M, in accordance with the relationship 


M=/10K (e.m.u.), (7.1) 


where J is the flux linkage in Maxwell-turns of 
the magnet with the solenoid, and K is the coil 
constant of the solenoid in oersteds per ampere.* 


TABLE I. 








Inner solenoid Outer solenoid 
(1) (2) 





Length (L) 2.54 cm 2.54 cm 
Diameter (2R) 1.59 cm 2.46 cm 
Total number of turnis (n) 24 10 








* The relationship, M=2/10K, is apparently not well 
known. It may be proved from the reciprocity relationship 
for mutual inductance. The magnet may be replaced by 
a small (secondary) coil of equal moment. The field, Hy, 
at the position of the magnet or small coil due to a current, 
I,, in the primary coil may be expressed in the form 


H,=Kp,lI> (e.m.u.), (7a) 
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This is still true if the long solenoid is replaced 
by any coil which, when energized with current, 
produces a uniform field over the region occupied 
by the magnet. An ideal type of coil for the 
measurement of magnetic moments of small 
magnets is one which produces a uniform field 
inside the coil and zero field outside. This makes 
it possible for the magnet under test to be 
removed a relatively short distance from the coil 
to a position of practically zero flux linkage; if 
this is done suddenly it is possible to measure 
the change in flux by connecting a fluxmeter to 
the coil. The double solenoid very nearly ap- ° 
proaches this type of ideal coil. 

8. The magnetic moments of small magnets 
were measured by this method using a double- 
solenoid coil 20 inches long supported by a tube 
made of non-magnetic and non-conducting ma- 
terials which held the test magnet in place. The 
tube was long enough so that the double solenoid 
could be easily slipped from over the magnet to 
a point 12 inches from the magnet where the 
linkage of flux was negligible. The change in flux 
linkage was measured with a fluxmeter, and 
values of magnetic moments were calculated 
from fluxmeter readings and the current constant 
of the double solenoid. 
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which equation defines the coil constant Ky. Suppose the 
secondary coil has m, turns and an area projected on a 
plane normal to the direction of the field of a, square 
centimeters. Then A,, the flux linkage in the secondary due 
to the field, Hy, is 


A. = Hyasn,, (7b) 
and the mutual inductance is 
Msp=),./Ip = Kyat. (7c) 


The mutual inductance is also equal to flux linking the 
primary for unit current in the secondary 


Mep=)p/I.=,0.0./ Me. (7d) 


In the last equality, the current has been expressed in 
terms of the component of the moment of the coil in the 
direction of the field: 


M,=I ah. (Je) 


Equation (7.1) may be obtained by combining (7c) and 
(7d) and expressing the coil constant K in oersteds per 
ampere instead of in e.m.u. (K =0.1K,). The moment of 
= magnet, of course, replaces the moment of the secondary 
coil. 
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The Class of Electron Lenses Which Satisfy Newton’s Image Relation 


R. G. E. Hutrer* 
Division of Electron Optics Stanford University, California 
(Received May 1, 1945) 


Newton's image relations are satisfied for all electron lenses if object and image are located 


outside the region of the electromagnetic fields. It will be shown that a s 


ial class of electron 


lenses exists for which these relations hold for any object and image position. The results will be 


compared with those obtained by W. Glaser. 





1. THE FOUR CARDINAL POINTS; 
IMAGE CONSTRUCTION 


T is known in light optics that a simple con- 
struction can be used to find the image of any 
object if a single set of the so-called “four cardinal 
points’’ of a lens is known. This construction is 
based on Newton’s image relation and the 
magnification formula. The fact that one and only 
one set of cardinal points is necessary to construct 
the image of any object is the justification for the 
introduction of the concept of the cardinal points. 
The four cardinal points are the two focal 
points zFo, zF;, and the two principal points 
ZHo, 2H;. The distances between the focal points 
and their corresponding principal points are the 
focal lengths fo, f;. The planes through the 
cardinal points perpendicular to the optical axis 
are called cardinal planes. There are two principal 
planes Hy, H;,** and two focal planes Fo, Fj. 20 
and 2; are axial object and image coordinates re- 
spectively. r(zo) and r(z,;) are points on the object 
and image respectively. r, z are considered to be 
cartesian coordinates. Focal lengths are positive 
quantities, and distances above the optical axis 
are positive, below the axis negative (see Fig. 1). 
Newton's image relation can then be written in 
the form: 


(2; —2F;)(29—2Fo) = — fofi. (1) 


The magnification formula is: 


M= r(z;)/r(Zo) = fo ‘(29 — ZF 0) = ( == 3 —2F;) / fi. (2) 





* This paper is part of a dissertation submitted to the 
department of physics and the committee on graduate 
study of Stanford University in partial fulfillment of the 
requirements for the degree of doctor of philosophy. The 
author is now a research physicist with Sylvania Electric 
_ Products, Inc., New York. The paper was completed while 
he was research associate in the Division of Electron 
Optics at Stanford University. 

** It is to be noted that in Fig. 1 the principal planes are 
crossed, that is, the object and image spaces overlap. The 
answer to the question whether all electron lenses have 
this property is not known to us. We will show later, 
however, that all lenses discussed in the present work 
must have crossed principal planes. 
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Based on these two formulae is the following 
well-known image construction (see Fig. 1). A ray 
(the first principal ray) entefing the lens from the 
object side parallel to the optical axis intersects 
H; at A. The first locus of the image point J of the 
object point O is the straight line (Azr,). The 
second ray (the second principal ray) from O 
through the focal point zr» intersects Ho at B. 
The second locus for the point I is the straight 
line through B parallel to the optical axis. J is 
the point of intersection of (Azr;) and the parallel 
line through B. 

This simple construction gives correct results 
in all cases where object and image are outside 
the region which the actual lens occupies along 
the optical axis. Within the limits of the first 
order imagery the behavior of the lens is com- 
pletely determined by a single set of cardinal 
points. 


2. THE LIMITED USEFULNESS OF THE FOUR 
CARDINAL POINTS FOR THE 
IMAGE CONSTRUCTION 


The significance of the four cardinal points is 
lost if these quantities become functions of the 
object position. We want to illustrate this for the 
case of a light optical system made up of infinitely 
many glass lenses arranged along an axis (see 
Fig. 2). : 

For any image only a finite number of lenses of 
the total system takes part in its formation. These 
lenses can be replaced by a single lens which has 
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the same focal points and focal lengths as the 
combination lens. However, the number of lenses 
changes every time the object or the image 
changes position by one lens in the whole system 
so that every time this occurs a new set of 
cardinal points must be determined. 

Any electron optical system whose axial field 
distribution is of infinite extent is analogous to 
the described glass lens arrangement. Object and 
image path which are the straight line portions of 
the total path cannot be distinguished anymore. 
It is hence impossible to speak of a one to one 
correspondence between the two sets of rays. The 
path of an electron is continuously curved. 

It is therefore a surprising fact that nevertheless 
electron lenses exist which can be described in 
terms of a single set of cardinal points. 


3. THE CONDITION FOR AN ELECTRON LENS TO 
HAVE A SINGLE SET OF CARDINAL POINTS 


In order to determine the conditions for an 
electron lens to have a single set of cardinal points 
we start out with the paraxial ray differential 
equation. The solutions of this equation namely 
have now to satisfy Newton’s image relation and 
the magnification formula, which in turn will 
impose conditions on the field quantities of the 
lens. 

The paraxial ray differential equation of a 
combined electrostatic and magnetic electron 
lens is given by: 


(2) -7’’(z) ++3¢'(2)r’ (2) 
+[i¢'"(2)+ (¢/8m)H*(z) y(z2)=0, (3) 


where $(z) is the electrostatic potential distri- 
bution along the optical axis, H(z) is the magnetic 
field strength along the same axis, r(z) is the 
function describing the distance of the electron 
from the optical axis* and e and m are the charge 
and mass of the electron respectively. The 


primed quantities are derivatives with respect 
to 3. 


* In case of the presence of a magnetic field the electron 
will describe a spiral motion about the optical axis. r, z, and 
a third coordinate—the angular displacement—would be 
necessary to describe completely the position of the electron 
everywhere in space. It is not essential in our problem, 
however, to consider the angular displacement of the 
electron. We can think of the r and z coordinates as a 
system of polar coordinates rotating with the electron 
around the z axis. 

1E. Brueche and O. Scherzer, Geometrische Elektronen- 
optile (Verlagsbuchhandlung, Julius Springer, Berlin, 1934). 

2L. M. Myers, Electron Opiics (Chapman and Hall, Ltd., 
London, 1939). 
31. G. Maloff and D. W. Epstein, Electron Optics in 
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Object 





ANAL A 


WV 


an image 


F 1G. 2. Light optical system for which focal points and focal 
lengths are functions of object and image position. 


To determine ali the electron lenses whose 
optical properties are fully described in terms of 
a single set of cardinal points means mathemati- 
cally to find all the functions ¢(z) and H(z) for 
which the solutions r=r(z) of Eq. (3) satisfy 
Eqs. (1) and (2), where 2 and 2; are any two 
successive roots of r(z). 

Equation (1) requires that a projective trans- 
formation shall hold between the successive roots 
of the path equation r(z). This restriction on the 
solutions r(z) of the differential Eq. (3) will lead 
us to a limited class of field functions ¢(z) and 
H(z). Equation (2) imposes a further restriction 
on the solutions of Eq. (3). This condition will 
put a restriction on the free choice of the focal 
lengths in case of purely magnetic or purely 
electrostatic lenses and will exclude the existence 
of combined lenses. 

A solution of the stated problem becomes 
simpler if we apply a standard transformation to 


Eqs. (1)—(3). Let: 


r(z) =y(z)- exp | —> f San 


=y(z)-p-*(z). (4) 
Then Eq. (3) becomes :5 


ywor{—— 4 a) ber=o, ¢ (5) 
8m ¢(z) 16 $(z) P 














or: 
y’’ (2) +A(z) -y(z) =0; (6) 
where: 
e H(z) 3 s¢'(z) 
Ms) =— +—( (7) 
8m (2) 16\ ¢(z) 


In case of a purely magnetic lens ¢’(z)=0, 
¢(z)= V=const., where V is the accelerating 





Television — Hill Book Company, Inc., New York, 
1938), first edition. 
*V. K. Zworykin and G. A. Morton, Television (John 
Wiley and Sons, Inc., New York, 1940). 
5A. Recknagel, AEG Forschungs- Institut, Jahrbuch. 
Sonderheft: Unleventivedion (Verlagsbuchhandlung Julius 
Springer, Berlin, 1940). 
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potential : 


e 
\(z) =An(z) = 
8m 





H(z), (8) 


r 


and in case of a purely electrostatic lens H(z) =0: 


3 /¢'(z)\? 
a) =n(2)=—( = ) : (9) 
16\ ¢(z) 


We see that if the roots of the solution r(z) 
satisfy Eq. (1) the roots of y(z) will do likewise 
since both functions have the same roots (r(z) 
differs from y(z) only by an exponential function). 

The condition (2) however becomes: 


y(z;) € fo aa - Y ¢ “4 














y(Z0) 2o—ZFo 4/2, o(u) 
0 ; 3 
a f (<*). (10) 
Zo — ZFo \o(Zo) 


Now the problem is to find all the functions 
\=A(z) defined by Eq. (7) for which the solutions 
y=y(z) of the differential Eq. (6) satisfy the 
conditions (1) and (10). 


4. DERIVATION OF A FUNCTIONAL RELATION FOR 
THE FIELD FUNCTION 


The first step towards a solution of the problem 
will lead to a functional relation which A(z) must 
satisfy. 

The general solution of the differential Eq. (6) 
can be written as: 


(11) 


where y:(z) and y2(z) are twolinearly independent 
solutions and ¢, and cz are two arbitrary con- 
stants. If y(z) has a root at z=29 we can express 
y(z) in the form: 


y(2) =C1-y1(2Z) +€2-y2(z), 








y(z) = C{yo2(Z0) -y1(2) —yi(20) -yo(z)}. (12) 
Let z=2; be the root following zo; then: 
Sofi 
y(ai)= C{ yas) -94( 20. ) 
. 20 — ZF o 
= ys(24)-y4( 2. fs )| =0, (13) 
20— 2Fo 


assuming that Eq. (1) holds for 2 and 2;. Hence: 
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fu 28) 
y1(Z0) >, Sep Zo— 2Fo 


Y¥2(Zo) ‘ ( Sofi ) 
Vol 2Fi— 
Zo — 2Fo 


If we assume (1) to be true for all values z= 2») we 


have: 
ies 
Vit 2Fi— 
yi (2) Z—2ZFo 


yx(2) ( fofi ) 
Yat SFi— 
Z—ZFo 


Differentiating both sides with respect to z, we 





: (14) 











, (15) 











get 
1 Sofi 
= ° ’ 16) 
y2"(z) Sofi \ (2—2Fo)? 
y3t( ers ) 
. 2—2Fo 
since: 
yo(z) -y1' (2) — vi (2) - yo" (z) = const. (17) 


for all values of z (Wronskian determinant). 

We obtain a similar expression for y,(z). 
Therefore y(z) must satisfy the functional 
relation : 


Z—ZFo 


° o( er —_ Sefs ) 
(fofi)* 2—2Fo 


y(z) must also satisfy the differential Eq. (6), i-e., 


( Sofi ) 
yl 2ri— + 
2—ZFo 








y(z) = (18) 


(Z—2Fo)* 


(fof)? 


(»-= 2) 
“yl 2Fi— =0. 
Z—ZFo 


We can therefore conclude that A(z) must satisfy 
the following functional relation : 
(2—2Fo)* 


(=) 
———-=\f{ 2F;-— . 
(fof i)? Z—2ZFo 


5. DERIVATION OF AN ANALYTIC EXPRESSION 
FOR THE FIELD FUNCTION 





d(z) 





(19) 


A(z) 





(20) 


In order to simplify the discussion of Eq. (20) 
we will first apply a number of transformations. 
If we let: 
(21) 


Z—ZFo=x and 2F;—2Fo=d, 
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then Eq. (20) becomes: 


Sofi oe 
Lid—- = -L(x), 22 
( ) Gare — a 























since 
A(z) =A(x+2F5) = L(x) (23) 
and 
fofi Sofi 
a( 20: ) =0( srot2r.—2r— ) 
Z2—Z2Fo bo 
=L(27.-2n—- ~~) =1(4—-™*), (24) 
x x 
Furthermore let 
=(fof;)*-t d =E; (25) 
x= (fof an ay 
then Eq. (22) becomes: 
1 
u(B-—)="-Mo, (26) 
since 
L(x) =L((fof:)*-t]= M(d), (27) 
and 
(a) =1a- nell 
x (fof :)*t. 








se. ty 1 
=H eso irae -u(z—-). (28) 


We see that a function satisfying the functional 
relation (26) is given by: 


N(t) (29) 


~ (1 —Et+e)?" 


which can: easily be verified. This N(¢t) is es- 
sentially the function representing the magnetic 
field distribution along the optical axis whose 
lens properties were discussed by W. Glaser.*® 
If M(t) is assumed to be another solution of 
Eq. (26) the ratio M(t)/N(t) defines a new 
function f(t): 
fQ™=M(b)/NO (30) 


and f(t) must satisfy the functional relation 
1 
f(z--)=103 (31) 


if f(t) are all the solutions of Eq. (31) the solutions 
* W. Glaser, Zeits. f. Physik 117, 285 (1941). 
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of Eq. (26) are then given by 


1 
M(t)=N(t)- 8 pes . 32 
MO=NO-JO=-— $0. (32) 


The fixed points of the transformation 


t'=E-1/t (33) 


are 


E /E? \} & we 
a-—+(—-1), e-—-(—-1). (34) 
2 4 2 4 ; 


Depending on the value of E there are either two 
real, one real or two conjugate complex fixed 
points. 

If: 


1. E>2, a and @ are real quantities; 
2. E=2, a=8 and real; (35) 
3. O<E<2, a and 8 are conjugate com- 

plex quantities. 


Let 
t= (8u—a)/(u-1) 
and 
t' = (Bu’ —a)/(u’—1), (36) 
or 


u=(t—a)/(¢—B) and u’=(t'—a)/(t’—8); (37) 
then 
u’ = (t’—a)/(t’—B)=p(t—a)/(t—B)=p-u, (38) 
where 

p=B/a (39) 


since t/= © for t=0, (a/8)-p=1. 
Then f(t) becomes a function F(u) which 
satisfies the functional relation : 


F(p-u) = F(u). (40) © 


If E>2, pis a real quantity and 0<p<1. For 
this condition the only solution of Eq. (40) is: 


F(u) =const. 
hence 
f(t) =const. 


and the only solution of Eq. (26) is given by 
Eq. (29) multiplied by a constant 


M(t)=const./(1— Et+#)*. (41) 


This function, however, does not lead to fields 
which can be realized in practice. 

The same is true for the case E=2; we there- 
fore can omit a detailed discussion of these two 
cases. 
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If 0<E<2, we can write E as: 


E=2-cos ¢ (42) 
where ‘ 2 
<g<4/2. 
Then a. 
a=cos g+1 sin g=e'*, 
. | (43) 
B=cos g—isin g=e-'*, 
and the value of p in Eq. (38) is given by 
p=B/a=e*", |p| =1 (44) 


u and u’ are now complex quantities. The ¢ axis 
is transformed into the unit circle and F(u) is a 
real and continuous function on the unit circle 
|u| =1. Since 


F(p-u) = F(u). (45) 
We expand F(u) in a Fourier series 
n=+0 
F(u) = F(e®)= > ane. (46) 
Then Eq. (45) is: 
n=+20 n=+O 
>. a,e'" @-2¢) _ Zz. a,e*™. (47) 
Hence ie; bate 
a,(e~?i"*¥—1) =0. (48) 


(i) If ¢/x is irrational, a, must be equal to 
zero for all n except »=0 and f(t) =const. There- 
fore N(t) given by Eq. (29) is the only solution. 

(ii) If g/x is rational, ¢/r can be written as 








b/q where p>0, g>0 and p and gq are relative 
prime integers. Then: 


e~tinepls = 1, (49) 
for all m=mq, m being an integer. Since 
u=(t—a)/(t—B) =e", 
eint = simet = ((t— a) /(t—B))™; 
and 
m= + 
F(u)=f(t)= DL cm((t—a)/(t—B))™, (50) 
with mri 
axzB=e'* and .c_.=m, 


where the c,, are any arbitrary constants. 
We therefore obtain the following complete 
solution of the functional relation (26) 


1 
M() =—__—_—__ 
(1—Et+f)? 


2 t—a\™ = t—B\™ 
XE e|(—) +(—) |, on 
m=0 t—B t—a 
which is a real function (two conjugate complex 
terms of the series (51) are combined to give a 
real expression). 
If we reverse the substitutions (21) and (25) 


the expression for \(z)—the function which 
satisfies Eq. (20)—becomes: 


(2—43(sri+2F0) —(i(2ri—2Fo)? mela pp. 


{2—43(2F;+2Fo) +4 (sri —2Fo)?— fof: _]*} 22 





@ Ce 
A(z)= 2D 
, m=() (fof)? 


We should like to remark that our distinction 
between rational and irrational values of ¢/z 
loses its meaning in physical applications and is 
mentioned only to give the complete mathe- 
matical solution of the problem. 


" 6. THE RESTRICTING INFLUENCE OF THE 
MAGNIFICATION FORMULA 


So far we have used only the lens equation as 
a restricting condition imposed on the solutions 
ef the paraxial ray equation. In terms of the 
solution of Eq. (6) we had written the magnifica- 
tion formula as: 


(zi) /y(20) = Lfo/(%0—2Fo) ](o(2:)/(z0))* (10) 


or 


674. 


—. ($2) 


[(s—2Fo)? —(zF;—2Fo)(2—2Fo) + fof: ]"**? 








2—2Fo (O(2Fi— fofi/(2—2Fo))\— 
y(z) = ( ) 


$(2) 


-y( =" Safi ). (53) 
2—2Z2Fo 


If we compare this equation with Eq. (18) which 
we derived by using the lens equation only 


0 





(2) Z—2Fo Sofi 
ee fefy »( +.) 


we see that 








(18) 


o( i—fofi/(2- 0) t 
(fofi)* =fo wont << any, (54) 
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or 





= [See (55) 


fo $(2) 


This expression must be independent of z, 
otherwise the pair of focal lengths would become 
a function of the object position. 

We will distinguish three cases: (1) Purely 
magnetic lens. ¢(z) =const in this case. Therefore 
fo=fi. (2) Purely electrostatic lens. H(z)=0 in 
this case. It is: 


A(z) = de(z) = v6 (9' (2) /$(2))?, 


and ¢’(z)/¢(z) must satisfy the following func- 
tional relation: 


¢'(2F,— fofi/(2—2Fo)) _ (—2Fo)? ¢’ (z) 
o(2F:— fofi/(2—2Fo)) fofi (2) 
or 
¢' (zri—fofi/(2—2Fo)) fof ¢’(z) 
= : =——. (56) 
o(2Fi—fofi/(2—2Fo)) (2—2Fo)? (2) 


Integration with respect to z gives: 











6( =r. Sofi ) =const. #(2), (57) 


Z2—2Z2Fo 


Hence fi/fo is independent of z. (3) Combined 
electrostatic and magnetic lens. 

The requirement that f;/fo must be constant 
leads to Eq. (56) by the differentiation of Eq. (55). 
This means that the term ;:(¢’(z)/(z))? in 
Eq. (7) must be of the form given by Eq. (52). 
But it is then impossible to satisfy the condition 
that the total expression for (z) is of the same 
form since ¢(z) appears in the denominator of the 
term (e/8m)(H?(z)/o(z)). A combined lens satis- 
fying the lens equation and the magnification 
formula cannot therefore exist. 


7. PHYSICAL SIGNIFICANCE OF THE CONDITION 
FOR THE PARAMETER E 


Solutions of the problem exist only for the case 
that 0<E<2 where E was defined by 


d/(fofi)t=E <2. (25) 
For this condition: 


d/4 <fofi. (58) 


It can be seen from Fig. 1 that fo+f;<d is the 
condition that the principal planes are not 
crossed. This inequality contradicts, however, 
Eg. (58) as we can readily prove by the following 
consideration : 
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It is fo+fi<d, hence 











SPF ISATIE 
4 4’ 
or 
a Sofi So’ fi Sofi Sofi 
‘é:.2:°4°6 ¢ = 
d? (fo—fi)? 
| FN aah tae 
since 


(fo—fi)? d? 
———=0, —Z=fofi, 
4 4 


- but Eq. (58) required d?/4<fof;. Therefore all the 


lenses given by the function A(z) have crossed 
principal planes. 

Within the range 0<E<2 it is possible to 
choose E arbitrarily. Since: 


d 2F;—2Fo 


(fof)® (fof)! 


d, 2Fo, and fo-f; can be chosen so that O0<E<2. 
The ratio fo/f; and thus the individual focal 
lengths are then determined by Eq. (55). The 
field function of the lens having these lens 
parameters 2Fo, 2F:, fo, f; is then determined by 
Eq. (52). 





8. SIMPLE EXAMPLES OF ELECTRON LENSES OF 
THE TYPE DISCUSSED 


It remains to be shown that certain of the class 
of field functions given by Eq. (52) are physically 
realizable. We will list a number of simple field 
functions constituting proper lenses which belong 
to this class. 


A. Magnetic Lenses 
(1) H(z)=1/(i+2*), (3) H(z)=2/(1+2")!, 
(2) H(z)=1/(i+2)', (4) HA(z)=2/(1+2*)’, 
(S) H(z) =(2*+1)*/(1+27)*. 


Sketches of some of these field functions are 
shown in Figs. 3 to 5. 


B. Electrostatic Lenses 
(1) (z)=exp [c arctan z], 


1 
(2) (s)=exp Carr 
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Fic. 3. Sketch of H(z) =1/(1+32*). 


f Hu 








Fic. 4. Sketch of H(z) =1/(1+:°)!. 


Hey 








Fic. 5. Sketch of H(z) =2/(1+2*)!. 





(3) o(z)=k-exp lee ah 


1 
(4) (s)=k-exp (« ). 
1+2? 


Some of these axial potential distributions are 
shown in Figs. 6 and 7. 





9..COMPARISON OF THE RESULTS OBTAINED 
HERE WITH THOSE FOUND BY 
W. GLASER’ 


The investigations presented here were con- 
cluded when a paper published by W. Glaser in 
1941 came to the attention of the author. In this 
paper Glaser obtained a solution of the same 
problem by a different mathematical approach. 


7™W. Glaser and E. Lammel, Ann. d. Physik [5] 40, 
367 (1941). 
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The form of his solution also has a different ap- 
pearance from the one given here so that it is 
worthwhile comparing the two results. Glaser’s 
solution is as follows: 

An image equation of the type of geometrical 
light optics: 


2, = (az+5) /(cz+d) (59) 


exists between object and image coordinates z 
and 2, if and only if the electromagnetic field 
along the axis has the general form 


3 so’'(z)\? +e A(z) 
a oe Yam (2) 

1 (a—d)?+4be 
ve [c2?+(d—a)z—b |? 











[i—J(s)], (60) 


where J(¢) is a periodic function of the argument 


2cz—a+d 




















¢=arc cot : (61) 
[—L(a—d)*+4bc]]}! 
with the period 
a+d 
+ =arc cot : (62) 
[—[(a—d)*+4bc}} 
We can write Eq. (1) also in the form 
ZF 2— it2PF2F 
25> Lief od (63) 
2—ZFo 
Comparison with Eq. (59) yields: 
a=2F; b=—[fofiterzro], c=1, 
(64) 
=-—2ZFo0, 
and 
E = (zr;—2Fo)/(fofi)'=(a+d)/(ad—bc)'; (65) 
furthermore 
22—(zFi+22Fo)+2Fo 
cot ¢= 
[—(eri—ar0)*—4fof.} 
t—E/2 
= , (66) 
[ —(E?/4—1) }! 
2F;—2ZFo 
cot y= 
[ —(sr;—2Fo)?—4 fof; }! 
E/2 
= (67) 





~ [-(E2/4-1)}" 
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E?/4-1 1 
A(t) = 
fofi (P—-Et+i1} 
t—E/2 
| 1-a(are cot amd} (68) 
(1—E?/4)! 
where J has the period: 
E/2 
- = arc cot ——————_-. (69) 
(1—E?/4)! 
It was: 
t= (Bu—a)/(u—1), (36) 
with 
E E? ; E E? 4 
a=—+(—-1) | e->-(—-1) ; (34) 
2 4 2 4 
hence 
t—E/2 u+1 
————_ aj, (70) 
(1—E?/4)! u—1 
therefore 
t—E/2 
3) arc cot al 
(1—E?/4)! 
u+1 
= (are cot (-i )). (71) 
u—1 
since 
1 %«@ 1-—i%; 
arc cot -=-— In - 
z 2 1+7%2 
( u+1 
1+ 
u+i 1 u—1 
Jar cot ( i ))-s —In 
u—1 2 u—1 
u+1) 
1 
=I(< u), (72) 
hence 
1 
f=-Inu. (73) 
From Eq. (42) 
. 2 cos g=E, (42) 
cos* g=E?*/4, sin? g=1—E?/4, 
cot g=cos g/sin g=(E/2)/(1—E?/4)! 
hence 
g=y- (74) 
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ee ee we ewe ee ee wee ee eee eee ee es 


“>> ,ose 








: 
Fic. 7. Sketch of $(z) =k-exp [3] . 
Therefore with Eq. (44) 
p= * (44) 
or 
y= (7/2) In p. (75) 


J(g) is, according to Glaser, a periodic function 
with the period 7; it can therefore be expanded 
in a Fourier series 


+00 +a 
> a,e*ety = > ae; (76) 
this means 
yn 
e**y¥=1=exp | -2=:(=-)] (77) 
x 2 
hence 
n=m-2q. (78) 
This leads to: 
J(S)=L damge™? “=P domgu™?= M(t). (79) 


Both results of the same problem are therefore 
identical despite their difference in formal 
appearance. 
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Rigorous Treatment of the Electrostatic Immersion Lens Whose Axial 


Potential Distribution Is Given by: ¢(z) =¢.e* *"* 


R. G. E. Hutrer* 
Division of Electron Optics, Stanford University, California 
(Received May 9, 1945) . 


Investigations have shown that the electrostatic immersion lens, given by the axial 
mero distribution $(z) =o exp(X arctan 2), is the simplest electrostatic lens for which 

ewton’s image equations are satisfied throughout the field. The general solution of the 
paraxial ray differential equation is derived. Based on this solution and the fact that Newton's 
image equations are satisfied, exact expressions are derived for all important optical quantities 
such as the focal lengths, the location of the four cardinal points, magnification, the object- 
image relation, and the ——— and chromatic aberration for any object position. The 
optimum positions of an object are determined so that the chromatic and spherical aberra- 
tions are reduced to a minimum. All quantities are represented graphically. A plot of the 


equipotential lines in space is given for a special value of the voltage ratio. 


LECTRON lenses of rotational symmetry 

are uniquely described by their axial field 
distribution. The field everywhere in space, the 
electron paths, as well as all optical quantities, 
can be expressed in terms of the axial field 
function and its derivatives. Special conditions 
imposed on lenses, such as minimum of aberra- 
tions or the validity of Newton’s image equa- 
tions, can be interpreted as conditions for the 
axial field function. We therefore have fre- 
quently, at the start of the actual investigation 
of an electrostatic lens, the axial potential as 
the basic quantity. 

In electrostatic electron optics it often occurs, 
however, that lenses given in terms of their axial 
potential function are hard to realize in practice 
by means of simple electrode configurations. 
Singular points of the potential cause this quan- 
tity to be multivalued in the normal cartesian 
plane, and it does not help the physicist that 
the function becomes single-valued in the Rie- 
mannian plane. No mechanic can build lenses in 
this kind of space. In such cases it is only 
possible to approximate to the desired potential 


* This paper is part of a dissertation submitted to the 
department of physics and the committee on graduate 
study of Stanford University in partial fulfillment of 
the requirements for the degree of doctor of philosophy. 
The author is now a research physicist with Sylvania 
Electric Products Inc., New York. The paper was com- 

leted while he was research associate in the Division of 
lectron Optics at Stanford University.- 
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distribution. The problem whether or not the 
axial potential function 


(2) = do eX arctan (s/a) 


(1) 


will give rise to a multivalued potential distribu- 
tion in space can be solved by a computation of 
the equipotential surfaces. By inspection of such 
an equipotential plot it is then also possible to 
determine an electrode configuration producing 
exactly or approximately the axial potential dis- 
tribution of Eq. (1). This distribution is plotted 
in Fig. 1 for some values of the parameters ¢o, « 
and a. It is compared with that for a special 
cylinder lens. 


1. THE EQUIPOTENTIAL PLOT 


A first glance at the function ¢(z+ir) reveals 
that it has a logarithmic singularity at r=1, 
z=0, since 


itu i 2+i(1+7) 
arctan “=— In —— =— In ——————_ 
t—-u 2 —s+2(1-—7) 
if 
u=z+ir. 


This shows that the series 
a) (—1)° r 2v 
o(8, 7) = ama (Z) 
= ©)? (2) , 


commonly used for the determination of the 





JOURNAL OF APPLIED PHYSICS 














-i - 





FOP POSES S @ OO RENE Ee HF RHEE EOOES © EE EEEE HOSE ESOS SEES EEO ERE EEE etesseee 


S?antan = 
——-—-— Dias 2445-e = 


A $7 arc tan £053 
——— -Gm-2008- ) 





4 . 4 rn ‘ nm 
“L8 -i6 -I4 -12 -10 -8 -& = 


a a a 


P Fic. 1. Axial potential functions. 


potential distribution in space, will converge 
only for values of r with r < ro, where 


ro= (1+27)!. 


For values of r near the limits of convergence a 
great many terms are necessary before it is 
permissible to break off the series for numerical 
computation. If r>vro the series is divergent, so 
that no information can be gained from it about 
the character of the potential function. 

Another method makes use of the solution of 
Laplace’s equation in the form of a definite 
integral: 





1 2r 
¢(z, r) =— f o(z+ir cos a)da 
2x Jo . a 
=— f Re\¢$(z+%r cos a) }da. 
T /0 
In our special case we have: 


‘ , Ki 12(1+7rcosa)+2 
K arctan (2+7r cos a) =— In \ ) 





Ps ’ 
i(1—r cos a)+z 
so that: 

eX arctan (z+ir cos a) 


“| 2?+(1+r7 cos a)?\! 
=exp > In ( ) 
2 2?+(1—r cos a)? 








—1 arctan 
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22 | 
(1—r? cos? a) —2? 4) 


The real part of this function is then: 


Re[eX arctan (z+ir cos @)'] 


K 2z 
= }exp | — arctan 
2 (1—r? cos? a)z? 
K_  2#+(1+r7 cos a)? 
Xcos | — In 
4 2+(1—r cos a)? 








; (2) 
or 


Re[eX arctan (z+ir cos «)] 


K 2 
= {exp [> (arctan > RRR. 


1—rcosa 





Z 
+arctan =...) | 
1+7 cosa 
K_ 2#2+(1+r cos a)? 
Xcos | — In , 
: 4 2#+(1—r cos a)? 
since: 
a+b 
arctan a+arctan b=arctan ———, 
1—ab 
2 
arctan ——————+- arctan ——————- 
1—rcosa 1+7 cos a 
22 
=arctan 





(i—?r?* cos? a) —2? 


The integral form of the potential function 
permits the application of the process of analytic 
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continuation, by which the potential distribution 
in the region inaccessible by the series can be 
obtained. The analytic continuation of ¢(z, r) 
into the region not covered by the series can 
readily be accomplished by properly observing 
the multivaluedness of the arctan-function in 
the exponent of the exponential function: 


K z 
exp (> Te ———— 
2 1—rcosa 








of Eq. (2). If z is positive and r>(1+32")! then: 


z 





arctan = 2x—arctan " 
1—rcosa |1—r cos a| 
and for negative values of z and r>(1+2?)!: 


zZ 
arctan = —2x-+arctan 
1—rcosa 








—?rcos a 


‘If s=0 it follows that for r>1: 


1+r7 cos a 





2 arc cos l/r K 
(+0, n== f gt els cos | = 
Fo 


K_ 1+rcosa 
In feo+ a cos =H a, (3) 
2 —itrese oe en We 1—rcosa 


1+r7 cos a 





2 are cos l/r K 
¢(—0, r)=- f e~**/2 cos E In 
0 


us 


This shows that the potential function is uni- 
valued in the plane the physicist is interested in, 
except for the portion of the r axis for values 
from r=1 to r= ~. It is through the cut along 
this line that the function g= ¢(z,r) disappears 
into the other sheets of the Riemannian plane. 
In the cartesian plane ¢(z,r) is double valued 
along this line. 

For the application of formulae (3) and (4), 
which determine the potential along the r axis, 
it is helpful to know the values of a for which 


E 1+7 cos a 
cos | — In 
2 +1+r cosa 





(5) 


becomes indefinite, has maxima and minima and 
zeros. For any z these values are determined by: 
K_ (1+7 cos a)?+2? w 
—In =n-—. 
4 (1—rcos a)?+2 2 


n=2m+1, m=0, 1, 2,--- 





zeros; 


n=2m ,m=0,1,2,--+ maxima and minima. 


Equation (5) leads to: 


| th —+ ( th? 2+1 
cos = —_ —_ — - 
1,2 co co ( )) 


For the case z=0 the zeros are therefore given by: 


cosh 


(2m+1)x 
—— +1 
K 





COS @2>=— 


h (2m+1)x 
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fiat ~{" cos Cage RON a. (4) 
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The integrand becomes indetermined for: 


a@-=arc COS -, 
r 


T 
for r=1, 0=a, <>: 
The general behavior of the integrand is illus- 
trated in Fig. 2. It is noteworthy that the zeros 
of the integrand appear to be clustered around 
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Fic. 2. The integrands in Eq. (3). 
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the indeterminate point (for z=0) as given by: 
cos a, = 1/r, 


so that the areas under the oscillatory part of 
the curve could very well be neglected when 
performing the integration by means of the 
planimeter. 

After we had obtained a general picture of 
the behavior of the equipotential lines in space 
another method was suggested to us by Professor 
W. W. Hansen of Stanford University. Since this 


method reduced the amount of computational ~ 


work considerably, it was used to determine the 
actual potential plot of Fig. 3. 

The solution of Laplace’s equation in spherical 
coordinates with rotational symmetry about the 
z axis can be written as: 


gi(z, p) =>. anp"P,,(cos 6) 


n=0 


for p <1, O=0=7; 


* 1 
g2(2, p)= 2) bx——-P, (COS 8) +94 
n=) p 
for p >1, 0S0=7/2; 


= 1 

¢3(2, p) =) Cn—-P,(cos 0) +o_« 
n=0 p 

for p >1, x/2S0=7; 


where p is the distance from the origin of the 
spherical system of coordinates, @ the angle 
between the positive z axis and p, and P, is 
the mth Legendre polynomial. (Fig. 4.) 

The coefficients a,, b,, ¢, must be determined 
so that the potential functions ¢, ge, gs; represent 
the axial potential 
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Fic. 3. Equipotential lines. 
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Fic. 4. Spherical coordinates (with rotational symmetry 
about the z axis). 


$(z) = @X arctan z (6) 


in their respective regions. 

The coefficients a, will be obtained by an 
expansion of our axial potential in a power series 
of z. We get: 


(n+1) !dngi = Ka,n!—n(n—1)ax_1(n—1)! 
w=Q0,1,2,...3 ae=l1. 
The first 10 coefficients (m!a,) are given by: 
Q=1, a=K, 2a,.=K*-2, 3!a;=K*—8K, 
4!a,= K*—20K*+24, 5!a;=K'—40K*+184K, 
6!ag= K*—70K*+- 784K? —720, 
7 !a;= K’—112K*°+-2464K*—8448K, 
8!a,= K*— 168K*+-6384K*—52352K?+40320, 
9 !ag= K®— 240K" + 14448K5 
— 249760K*+648576K. 


In order to determine the coefficients 0, we 
expand 


1 
o(z) =exp |x; — arctan -) | © >s=1 


b4 


in a series of powers of 1/z. We get for the coeffi- 
cients in: 


1 1 1 
$(z) = 40 tb0-—+bi-—+b2-—-=+ eee 
Z z 3 





the following general relation to the coefficients 
re. 


bn =(—1) "day -exp (<>). 


Furthermore: 


T 
Pix =exp (x7) = Ve. 


Similarly we obtain the coefficients c, by an 
expansion of: 


1 
$(z) =exp |x( —5-+aretan -) | 
2 


—x<s<-l, 


where 


i.e., 


oo 


1 1 1 
$(z) = 0 tlo*— Fi +2 — + eee 
Pid 2° z , 


ghtl 


+02: 





The relation between the a, and the c, is found 
to be: 


= _ r 
Cn =An41 EXP (-x*), o_.=exp (-x) =). 


For the special voltage ratio V2/ V,=6, K is equal 
to 0.571. The coefficients a, then assume the 
following values: 


ao = 1, a5= 0.0836, ay = 0.0521 ’ 

a,=0.571, ag=0.0775, ay, = 0.02821, 
a,= 0.163, a;= 0.0535, ay = 0.04419, 
a3=0.1593, ag=0.0618, a;3=0.02192, 
a4=0.10437, ag=0.0377, ay=0.03735. 


We did not attempt to prove generally that 
the solutions join on the surface of the sphere of 
radius 1. The actual calculation showed within 
the limits of computations with finite parts of 
infinite series that the solutions assume the same 
values on the respective parts of the sphere. 

All series converge, though the convergence is 
better, of course, for values of r sufficiently far 
away from the value 1. 


2. THE ELECTRODE CONFIGURATION 


Due to the fact that the potential function is 
multivalued along the axis z=0 and r>1 it 
appears to be impossible to produce the axial 
potential distribution (6) by means of a simple 
electrode configuration. There is, however, a 
possibility of realizing the exact distribution in 
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either one of the following two ways: (a) by means 
of an electrode configuration consisting of re- 
sistance material, (b) by the use of metallic 
rings arranged along the optical axis and charged 
to proper potentials. In both cases there is an 
infinite number of possibilities for the choice of 
the geometry of the arrangement, owing to the 
fact that any surface of rotational symmetry 
about the z axis can be chosen as basis. The axis 
along which the potential is double-valued must 
be excluded. According to method (a) one has 
to construct an electrode of resistive material 
whose inner surface is the chosen surface of 
rotation. If the circular ends of this surface are 
then connected to a potential source a current 
will flow along the walls of the surface. The 
driving potential and the resistance distribution 
can be chosen in such a way that a potential 
distribution is obtained which can be determined 
from the equipotential plot.: The method (b) 
also uses a surface of rotation as a basis, sub- 
dividing it into a number of rings whose radii 
are equal to the radius of the surface at the 
coordinate z. Each metallic ring is connected to 
a tap on a potential divider whose potential 
value can again be determined from the equi- 
potential plot. In this manner the axial potential 
distribution could be achieved to any degree of 
accuracy. Figure 5 shows the potential distribu- 
tion along a cylindrical surface of radius r=1, 
as well as the axial potential distribution. 


3. THE SOLUTION OF THE PARAXIAL 
RAY EQUATION 


While for most fields the solution of the par- 
axial can be achieved by approximate methods 
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Fic. 5. Axial potential distribution and potential variation 
along a cylinder of radius r=1. 


11. G. Maloff and D. W. Epstein, Electron Optics in 
Television (McGraw-Hill Book Company, Inc., New York, 
1938), first edition. 
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only, its solution for the special function: 
4 2 
(2) = do exp (5.3 arctan *) (7) 
a 


can be given in general form in terms of ele- 
mentary functions. It is: 





. 
’ 


(e) 4 5b 1 (; ab “) 1 
2) =—V3koo- exp [ —v3k arctan— } —— 
. 3 oe" wd 3 — a/ 1+(z/a)? 


3 /9'(z)\? = (k/a)? , 
(aa), 
16 \ $(z) 1+(z/a)? 
The paraxial ray differential equation for the 
case of a purely electrostatic lens is:'!~ 


hence: 





o(2)r’’(z) +30’ (z)r’(2) +40" (z)r(z)=0; (9) 
and with: 
1 ¢ ¢'(u) _ 
r(z) =y(z)-exp (-; f av) =30)-6 ‘(z), 
4J $(u) 


Eq. (9) becomes: 
3 oe) 
x —{ —— )=0. 10 
y w+, ba J? (10) 


Substituting (8) in Eq. (10) gives: 


9 


(k/a)? 
[1+ (c/a)? 





y'"(2)+ (2) =0. (11) 


Let: 
w= (1+). 


It has been shown by W. Glaser that the solution 
of (11) can be written in the form: 





y(z) = (1+ (<) )| C, sin (« arccot *) 


z 
+C, cos (« arccot *) |, 


a 
and because of: 


¢*(z) =exp ( — 4v3k arctan *) A 


a 


z\?\3 z 
r(z)= (1+ (<) ) exp (- 4v3k arctan *) 
a a 


-{C, sin [w arccot (z/a) ] 
+ C2 cos [w arecot (z/a)]}. 


This is the general solution of the paraxial ray 
equation describing the path of an electron 
through the lens whose axial potential is given 
by Eq. (7). The two constants of integration 
C, and C2 (or K, and Kz in the second form of 
the solution) are available to satisfy special con- 
ditions such as r(zi)=n, r’(a:)=n' or r(a) =n, 
r(%)=r2 for a given lens characterized by the 
two constants p and a. 


It is useful to write the path equation in 
terms of the variable ¢; where: 


cot g=2/a. 
Hence: 
tan ((7/2)— ¢)=cot g=2/a 
or 
arctan (z/a) =(x/2)— 9, 


and, therefore: 


R(o) =Ki exp| 3e(o—") 


sin w(¢+Ke2) 


sin 9 





(12) 


4. THE DETERMINATION OF THE CARDINAL POINTS; THE EXPRESSIONS 
FOR THE FOCAL LENGTHS 


In order to determine the location of the focal points, the constant K» in Eq. (12) is chosen so 
that the path remains finite at z=+ © or g=0. It is then: 


r\ |sin we 
R(¢)=K;, exp| 3#( o—7) (13) 





2 sin ¢ 


* E. Brueche and O. Scherzer, Geometrische Elektronenoptik (Verlagsbuchhandlung Julius Springer, Berlin, 1934). 
: 7 M. Myers, Electron Optics (Chapman and Hall, Ltd., London, 1939). 
aaa 


K. Zworykin and G. A. Morton, Television (John Wiley and Sons, Inc., New York, 1940). 
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for Kz must be equal to zero. The slope of this path is equal to zero at z= +0: 


dr dRd in “CO —sin -cOos 
eg feeramernar en) 





de és sin? ¢ 








sin wy . T sin? 9 
+— -4y3k exp| 1v3e( o—7) |} (- ). 
sin ¢ 2 a 


=0.. 


Hence: 
dr(z) 


dz 





| sa 


That is, the electron enters the lens parallel to the optical axis. Its distance from the axis is at that 
point: 
Tr 
r(+ 0©)= Kw exp (=n). ‘ (14) 


For a ray entering the lens at z= — © parallel to the axis the corresponding expressions are: 


i | 93*(e—*) iste) oe 


sin ¢ dz 





=0, — 22) =—Ku exp (1385). (15) 


z=—@ 
According to the definition of the focal points their location on the optical axis is given by: 
sin [wero ]=0 


from Eq. (13). The index (0) shall indicate the object side of the lens, i.e., the side of the lens with 
negative abscissae. The other side is called the image side and is characterized by the index (i). 
The focal points on the object side are then: 


woerg=n-e, n=1,2,3,--- or gro=n-x/w, n=1, 2, 3,-- 


n=1 is the first focal point through which the ray coming from z= + © or g=0 passes. 
The focal points on the image side are determined by: 


sin [w(gr;— 7) ]=0, 
hence 
w(gr;—m)=—nx; gFi=x—nr/w; n=1,2,3,--- 
or in terms of the axial coordinate z: 


2Fo=a cot (mx/w), 2F;= —a cot (nx/w). 


Knowing the location of the focal points, it is now possible to determine the corresponding focal 
lengths. These focal lengths are defined by: 








r(+ «) ae 
= = focal length on the object side, 
r’ (2Po) 
r(— #) ' 
fi=— = focal length on the image side. 
r' (zF;) 


Using Eqs. (14) and (15) and the following expressions: 


1 
r’ (2Fo) = K, sin - exp | x34( ==) |-o-(- 1)*-L, 
w w 2 a 


. T x 1 
r’ (z¥s) = K, sin (x27) exp | 34(r-n=—") |-0-(— y*.(-), 
w w 2 a 
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fo and f; become: 


tg 
a exp (- {N3kn™ ) 
w 








fo=(-—1)"" ’ 
ae 
sin n— 
Ww 
(16) 
ns T 
a exp (+3:3en") 
Ww 
fi=(-1)" 
a T 
sin n— 
Ww 


There will be a value of k for which fo/a as- 
sumes a minimum. This value of Rk will, how- 
ever, not give a minimum for f;/a. In case of 
the corresponding magnetic lens investigated 
by W. Glaser® the two focal lengths are equal . Tres 
and assume their minimum value therefore for 
the same value of the lens parameter. 

The condition for the minimum for f/a is: 











| . 





Fic. 6. Graphs for the determination of the minimum 
values of the focal lengths. 


[d(fo/a) |/dk =0, 


which yields the transcendental equation: 


tan x/(1+k?)!=v3k. (17) 


For the other focal length f;/a the condition: 
(d(f,/a) ]/dk =0, 
yields: 


tan 4/(1+k?)!= —v3k. (18) 


The two values of k satisfying the two Eqs. (17) and (18) can be obtained graphically. These 
values are: 


ke;=1.185 or (V2/Vi)s;=5404, 


k~=2.178 or (V2/V1)fo27,275,000 


(see Fig. 6). . 
If the voltage ratio V2/V, increases, the focal length on the image side assumes its minimum 
before the focal length on the object side does. These minimum values are: 


fo/a= 0.199 
fi/a= —4.45. 


The values of V2/ V; are very high so that it seems almost impossible to obtain the minimum focal 
y nig po 


~ lengths in practice. 
The focal lengths and the location of the focal points as a function of the voltage ratio are plotted 


in Figs. 7 and 8. 
° W. Glaser, Zeits. f. Physik 117, 2859(1941). 
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Fig. 7 — The focal Lengths 


t te, ia 


+10 ——— 




















eo 


4 10 eo’ 10000 ~"% wal 


Fic. 7. The focal lengths. 























5. COMPARISON WITH DESIGN FORMULAE FOR THIN LENSES 


Due to the fact that it is impossible to obtain accurate formulae for the focal lengths of electron 
lenses approximate formulae are often used for design purposes. It is therefore worth while to com- 
pare these formulae with the accurate ones we obtained for the special lens discussed here. 

These approximate formulae for the focal lengths are: 


a 7 [¢ (s)F 

















1 
—_— z, (19) 
fo S[b(+2)PJ_.  ¥(z) 
ee +# [9'(z) 
nee f [¢’(z)] Qs (20) 
e —8[¢(—2) PJ.  gl(z) 
With ‘ 
2 
$(z) =exp | 53 arctan (=) |-e [=a (5- e) | 
3 a 2 
ayes! E os (:)] a E u(* )] rm 
=—V3k- exp | —V sctan § — 3 )-————— *-VS- Vv _—— ‘sin’ ¢, 
"e 3 . 3 7 a 1+(z/a)? ee . 3 2 . antes 
a 
dz= ———_d¢ 
sin? ¢ 
the expressions (19) and (20) take the form: 
2(k?+3 1 ; —2(k?+3 1 
fo_2(k*+3) fe_—2(+3)_ sd 


= V3k ‘exp (2v3km)—1' a v3k ‘1—exp (—3V3kr) 


The focal lengths are represented graphically in Fig. 7. It can be seen from these graphs that for 
a fairly wide range of voltage ratios the formulae for thin lenses give a relatively good approximation. 
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This is a remarkable fact considering that the comparison between the corresponding formulae for 
the magnetic lens proved to be very unfavorable for the approximate expressions. 
Another expression in the theory of thin lenses relates to the ratio of the focal lengths. It is: 


; oo) 
fet ” 


From the approximate expressions (21) this ratio is indeed given as the square root of the ratio of 
the asymptotic voltages: 


f =exp (3v3kr). 


0 


Our accurate formulae on the other hand yield: 


fi 7 
—=exp | se | 
So (1+k?)! 


Both expressions are approximately equal for small values of k, as can be seen from the equations 
or from Fig. 9, where the ratio of the focal lengths is plotted as a function of the voltage ratio. 
6. NEWTON’S IMAGE RELATIONS 


In the preceding paper® we have shown that Newton’s image relations hold for the special lens in- 
vestigated here. This relation, however, can be derived independently of these considerations. 
It is known that object and image coordinates for any lens are related by: 


r1(Zo) - r2(;) — 71 (24) - r2(Zo) =0, (23) 


where 7;(z) and r2(z) are any two linearly independent solutions of the paraxial ray equation.§ 


ye 








Fig 8 = The location of the focal points. 
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Fic. 8. The location of the focal points. 
*R. G. E. Hutter, J. App. Phys. 16, 670 (1945). 
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Fic. 9. The ratio of the two focal lengths. 
We take 
. | r\ ] sinw¢ ; 
ri(y) = Ki exp | $v3k{ ¢—— } |-——., (24) 
4 2/ 1 sin ¢ 
2 j r\ |] cosw¢e 
ree) = Kz exp | §v3k{ e¢—— } |: ——— (25) 
2/1] sing 








as such solutions. Substituting (24) and (25) in Eq. (23), this equation becomes: 


sin [w(¢go— ¢i) ]=0, 
hence: 
Yo— gi=n- r/w. , (26) 
It was: : 
cot g=2/a. 
Equation (26) then becomes: 
2% 2,/acotnr/w—1 


, ’ 
a 2;/a+cot nx/w 


Zo ri 2; 7 7 1 
|= —cot n= | +c0t n= | -|1 +cot? n | tah science, 


a wiLa w w sin? nx/w 





or: 


Equations (16) show that: 


a? 
jeh= ~ gin? nr /w 
It follows therefore: 
(29 —2Fo)(2:—2F;) =fo-fi (27) 


which is Newton’s image relation. 
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In order to derive the magnification formula a general formula for the magnification of any elec- 


tron lens is used as the starting point. If r(z) is any particular solution of the paraxial ray equation, 
the magnification can be expressed by:°® 


M=r(2;)/r(z). (28) 
If we take the particular solution given by Eq. (25) the expression (28) becomes: 


cos (w¢y;) sin go 














cos —n sin 
M=exp [4\3k(g:— 0) ]-— weap (— Whale eee 
sin (wi) cos (w¢po) sin (go—mm) cos (w¢p) 
exp (— 3vV3kn2/w) 1 fo 1 
=(—1)" : poe (= 84 et 
cot nmr/w—cot go sin nx/w a 2Fo/a—%/a 
hence: f 
2;—-2F; 
M=———= ' (29) 
2—2F =f 


The simple geometrical image construction based on the two expressions (27). and (29) is illus- 
trated for a special choice of the voltage ratio in Fig. 10. The two principal rays used have nothing 
in common with the actual electron path except the end points. Within the limits of the Gaussian 


dioptrics (i.e., neglecting aberrations) the influence of the field can be replaced by the two principal 
planes and the focal points. 


The principal planes are those planes with unit magnification. With M=1 in Eq. (24) these 
planes have the following axial coordinates: 


(—1)"-' exp (— 4v3kn2/w)+cos (n2/w) 
sin (nr/w) 
(—1)" exp (+ 3v3kn2/w) —cos (nx/w) 


sin (n/w) 


2Ho = fot2Fo =a 
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2H; =f;+2F,;=a 
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Fic. 10. Image construction for the special lens with the axial potential distribution ¢(z)e* s°%* *, 
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Fig 11 The location of the principal pounts 
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Fic. 11. The location of the principal points. 





If n=1 " 
exp (— 4v3k-2/w)+cos (x/w) 
L£Hop =a . Sigg ’ 
sin (2/w) 
exp (4vV3k-2/w) —cos (4/w) 
2h j= es, 
sin (2/w) 
where 
w= (1+k?)!. 


ZHo, 2H; are plotted in Fig. 11 as a function of the voltage ratio V2/ Vij. 


7. CHROMATIC ABERRATION 


If a beam of electrons entering the field of an electron lens is not homogeneous with regard to 
speed the curvatures of the individual paths will be different. This means that electrons leaving one 
point will not again go through a single point after having passed the focusing field. This effect is 
called chromatic aberration. 

Expressions for the magnitude of this effect have been derived by W. Glaser’ and A. Recknagel.® 
Glaser’s expression can be written as: 


1 zi wey r2(z) 
be r= %a-AV——_— P dz. 30 
[o(zo) }! J. ( o(z)/ [¢(z)}! . (30) 


where ¢(z) is the sum of the initial accelerating potential of the beam and the potential up to the 
coordinate z in the lens. AV is the voltage inhomogeneity occurring in the electron beam satisfying 





7W. Glaser, Zeits. f. Physik 116, 56 (1940). 
$A. Recknagel, AEG Forschungsinstitut, Jahrbuch. Sonderheft: Uebermikroskop (Verlagsbuchhandlung Julius 
Springer, Berlin, 1940). 
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the condition: AV&¢(z). r=r(z) is a solution of Eq. (9) which satisfies the following conditions: 
r(zo) =r(z;)=0; r’(zo) =1. (31) 


a is the angular aperture of the electron beam at the object position. 
For the case of the electrostatic lens discussed here the expression can be evaluated. We obtain: 


[5*(»-5) 
73 Nk-a “i t 6 


‘4 7kh2+3° 


The dimensionless quantity: 


[s3*( 5) | 
>. 4 —-V au a= 
Go 32: ut UO 


a 47k+3 


chr = @ 








sin? gp 


4 © 
{1—exp| -3.6k (32) 











{1 inne 
exp 3 scam 


characterizes the property of the lens as regards chromatic aberration. C.»;/a is a function of k as 


well as of the object position go. We will show that an optimum position ¢o opt exists so that C.,;/a 
is a minimum. From: 


sin? gp 


d(Cenr/a 
(Cenr/ - 
d go 


’ 


it follows that: 


4 
rhe sin @p = 2 COs gp; 
or 


20 opt 





cot Yo opt — 2v3k. 





-2 Fig. 12 The variation of the chromatic aberration as a 


function of the voltage ratio % : 
(abject located at the optimum position ) 





















































Fic. 12. The variation of the chromatic aberration as a function of the voltage ratio V2/V; 


(object located at the optimum position). 
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{& Fig. 13 The vanalion of the chromatic aberration asa function of the valtage raho  . 
(object located at the focal point of the object side ) 
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Fic. 13. The variation of the chromatic aberration as a function of the voltage ratio V2/Vi 
(object located at the focal point of the object side). 


It is interesting to note that this point is the same as that for which the axial potential distribu- 
tion, has its point of inflection. For such a point the second derivative must be equal to zero, i.e., 


4 Zo 
d? exp (<3 arctan ~) 


a 





4 
=0, or —vV3k—2—=0. 
d(z/a)? 3 a 


20 infl 


Hence: 





= $v3k. 


a 


The value of the constant for the chromatic aberration at this optimum position is given by: 
Cenr opt V3. 3+4k 4 4 2 
ional I: —exp (-33 v4 )| exp( ——v3k arctan “Gk. 
a 4 3+7k? 3 1+k? 3 3 


Cenr opt/@ is plotted in Fig. 12 as a function of V2/ V3. 
For high magnifications the object is near the focal point 279. Substituting the value for go= 2/w 
in Eq. (32), the chromatic aberration constant becomes for this object position: 


(1+?)*—2 
exp | —3v3kr 
Cer focal 3 v3k (1 +k?) 4 T 
=- . {1—exp ——V3k——— |}. 
a 43+7k? sin? r/(1+k?)! 3 (1+k?)! 
Figure 13 shows a plot of the constant for this last case. Figure 14 finally gives a graphical repre- 
sentation of the dependence of (29 op:/a) on V2/ Vi. 


The dependence of the chromatic aberration constant on the position of the object is shown in 
Fig. 15 for the special case V2/Vi=6. 
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8. SPHERICAL ABERRATION 


It has been pointed out before that perfect imagery could only be obtained if the rays through 
an electron lens would satisfy a differential equation of the type of the paraxial ray equation.® 
Since this is the case approximately only for the so-called ‘‘paraxial” rays the imagery of real elec- 
tron lenses is never perfect. If the angular aperture is increased the image of a point becomes a 
small disk whose radius increases with the third power of the angular aperture. This radius divided 
by the magnification M serves as a measure of the aperture defect termed ‘‘spherical aberration.’” 


Glaser and Scherzer" have derived expressions for: 
5.» =Ar/M 
in terms of the field function along the optical axis: 
o=9(2), H=H(sz) 
and the paraxial ray path function: 
r=r(z). 
For a purely electrostatic lens Scherzer’s expression becomes: 


— [io F 5% 149% 7 3 # 
Pn Sw 
16[ (zo) } Jz 


dsp = 


4 “ae 3S é6r 2 # 
The function r(z) again satisfies the conditions (31): 


r(Zo) =r(z,)=0; r’(2o)=1 


—¢//34.— —+-— — —~—- ¢/2— r'dz. 





Fig. 14 The variation of the optimum pesition for minimum chromatic aberration 
as a function of the voltage ratio } ' 
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Fic. 14. The variation of the optimum position for minimum chromatic aberration as a function 


of the voltage ratio V2/ Vi. 


*L. Marton and R. G. E. Hutter, Proc. I.R.E. 32, 3 (1944); Proc. I.R.E. 32, 546 (1944). 


‘© H. Busch and E. Brueche, Beitraege zur Elektronenoptik (Johann Ambrosius Barth, Leipzig, 1937). 
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Fic. 15. The variation of the chromatic aberration as a 
function of the object position (for the voltage ratio V2/ V1). 


Making use of the fact that the function r(z) vanishes at the end points of the integral, Eq. (33) 
can be transformed into: 


5 piece [ler ee ridz (34) 
160 6(z0) }# Jeo 2 @¢ 2¢ 4 
by repeatedly applying the process of partial integration and substitution from Eq. (9). 
For the case of the electrostatic lens discussed here we can evaluate the expression (34). 


Let: 
K =(4/3)./3k. (35) 


We then obtain the following formula: 


He (-2)-bn (ade 
sp — ex a a = exp a> sin am 
> Sidef sin’ os Pl on , Quo ik 


Kr T 
—Cexp (-=) cos 2( 7) +B sin 2¢9+C cos 2enf (36) 


w Ww 





Again as in the case of chromatic aberration a constant can be defined which is characteristic 
of the particular field function. C,,/a is a function of the lens parameter & as well as of the object 
position go. This constant is defined by: 


Cw K? Kr Kr\ . ™ : 
= A| 1—exp { —— } |+ B] exp | ——} sin 2{ go—— } —sin 2¢ 
2w 2w w 
Kr 7 
+C] exp (-=) cos 2( =) —cos 20 
2w w 











@ 512w*sin* gp 


» (37) 
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where the coefficients A, B, C of Eqs. (36) and (37) are: 


1024 K (64+13K?)’ 
-2|——*- 131K4+512K?—768 
K?+16 169K*+1312K?+2304)’ 
cox{ 3K*— 64 n 39K*—316K?—2560 f 
K?+16 K?(K?+4) 169K*+1312K?+2304 
Equation (37) can also be written as: 


c. K? 
°= ——{L+M sin 2¢+N cos 2¢0}, (38) 
a 512‘ sin* gp 














where the coefficients L, M, and N are related to the A, B, and C by the following equations: 


: Kr 
L=A|1-—exp (-=)| 
a 2w 


Kr 2x ! Kr\ . 2x 
iis ‘ied (-=) cos ——1|+C exp -—*) sin —, 



































2w w 7 w w 
Kr 2a ‘ Kr\ . 2x 
N=C\ exp | —— } cos ——-—1]—Bexp { —— } sin —. 
2w w : 2w w 
Just as was the case for the chromatic aberration constant there exists an optimum position ¢, or 
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Fic. 16. The variation of the optimum position for minimum spherical aberration as a function of the voltage ratio V2/ V. 
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The variation of the Spherical aberration as a function of the valtage ratio 
(object located at the optimum position ) 
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Fic. 17. The variation of the spherical aberration as a function of the voltage ratio V2/ Vi 
(object located at the optimum position). 


Xm, SO that the spherical aberration constant be- 
comes a minimum. From the condition: 


d(C,p/a)/dxo=0, 


for the minimum of the constant C,,/a it follows 
that the position coordinate of this minimum is 
determined by the cubic equation: 


Xm? (2L+- M)+3Mxn?+2Lxm+ M =0. 


The solution of this equation is plotted in Fig. 
16 as a function of the voltage ratio V2/V;. 

Figure 17 shows the variation of the spherical 
aberration constant under the condition that 
the object is always located at the optimum 
position. 

Figure 18 represents the spherical aberration 
as a function of the voltage ratio if the object 
is located at the focal point of the object side of 
the lens. 

In Fig. 19 finally the spherical aberration con- 
stant is plotted as a function of the object posi- 
tion for the special voltage ratio of V2/Vi=6. 
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9. A GENERALIZATION 


In order to obtain a greater degree of freedom 
in matching axial potential distributions of simple 
lens configurations to the one treated here the- 
oretically, a slight generalization is helpful. We 
shall investigate a lens which is defined by: 


4 z 
$1(2) =o exp (<5 arctan =), 2=0 


ay 
(38a) 
4 z 
$2(Z) = 02 exp (534. arctan =). 2=0 


a2 


The formulae for the optical quantities of this 
lens will only be slightly more complicated than 
those derived previously for the simple lens. 
The constants oq: @;, ki, and @o2 de, ke must be 
determined so that: 


¢1(0) = ¢2(0), 
1’ (0) = ¢2'(0), 
i.e., the potential ¢ and the field strength along 


the optical axis are continuous functions of the 
axial coordinate z. 


(39) 
and 
(40) 
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Equation (39) requires that: 


oi = oz 
and from Eq. (40) it follows that: 


k;/ay =k»/de. 


If only one expression of the form (1) is used 
to describe a potential distribution which has 
two asymptotic potentials V2 and V; the con- 
stants k and ¢o are completely determined. If, 
on the other hand, two expressions of this form 
are chosen, k; and ke are only determined after 
a convenient choice of @9 has been made. a; and 
dz, are again related to the slope of the axial 
potential distribution at z=0 or at any other 
point. 

A comparative plot using various funcions is 
given in Fig. 20. 

The relation between the object and image 
coordinates is obtained in the following manner: 
An electron path intersecting the optical axis at 
the object coordinate go and whose slope at this 
point is: 

ry’ (zo) = 1, 
can be expressed as: 


ay 
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Fig. 19 | The variation of the spherical aberralion 
as @ function of the cbyct position 


(for the vellage ratio $+ 6) 
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Fic. 19. The variation of the spherical aberration as a 


function of the object position (for the voltage ratio 
V2/Vi=6). 


The ordinate of this path at z=0 is then: 


ri \20= — > ———— exp | avses(5— «)| 
@1 SIN go 


L(+) 


, 





ie 
sin — 


and its slope at the same point is: 
Y 1 7 

ry’ | 20 =—— exp | 4V3kil —— go 
@1 SIN Go 2 


| 5h sin 1 (=- «) +o cos w(=- «) | 


An electron path intersecting the optical axis at 
the image coordinate ¢; is given by: 


= Cexp| indhs(e-7) | sin [we( e— e)] 


| sin ¢ 





The constant can now be determined so that 
both paths have the same ordinate and the same 
slope at z=0, if go and ¢; are really object and 
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Fic. 20. Axial potentials of the generalized lens. 


image coordinates. The matching conditions are: 
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Division of the second equation by the first gives: 


} v3 k2+w2 cot | »(=— ) | 





a2 us 
=— 4v3ki +o: cot | »(=— «) ||: 
a 2 
or since: 
a2 
ke =—k, 
a 
[-(5— ee) Jc [(F-) | 
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This is the relation between the object and 
image coordinates. For the case a;=a_2 and 
ki=k2, it becomes identical with Eq. (26) as 
required. 

Using Eq. (38) the matching constant C can 
be determined; it is: 


Tv 
~~ | nan (F- ») | 





C=— - 
@1@2 SIN go 
' | ae?” sin? | »(=— »)| 
| 2 
~ 4 
+2*w,? cos? | ox(=— «) |} ‘ 


W. Glaser® has given a graphical method for 
the determination of the image coordinate if the 
object coordinate is given. This construction is 
illustrated in Fig. 21. The triangle AOC has the 
base diw2, and the angle at A is made equal to 
wi(7/2—¢o). If, then, OB is equal to daw, the 
angle at B will be equal to we(r/2—¢,). The 
image coordinate is then given by: 


T 
2;=d2 arc cot EG e) | 


The magnification can be determined by the 
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same method which we used previously for the 
single function lens. A solution of the paraxial 
ray equation for the left-hand field which is 
different from zero at the object position and 
which is zero at z=0 is given by: 





in —1/2 
ry, =exp [3v3ki(e—2/2)]-— te */ 4 
¢ 


The path function in the second part of the field 
is a continuous extension of that in the first part 
if the constants in: 

sin [we( ~—C2) ] 


r2=C, exp [$vV3k2( e— 2/2) ]- : 
sin ¢ 


are: 





Ci =d201/dy2, C2= a/2 


The magnification is then nothing else but the 
following ratio: 


M= 1 (2) /r2(2;). 


Substitution of the particular solutions r; and r 
yields: 


ko w 
M=— — exp [4v3{he(o:—#/2) —ki( go —27/2)}] 


1 We 


sin [we(gi—2/2) ] sin go 
sin [wi(go—2/2)] sin 9 





The two focal points can be determined by 
means of Eq. (38). If 9,=0, ie., z=+o and 


go= 7, i.e., z= —« they are given by the two 
following equations: 
F T . Ry we aa 
tan n(5- or) | = > a tan ord 
tan wn(=- en) = o.. = tan Phat 
‘4 2 d Ry we ies 














We should like to emphasize that the notion 
of focal length loses its meaning for the lens 
described by the two expressions: 
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Fic. 21. Construction of the image coordinate. 


4 Zz 
$1(Z) = do exp (<3, arctan =) . 


ay 


4 z 
$2(2) = do exp (<3. arctan — 


a2 


since the latter does not belong to the class of 
lenses which satisfy Newton’s image relations. 
In this case the simple image construction can- 
not therefore be used. 

Spherical and chromatic aberration can also 
be computed for the lens given by Eqs. (38a). 
This can be done by integrating the expressions 
for C.nr and C,, in two steps. The first integral 
has the limits go and 2/2, the second integral is 
taken between 2/2 and g;. The influence of 
different ratios k2/k; on the magnitude of the 
aberrations could be investigated. This investi- 
gation, however, would be beyond the scope of 
this work and shall therefore be omitted. 
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The Methods of Specirying the Properties of Viscoelastic Materials 
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Seven methods (Voigt model, Maxwell model, operator equation, mechanical impedance 
function, creep curve, relaxation curve, and dynamic modulus function) of specifying visco- 
elastic behavior are discussed. A number of exact relations between these methods of speci- 
fication are worked out in detail. The majority of these relations are simple enough to be of 
practical value although a few are too cumbersome. Approximate relationships between the 
creep curve, the relaxation curve, Maxwell model, and Voigt model are discussed; and nu- 
merical examples show the magnitude of errors introduced by the approximation to be small 
even in quite unfavorable cases. A consideration of the practical utility and physical meaning 
of the various methods of specification distinguishes between (1) those of general descriptive 
value and those of direct experimental value; (2) those useful in a phenomenological study of 
mechanical behavior and those more suited to a formulation of molecular theory. A summary 





of the present molecular theories is presented together with their interpretation in terms of 


the Voigt and Maxwell specifications. 





INTRODUCTION 


URING the last decade the investigation of 

the mechanical behavior of viscoelastic 
materials has greatly expanded. Accompanying 
this development there has been a corresponding 
increase in the number of methods used by 
various investigators to summarize their data 
and to specify the properties of the materials 
studied. In many cases it has not been clearly 
recognized that the methods used were actually 
equivalent and this has led to an incomplete 
utilization of the data available. No less than 
seven distinct methods have been employed to 
define the properties (in shear) of viscoelastic 
materials. It is the purpose of this paper to 
examine these methods and show that they are 
all mathematically equivalent and that they are 
all interconnected by definite relationships. 

We are only concerned with incompressible, 
isotropic, viscoelastic materials. The term incom- 
pressible indicates that the material is much 
more easily deformed in shear than in com- 
pression; thus, the viscoelastic deformation is 
essentially a constant volume process. The term 
isotropic requires that in the unstrained state 
the properties of the material are the same in all 
directions. The term viscoelastic is here restricted 
to designate those materials whose response to 
stress obeys the superposition principle;' that is, 
_materials which can be represented by a model 
constructed from elements which obey Hooke’s 
elastic law and elements which obey Newton’s 
viscosity law. The mechanical properties of such 
an incompressible, isotropic, viscoelastic material 





tJ. C. Maxwell, Phil. Mag. 35, 134 (1868). 
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are completely defined by its response to homo- 
geneous shearing stresses. A rigorous method is 
available for calculating the behavior in the pres- 
ence of combined and inhomogeneous stresses.” 
The mathematical theory of non-linear strain 
has been developed by Murnaghan.* 

The seven distinct methods of specifying the 
properties of viscoelastic materials may be 
divided into two classes. In Class I may be listed 
the four that are more general and fundamental 
in nature. 


A. The Voigt model, consisting of a set of Voigt ele- 
ments (retarded elastic elements) coupled in series. 

B. The Maxwell model, consisting of a set of Maxwell 
elements coupled in parallel. 

C. The operator equation, Ps=2Q«, where P and Q 
are linear differential operators. 

D. The mechanical impedance function, Z(w). 


In Class II may be listed the obvious experi- 
mental curves which serve to “map out’’ the 
viscoelastic character of the material. 


E. The creep curve, consisting of the strain as a function 
of time at constant stress. 

F. The relaxation curve, consisting of the stress as a 
function of time at constant strain. 

G. The dynamic modulus curve, consisting of the elastic 
modulus as a function of frequency. 


The Class I methods are more satisfactory as 
specifications of properties since they are of the 
nature of differential equations which can be 
solved for a wide variety of transient conditions. 
Class II methods on the other hand represent 
special solutions of the fundamental differential 


2 T. Alfrey, Quar. App. Math. 2, 113 (1944). 
3F. D. Murnaghan, Am. J. Math. 59, 235 (1937). 
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equations for very particular cases. Conse- 
quently, these latter functions are not so well 
adapted to the treatment of a wide range of 
transient conditions. 

It is obviously important to be able to trans- 
form readily from one fundamental mode of 
description to another, to be able to derive the 
basic properties of a material (Class I) from 
experimental data (Class II), to be able to 
predict behavior in special conditions (Class I1) 
from the properties (Class I), and to be able to 
predict behavior in one special set of conditions 
from the known behavior under other conditions. 
If Class II methods are accepted as specifying 
the properties of a viscoelastic material, then all 
of these problems are concerned with the trans- 
formation from one of the seven modes of de- 
scription to another. 

Some of these transformations have been 
worked out and are of a form simple enough to 
be of practical value. In Fig. 3 these transforma- 
tions are indicated and labeled as a, b, ---, m. 
Several other transformations have been formally 
developed but are perhaps too cumbersome to 
serve as practical working schemes. These are 
indicated in Fig. 4 and labeled n, 0, ---, ¢. 
Finally, several transformations can be made by 
simple approximation methods. These are also 
summarized in Fig. 4 and labeled u, v, w, and x. 

This paper is divided into four parts: 


I. Brief, individual discussions of the seven methods of 
specifying viscoelastic properties. 

II. Development of exact relationships among the 
seven methods. 

III. Development of approximate relationships among 
the methods. 

IV. Discussion of the relative utility and physical 
meaning of the different methods. 


I. METHODS OF SPECIFYING VISCOELASTIC 
PROPERTIES 


A. The Voigt Model 


A Voigt element‘ consists of a spring (repre- 
senting Hooke’s law behavior) and a dashpot 
(representing Newton's viscosity law behavior) 
in parallel. This element expresses the simplest 
form of retarded elastic behavior. Such a single 
Voigt element is characterized by two param- 
eters—the modulus of the spring, G, and the 
viscosity of the dashpot, ». The strain, 2e, 
associated with the element follows the differ- 


*W. Voigt, Lehrbuch der Krystallphysik (Leipzig, 1928). 
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Fic."1. The Voigt 


model. 


rs 


$=2G; 4, +27; é 


ential equation. 
s=2Ge+2n¢€, (1) 


where s is the stress. The factor 2 in this and 
following equations arises from the fact that our 
€ represents a shearing component of the strain 
tensor, and is thus equal to 3 of the corresponding 
“classical’’ shearing strain. It follows that s 
represents the corresponding component of the 
deviatoric stress tensor.* In the case of simple 
shear, s is the shearing stress and 2¢ is the 
shearing strain. 

An important quantity is 7, the retardation 
time of the element, which is equal to 7/G. 
Another useful quantity is the compliance, J, 
which is simply the reciprocal of the modulus G. 
It will be convenient to use G and J inter- 
changeably as a measure of the equilibrium 
elastic properties. 

Figure 1 shows a Voigt model consisting of 
four Voigt elements, one of which is degenerated 
into a spring with no dashpot and one into a 
dashpot with no spring. The material is thus 
defined by the values of six parameters, Ji, J, 
nz, Js, n3, and m4. The total strain, e, is made up 
of the contributions: from the four elements: 
€=€:+6€2+6¢3+¢,. The retarded elastic elements 
obey the differential equations: 


s= 2G2€2+2n2€2, (2a) 
s= 2G3€3+ 2n36€3. (2b) 


* The stress tensor o;;, can always be resolved into two 
parts—the mean normal stress ¢, and the deviatoric stress 
tensor $;;. The mean normal stress is equal to $ (o11+o22 
+33) and s;; is given by oi;—4;;, where 6 is the Kroneker 
delta. In the case of all isotropic materials this is in fact 
a very useful resolution since it permits the isolation of 
pressure-volume effects and shear effects. The deviatoric 
stress tensor includes all shearing stresses and only these. 
In an incompressible isotropic material the role of the 
deviatoric stress tensor is even more important since no 
strain results from the mean normal stress. 
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For the degenerate elements these of course 
reduce to 

s= 2G, (3a) 
(3b) 


More generally, the behavior of a material 
which is represented by m retarded elastic ele- 
ments is governed by a set of n-differential 
equations of the form 


s= 24€4- 


$= 2G ye; + 2ni€;, (4) 


(of which one or two may be degenerate) plus 


n 


the relation e=}>- ¢;. In the limit, the set of 
discrete Voigt elements may be replaced by a 
continuous distribution of elastic compliance as 
a function of retardation time. That is, the finite 
set of physical constants—J;, mm, Jo, m2, °*°; 
Jn» In—is now replaced by a continuous function 
J(r)dr which tells how much elastic compliance 
has the retardation time, 7, in the range of dr. 
The total retarded strain is given by the integral, 


e(t) = f “&(t, 2dr, 


where ¢* is governed by the differential equation: 


2e*(t, 7) 


2r de*(t, 7) 
t)= 
s(t) Tn) 


J(r) at 





(5) 


B. The Maxwell Model 


Figure 2 shows a Maxwell model! consisting of 
four Maxwell elements, one of which is de- 
generated to a spring and one to a dashpot. In 
this model the strains associated with all elements 
are equal, but the total stress is divided among 
the elements.'* The stress on the ith element is 
related to the strain by the differential equation: 


2¢ ve 
em——§-7—"3 5 
Gin 


(6) 


which takes on simpler form for the two de- 
generated elements. As in the previous method, 
the behavior of the material is governed by a set 
of first-order differential equations. 

In the limit, this behavior goes over into a 
continuous distribution’ of elastic moduli as a 


5K. Bennewitz and H. Rotger, Physik. Zeits. 40, 416 
(1939). 


* W. Holzmuller and E. Jenckel, Zeits. f. physik. Chemie 
A186, 359 (1940). 


7™W. Kuhn, Zeits. f. physik. Chemie B42, 1 (1939). 
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function of relaxation time. The finite set of 
physical constants Gi, 71:--Gn, 7, is replaced by 
a function G(r)dr, telling the amount of elastic 
modulus with the relaxation time 7 in the range 
dr. The total stress is given by the integral: 


= fst t)dr, 


where s*(t, r) follows the differential equation , 
1 ds*(t, 7) s*(t, 7) 
G(r) at 7G(r) 





2e(t) = (7) 


C. The Operator Equation Ps—2Q¢ 


Method A involved a set of m differential 
equations of the first order, each relating a par- 
ticular strain contribution with the total stress. 
Method B involved a set of m differential equa- 
tions of the first order, each relating a fraction 
of the stress with the total strain. Method C*§® 
involves a single differential equation of the mth 


order, relating the total stress with the total 
strain: 


Ps=2Q« (8) 
where 
Pp o™ + o"—! 4 ie 
Gem Tape es 
oe" g"-l 





Q pt Cr Le + hc +o. 


ot"—1 


The parameters /;, gi, fix the properties of the 
viscoelastic material. It will be demonstrated in 
Part II that the same number of independent 
parameters are required to characterize a ma- 
terial, regardless of which method (A, B, or C) 
is used. 

In a material which is characterized by a con- 
tinuous distribution of elastic compliance J(r) 





; 
= 


2é= 


HW 





| 

= ..+ > ° 
G; sj Ni sj 
Fic. 2. The Maxwell model. 


8’ T. Alfrey, Mechanical Pro 
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instead of a model with a finite number of dis- 
crete elements, the operators P and Q likewise 


must take the forms of integrals rather than 
sums. 


oe or 
P= [ 9)—ar, 
0 ot” 


@ ot 
Q= { a)—ax. 


The properties are now fixed by the continuous 


functions p(x) and q(x) instead of a set of con- 
stants pi, qi. — ' 


D. The Mechanical Impedance Function 


The analogy between electrical networks and 
the mechanical models designed to describe 
viscoelastic behavior seems quite obvious once 
it is recognized. Voltage corresponds to stress, 
resistance to viscosity, current to the time 
derivative of strain and the reactance to —G/w. 
Consequently, the steady-state response of a 
viscoelastic material to a _ sinusoidal stress 
s=socoswt can be related to a (complex) 
mechanical impedance Z. Z = a + ib, where a is 
the ‘‘viscous resistance,” 7, and b the “elastic 
reactance,” —G/w. Both a and b depend in mag- 
nitude upon the frequency. The complex me- 
chanical impedance function Z(w) fixes the 
properties of a viscoelastic material in a manner 
completely equivalent to those of A, B, and C. 
Although the most obvious use of the impedance 
function is in connection with steady-state 
response to sinusoidal stress, the Fourier and 
Laplace transformations make it as general as 
A, B, and C in its definition of viscoelastic 
properties. 


E. The Creep Function 


If a relaxed specimen of a viscoelastic material 
is suddenly subjected (at the time t=0) to a 
stress s which thereafter is held constant, the 
strain will increase with time according to some 
function 2¢(t). Since the strain at any time is 
proportional to the stress s, dividing through by 
S gives a creep function ¢(t) which is independent 
of the stress. That is, 


(t) = 2e(t)/s. 


A given creep function is consistent with only 
one model of either type and with only one 
operator specification. Although the creep func- 
tion can be said to define the viscoelastic proper- 
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ties of a material, as has been pointed out, it is 
less satisfactory than A, B, C, and D for such a 
specification. 


F. The Relaxation Function 


If a sample is quickly forced to take on the 
strain ¢e, and is then held fixed in such a fashion 
that the strain must remain constant, the stress 
will decay with time according to some function 
s(t).!° Since the stress at any time is propor- 
tional to the strain, dividing through by the 
strain gives a relaxation function y(t) which is 
independent of the strain. That is 


¥(t) =s(t)/2e. 


The relaxation function of a material defines 
its viscoelastic properties in precisely the same 
sense as does the creep function. 


G. The Dynamic Modulus 


When subjected to a stress which varies 
sinusoidally with time (s=s9 cos wt), a visco- 
elastic material will respond" with a strain 2e(¢) 
which (in the steady state) will also vary sinu- 


- soidally with time: 2e= 2¢€9 cos (wt —a). The ratio, 


So/2€o, can be called the- dynamic modulus, G(y), 
and will in general be a function of the frequency. 
A given dynamic modulus function will be con- 
sistent with only one model of either type, and 
with only one operator specification. Hence, 
G(v) can be said to specify the viscoelastic proper- 
ties of a material. 


Il. EXACT RELATIONSHIPS AMONG 
THE METHODS 


While any one of the above methods is in 
principle sufficient to define the properties of a 
given viscoelastic material, it is desirable to 
understand how the various methods are related 
to one another. This is particularly true at the 
moment. The wide differences in choice of 
method to be found in the current literature often 
make it difficult to see the precise connection 
between the results of two different investigators 
in the same field. This lack of agreement on the 
language to be used is a healthy, if confusing, 
index of the present status of the theory of vis- 
coelastic behavior; it makes imperative, however, 


9 J. J. Thomson, Applications of Dynamics to Physics and 
Chemistry (London, 1888). 
10 E. Wiechert, Ann. Chim. Phys. 50, 335 (1893). 


1 A. P. Alexandrov and Yu. S. Lazurkin, Acta Physico- 
chimica 12, 647 (1940). 
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Fic. 3. Diagram of practical relationships. 


a ready means of translation from one actively 
used language to another. 

In this section, a number of exact relationships 
will be discussed. These may be divided into two 
classes—those which are simple enough to be 
used as working methods, and those which are 
chiefly formal in nature. The diagram in Fig. 3 
shows the translations which can be made 
exactly and yet fairly simply; in Fig. 4 the formal 
relationships, more difficult to utilize, are in- 
dicated together with approximate relationships. 
The lettering in these figures will correspond to 
that in this section. 


A. Voigt Model to Creep Function 


If the viscoelastic properties of a material have 
been specified by a Voigt model, it is always easy 
to compute the response of the material, 2¢(¢), 
to any simple stress sequence s(t), since the 
response of each Voigt element can be calculated 
separately, and then the individual contributions 
to the strain can be added together. The creep 
curve is a particularly simple case. A constant 
stress s, imposed upon a single relaxed Voigt 
element at time zero, causes a strain response, 
2¢(t), which must obey the differential equation 


2Ge+2né=s, (9) 


where s=constant. The general solution of this 
equation is 


2e(t) = Ae~'/*+(s/G). (10) 


Inserting the initial condition that «=0 at time 
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t=0 for a sample containing no residual strains 
we can evaluate the arbitrary constant, A, and 
obtain the particular solution 


det) == [1 et J=sJEL et") (11) 


If the model consists of m retarded elastic ele- 
ments, each with its own compliance and 
retardation time, the total response of the model 
is given by 


2e(t) =s > J{1—e~"**]. (12) 


If in addition to complete Voigt elements the 
model includes an unretarded compliance J; and 
a true viscous flow governed by the viscosity 7, 
the response becomes 


2e(t) =f n+ J (1—e—**) +1] (13) 
n 


i=2 


Finally, if the series m Voigt elements is replaced 
by a continuous distribution of elastic com- 
pliance, J(r), over a range of r, Eq. (13) becomes 


* 1 
26) =4 It f Hay(1—eydr+-t] (14) 
0 n 


Each of these functions can be changed to the 
standard creep function ¢(¢) simply by cancelling 
out the constant s from the right-hand side. 
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B. Maxwell Model to Relaxation Function 


While the Voigt specification allows a simple 
prediction of the strain as a function of time for 
a known s(t), the Maxwell model is designed to 
allow an equally simple prediction of the stress 
as a function of time when a known strain 
sequence, 2e(¢), is forced upon the material. 
s(t) for the known 2¢(¢) can be computed for 
each Maxwell element separately and the results 
added together. The relaxation of stress at 
constant strain is a particularly simple case. (It 
is here assumed that the strain 2¢ is suddenly 
imposed upon a relaxed sample at t=0.) The 
stress in a given Maxwell element must follow 
the differential equation: 


2€=(1/G)s+(1/n)s. (15) 


This leads to the explicit relation: 
S(t) =sye~9/* = 2Gee"!", (16) 


For a set of m Maxwell elements in parallel and 
a spring of modulus G, with no dashpot for 
relaxation, Eq. (16) becomes 


()=24 G4 E Ge (17) 


t=2 


We need not consider here the case where an 
isolated dashpot is present. Finally, if the set of 
n Maxwell elements is replaced by a continuous 
distribution of relaxation times, G(r), the stress 
relaxation at constant strain follows the relation 


s(t) =2f- G(r)e—""dr. (18) 


C. Voigt Model to Dynamic Modulus 
Function 


The steady-state relationship between sinu- 
soidal stress and strain is unique in one respect. 
Since this relationship can be reduced to a 
vectorial relation, between the two members of 
a rotating rigid pair of vectors, either s(¢) or 
2e(t) can be taken as the “known” or “im- 
pressed’”’ function, and the relative amplitude 
and phase of the other computed. If the visco- 
elastic properties are specified by a Voigt model, 
we obviously must choose s(#) as the “known’’ 
function. S(t) = So(cos wt). A single Voigt element 
must respond to this stress according to the dif- 
ferential equation 


So COS wt = 2Ge+2ne. (19) 
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The steady-state solution is 


=] 
sin (w- tan-! —) 
WT 


G(w*t?+ 1)! 





2e=So (20) 


A set of m Voigt elements, plus an unretarded 
compliance J; and a viscous flow term 7, will 
therefore have a steady-state response equal to 


—1 
sin (ot—tan- —) 
WT;i 


G(w?r?+ 1)! 





2€=50| Ji cos wt+)>- 


1 
+—sinwt}. (21) 
nw 


The extension to a continuous distribution of 
relaxation times is obvious. 


D. Maxwell Model to Dynamic Modulus 
Function 


When a Maxwell specification is used, we 
must choose the reverse formulation of the 
problem, and let 2e=2¢9 sin wt. By obvious steps 
we obtain, in the case of a set of » Maxwell 
elements, 


2ed 5 w*r*Gi ; 
s=2eq >) ——— sin w 
1+?7;? 


$5 = COS atl (22) 


and in the case of a continuous distribution of 
elastic relaxation times. 


s=A sin wt+B cos at, (23) 
where 
© 7°G(r) 1G(") 
A =2ew? f dr, B= dew fa 
Jo 1+w?r 1+w?r? 


The dynamic modulus function is then given by 
G(w) = (A?+B?)!. 


E, F. Maxwell Model to and from Voigt Model 


The translations from a Voigt model to a 
Maxwell model, and vice versa, have been indi- 
cated in Fig. 3 as being relatively simple. This is 
only partly true. If we are given a Voigt model, 
and must construct the equivalent Maxwell 
model, the first step is to find the correct struc- 
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Fic. 5. Equivalence of Voigt and Maxwell models. 


tural form of the unknown model, and the second 


is to compute the numerical values of the un- 
known parameters. The first, or topological, step 
is controlled by rules which can be formulated 
simply, but the second step can be easily made 
only in very simple cases. 

The following rules govern the topological rela- 
tionships between the two kinds of model: 


1. The Maxwell model will contain the same number of 
springs and the same number of dashpots as the equivalent 
Voigt model. 

2. A standard Voigt model (no incomplete elements) 
corresponds to a ‘‘doubly degenerate’’ Maxwell model (one 
isolated spring and one isolated dashpot). 

3. A standard Maxwell model corresponds to a “doubly 
degenerate” Voigt model. = 

4. A singly degenerate Voigt model with one isolated 
(unretarded) spring corresponds to a singly degenerate 
Maxwell model with one isolated (non-relaxing) spring. 

5. A singly degenerate Voigt model with one isolated 
dashpot corresponds to a singly degenerate Maxwell model 
with one isolated dashpot. 

6. The rather involved general rule which embraces 2, 
3, 4, and 5 is the following: The presence of an isolated 
element of one kind in a model corresponds to the absence 
of an isolated element of the opposite kind in the conjugate 
model (e.g., the presence of an isolated spring in a Voigt 
model implies the absence of an isolated dashpot in the 
équivalent Maxwell model). 


Figure 5 shows a number of specific illustra- 
tions of these generalizations. 

The numerical computation of the parameters 
‘of one model from those of its equivalent is 
much more involved. It is, indeed, out of the 
question to write down general functional rela- 
tionships for one set in terms of the other. The 
complexity which arises can be estimated from 
the following explicit equations for the relatively 


706 


simple case of two Maxwell elements in parallel, 
corresponding to one unretarded spring, one 
Voigt element, and one dashpot in the other 
language. Let Gi”, m™”, Go”, and 2™” be the 
parameters of the Maxwell model, and G,’, G2’, 
2", ns” the parameter of the equivalent Voigt 
model. The transformation equations are: 








Gi" =G."+G2™, (24a) 
v= GG" (mi +02")?(Gi"+G2") i 
ile (m™G_¥ — no@G,") = (24b) 
Ton ni™ ne” (Gi! +G2™)?(n1% + 2”) 
aie (m1%Go" — no“Gy") ee (24c) 
ns* = m+n, (24d) 


1 (GY)? (G,")? 

Gy“= (6 atG, PO 8 <<), (25a) 
28 no" ns" 
=} (Gy)? (G,")? 

G2 = — (Giva—Gira— aie —— vif -), (25b) 
23 ne! ns' 


1 
m™ -;,(4 ‘at+G,',8 


(G V)\2 G V)2 1 
eS) (=). 050 
ne" ns" a+B 
—1 


aq"! = (ae —G,",8 


(Gi)? (G,")? 1 
a aeemenene wae <<) (—) » (25d) 
n2' ns" a+, 








where ? ; 
($4545) 26a) 
a= , a 
2ns’ 2n2" 2ne" ' 
G2"G," 
Singh aimee, (26b) 
no ns’ 


Because of this complexity, it is suggested that 
numerical computations of this sort be made in 
two steps, via the operator formulation, rather 
than directly. 


G. Operator Equation to Voigt Model 


The first step in translating an operator equa- 
tion Ps=2Q¢ into the equivalent Voigt model is 
to establish the correct structural form (topology) 
of the model. The following rules govern the 
form of the model. 


1. The total number of ultimate elements of the model 
is equal to the total number of parameters involved in the 
operators P and Q. 
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2. If the Voigt model contains no isolated springs or 
dashpots, the order of P is one less than that of Q. 
(n=m+1.) The “standard” operator equation which is 
equivalent to a set of m complete Voigt elements can there- 
fore be written: 

Pett bmg +--+iee On 
pre Prati pos = an Qo€. 

3. If the coefficient g, vanishes, then one Voigt element 
of the equivalent model degenerates to a spring with no 
dashpot. 


4. If the coefficient go vanishes, one Voigt element of the 
equivalent model degenerates to a dashpot with no spring. 

5. If go and g, both vanish, the corresponding Voigt 
model is doubly degenerate. 


The second step of the problem is the deter- 
mination of numerical values for the parameters 
of the model. This step is discussed in detail in 
another paper.” The standard, or non-degenerate, 
case will be briefly outlined here. 

It is clear that in order for the operator equa- 
tion and the Voigt model to be mathematically 
equivalent, they must give the same prediction 
for the strain, e(¢), when any given stress sequence 
s(t) is imposed. It is sufficiently general to equate 
the two responses to a stress sequence s(t) =¢?-". 
The general solution of the equation Ps=2Q« 
is the sum of the general solution of the asso- 
ciated homogeneous equation s(¢#)=0 and the 
particular polynomial solution of the complete 
equation. In the same way, the response of the 
model to the stress s(t) is the sum of the general 
response to a zero stress plus a particular poly- 
nomial response to the stress ¢”—'. In order for 
the response of the model to be identical with 
that predicted by the operator equation, the 
parameters of the model must fit the following 
rules. First, the m retardation times of the 
Voigt elements are the negative reciprocals of 
the m roots of the characteristic equation of the 
operator Q; that is r;= —1/x;, where x; are the 
roots of the equation g,x"+---+q.=0. Thus, 
the m retardation times of the model are deter- 
mined by the general solution of the homogeneous 
differential equation. 

In order to complete the specification of the 
model the particular polynomial solution must 
now be used. The particular polynomial solution 
of the equation P(t"—') = 2Q¢ will be of the form: 


2¢(t) =Apt+ayt+aot?+ eee +a,-1/""'. (28) 


The coefficients ao, a1, - - -@,-1 can be determined 
by ordinary methods. The particular polynomial 


2 T. Alfrey, Quar. App. Math. 3 (July, 1945). 
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solution derived from the model will be of the 


form: 
n—1)! > ' - 
i, (Ociiy tii 


(Be? Per Te. (29 
— =) +E Te. (29) 








2e(t)=(m—1)! E— 


It follows that the ” compliances J;, Jo, ---Jn 
must be the roots of the following set of linear 
algebraic equations: 











Qn-1> p Ji, 
1 Ji 
dn2=> i 
—1 x; 
1 > J; (30) 
6,—-3= ; 2. 
(n—1)(n—2) ’ x;? 
1 > Ji 
(n— TY es xn 


The Voigt model is completely specified when the 
compliance |and retardation time of each Voigt 
element are determined. 

The degenerate cases where g, or go vanish can 
also be simply handled." 


H. Operator Equation to Maxwell Model 


A Maxwell model with no degenerate elements 
corresponds to an operator equation in which 
both g, and go vanish. The m relaxation times of 
the m Maxwell elements are equal to the negative 
reciprocals of the m roots of the characteristic 
equation of the operator P. The specification of 
the model can then be completed by the solution 
of a set of m linear equations in m unknowns. 


“Ge = . Nis 
m 











1 => n 
m(m— 1) iret” 
, (31) 
m Ni 
~“aiten = x2) 
1 i 
—ao= 2D vo 
m x;™ 1 
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I. Operator Equation to Impedance Function 


In electrical network theory it is shown that 
the complex impedance, Z(w), is equal to the 
ratio of the Fourier transform of the potential to 
the transform of the current. Since the current 
corresponds to the time derivative of the strain 
it is necessary to note the equivalent expression 
for the Fourier transform of €; that is 


f eritat= (ia) [eerie (32) 


Z(w) can now be calculated from the operators 
P and Q as follows: 


Ps =2Q<«, (33) 


§o+p$i1+ pos = 2[qse+qré+qiét+qoe ]. 


The Fourier transforms of each term are found 
by integrating over ¢ after multiplying by e**. 
Using (32), this procedure yields after collecting 
terms 


[ (iua)?-+ ps (iw) + po] f sets 


= 2 ga(ico)® +g (ies)?-+ qu (ies) +901] f ee-iwidt, (34) 


Thus the form of the impedance function is found 
to be: 


s(t)e-*“'dt 
Z(w) = 





(tw)2 f e(t)e—**'dt 


_ qa(tw)* +g2(tw)? +i (tw) +90 

(Ga) [ Ge)? +- pili) +0] 
_ Qliw) 
~ (tw) P (iw) 





(35) 


J. Impedance Function to Dynamic Modulus 


The relation between the impedance function 
and the response to sinusoidal stresses is, of 
course, particularly simple. At a given frequency, 
w/2mx, where the mechanical admittance, 1/Z(w), 
has the value A+B, the steady-state response 
to the stress, s=5 9 cos wt, will be: 


2é(t) =50A cos wt—soB sin wt, (36a) 


s 
2e(t) tl sin at+—B cos wt. (36b) 
w w 
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In this connection, the relationship between 
sinusoidal stress and strain (in the presence of a 
lumped mass) is sometimes written as: 


Mé+Aé+Be=Sy cos wt, (37) 


or as some obvious variant of this equation. It is 
quite true that such an equation will hold for a 
particular frequency. In order for it to hold at 
any frequency, A and B must become functions 
of frequency, rather than numerical constants. 
This has, of course, been recognized by those 
using Eq. (37). It is, however, possible to write a 
differential equation relating stress and strain 
in which all coefficients are true constants. This 
is the operator equation, Ps = 2Qe, which for this 
problem (after the introduction of a lumped 
mass) becomes: 


Mé+2Qe=P(5so cos wt). (38) 
K. Dynamic Modulus to Impedance 
Function 


The reverse transformation, by which the 
impedance function is derived from experimental 
data on the steady-state relation between sinu- 
soidal stress and strain, is also simple. If the 
amplitude and phase of the strain relative to the 
stress, s=sycoswt, is known, the strain can 
always be written as 


2¢(t) = Cso cos wt+Dsy sin wt, (39) 


where C and D are known functions of w. A, the 
real part of the admittance function, is equal to 
Dw; and B, the imaginary part, is given by Cw. 
Thus the impedance function is 


Z(w) =1/(w(D+10)). (40) 


L. Voigt Model to Impedance Function 


It follows immediately from (C) and (K) that 
the mechanical impedance function, Z(w) = A (w) 
+iB(w) for a material of known Voigt speci- 
fication, will be given by the relation: 


1 1 
emi aamebite 
Zw) 9 5 2j—1€;/w 


In the case of a continuous distribution of elastic 
compliance, J(r), this becomes: 


® 
+—. 


iG aie 





1 1 © I(r)dr 
+f ) (42) 


Zw) 9 Jo r—i/o 
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M. Maxwell Model to Impedance Function 


Similarly, from (D) and (K), it follows that 
the impedance function for a material of known 
Maxwell specification will be given by: 


1 
Z(a)==——, (43) 


nj 1G; 


or, in the case of a continuous distribution G(r), 


om G(r)d 
Z(w) = I _. (44) 


1 





N, O. Voigt or Maxwell Model to 
Operator Equation 


To go from a Voigt model or a Maxwell model 
to the corresponding operator equation is a 
little more involved than the reverse trans- 
formations (G) and (H). A general method for 
these conversions has been developed, however, 
and is described elsewhere.'* There exists in 
addition an obvious method which can easily 
be carried out for models which contain only a 
small number of elements. To illustrate this 
simple method, let us consider a Voigt model 
consisting of two Voigt elements—Gi, 71, Ge, ne. 
The material under consideration is governed by 
the two differential equations: 


$=2Gie1.+2mé1, (45a) 
S= 2G2€2+2n2€2. (45b) 


The second of these can be rewritten in terms of 
the total strain, e=€,+ 2 as follows: 


$=2Gi61. +261, (46a) 
s= 2Gre— 2G2e,+ 2n2éE— 2n2€1. (46b) 


Let us now differentiate Eqs. (46a) and (46b), 
obtaining : 


§=2G16:+2m 61, (47a) 
§ = 2G2é— 2G2€, +2 2€ — 2261. (47b) 


Equations (45)-(47) can now be combined in 
such a way as to remove the strain contribution 
€,, and its derivatives, by multiplying each 
equation by the proper factor and adding. The 
resulting differential equation of the third order 
involves only the total strain and stress and is the 
desired equivalent operator equation. 
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P. Creep Function to Voigt Model 


Since the creep curve of a viscoelastic material 
is a known function of the parameters of the 
specifying Voigt model, it is obviously possible 
to derive the values of these Voigt parameters 
from the shape of the creep curve in those cases 
where the model consists of a discrete set of Voigt 
elements with markedly different +r values. In 
general, however, this will not be the case. It is 
still possible to assume a topological form for the 
model, and then compute a reasonable set of 
parameters for the model. The result is an 
approximation which may be quite good if 
enough adjustable parameters are used, but the 
choice is rather arbitrary. Two different inves- 
tigators may well obtain entirely different 
numerical values for the individual elements of 
their approximating models—even though the 
over-all behavior of both models is quite similar. 
This is an obvious disadvantage. A specification, 
to be most useful, should be unique as well as 
accurate. It is for this reason that method (U) 
will be advanced. 


Q. Relaxation Function to Maxwell Model 


The construction of a Maxwell model from a 
known relaxation curve is precisely equivalent 
to the construction of a Voigt model from a 
known creep curve. Method (V) will be advanced 
to handle this problem in an accurate, unique 
fashion. 


R. Impedance Function to Operator Equation 


Since the impedance function is related to the 
parameters of the operator equation by Eq. (35), 
it is possible in principle to derive the operator 
equation by factoring the known impedance 
function into a ratio of two polynomials in (iw). 
In practice, using an empirical Z(w), this will 
usually be impracticable. The authors see no 
simple means for this conversion. 


S. Creep Function to Maxwell Model 


Simha™ has shown that the distribution 
function G(r), which defines a_ continuous 
Maxwell model, can be derived from the creep 
curve of the material. The method, which is 
based upon the use of Fourier transforms, is 
rather involved and will not be described here. 
Instead, a much simpler, although approximate, 
method will be developed in the next section. 


13 R. Simha, J. App. Phys. 13, 201 (1942). 


709 













































































r ! 
100 F- | —3 
a = 
% 60 f 
i / 
40 
/| 
20 + _ 
: ea } | 
seis | | 
-§ 4 -3 -2 -| 1@] i 2 
log ¢ 
ae ee 
: | 
: | /\ 
oS ee Eee 
be] 
g bs] J 
al 0 
-5 -4 -3 -2 -1 0 { 2 


log T 


Fic. 6. A creep curve and its distributions of compliance. 


T. Relaxation Function to Voigt Model 


Because of the symmetry of the relationships 
among the methods of representation, the 
problem of deriving an equivalent Voigt model 
from a relaxation function is mathematically 
very similar to that of deriving a Maxwell 
model from a creep function. The method of 
Simha, with slight alterations, could be used for 
an exact transformation ; an approximate method 
will be developed in the next section. 


Ill. SOME NEW APPROXIMATION METHODS 


The conversions described in the previous 
section were exact. This advantage, however, is 
offset by the fact that they involve complicated 
computations for materials whose viscoelastic 
models contain a large number of elements. For 
this reason a set of approximation methods will 
be developed in this section. 


U. Creep Curve to Voigt Model 


A creep curve shows how strain varies with 
time at constant stress.’ The strain at any time 
will be made up of two parts—elastic deformation 
and flow. If plotted as 2e versus logit, the 
(retarded) elastic part of the strain will yield a 
sigmoidal curve. The total strain can easily be 


“ H. Leaderman, Ind. Eng. Chem. 35, 374 (1943). 


710 


resolved into the two parts, if the creep data 
cover a wide range of time scale. The viscosity 7 
can be immediately calculated from the flow 
component. The computation of the viscoelastic 
spectrum J(r) from the sigmoidal elastic creep 
curve, however, is not so simple. The following 
method gives an approximation to J(r). We will 
actually derive an approximation to the function 
Lilog +), where L(log r)d log r= J(r)dr. 

Plot the elastic part of the creep curve as 
(2e/s) versus (log t). Determine the slope of this 
curve and plot against (log 7). This curve which 
we may call L’(log r) is a fair approximation to 
the actual distribution function L(log 7). L’ (log 7) 
fails to reflect sharp discontinuities in L (log r), 
but reproduces the general contour. If L(log 7) is 
a smooth distribution, extending over a wide range 
of time scale, then L’(log 7) is a very good ap- 
proximation. If L(log 7) contains discontinuities, 
then L’(log 7) is not so good an approximation. 

Let us investigate the closeness of the ap- 
proximation in the most unfavorable case—that 
of an extremely sharp distribution function. 
Consider a material whose Voigt model consists 
of an instantaneous compliance of 10-* cm?/dyne 
and a retarded elastic compliance of 10-7 
cm?/dyne whose retardation time is 1 second. 
(The viscosity of the retarding dashpot is thus 
10’ poise.) 

The upper part of Fig. 6 shows the exact 
creep curve for such a material. The lower part 
of Fig. 6 shows the derived distribution of com- 
pliance L’(log 7). This can be compared with the 
true distribution L(log r), which is the discon- 
tinuous function. The area under the hump in 
L' (log r) is exactly equal to the correct value of 
Jo, and the maximum of L’(log r) comes at 1 
second. Furthermore, the correct value for the 
instantaneous compliance is obtained. The ap- 
proximation, however, blurs the sharpness of the 
distribution. Actually many viscoelastic mate- 
rials possess distributions which spread over as 
many as twelve decades of log ¢. 


V. Relaxation Curve to Maxwell Model 


In the same way, the distribution of Max- 
wellian relaxation times can be approximated 
from_the relaxation curve. Plot the relaxation 
curve as s/2e vs. log ¢t. Determine the slope at 
each point and plot the negative of the slope 
versus log r. The resulting function, which may 
be called K’ (log 7), is a good approximation to the 
true distribution function K(log r)—where 


K (log r)d log r=G(r)dr. 
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W, X. Creep Function to and from Relaxation 
Function 


The creep function, ¢(¢), and the relaxation 
function, Y(t), are approximately the reciprocals 
of one another at any time. (This is not an exact 
relationship, in spite of the superficial appearance 
of these two functions.) The approximate reciprocal 
relation leads us to a very simple graphical 
method for going from one function to the other: 

Plot the creep function as log (@) versus log t. 
The mirror image of this curve can be obtained 
by reflecting through the horizontal zero-line, 
thus replacing every ordinate by its negative. 
The reflected curve will be an approximation to 
the relaxation function plotted as log (W) versus 
log t. This method obviously works equally well 
in both directions. 

Let us now investigate the closeness of this 
approximation for a very concentrated distribu- 
tion of compliance. We already have the creep 
curve for a Voigt model with J;=10-°, J.=10-’, 
n2= 107. We can also compute, exactly, the re- 
laxation function for this model. This is shown 
in the upper part of Fig. 7 (solid curve). Using 
the known creep curve and the approximation 
method here advanced, we obtain the approxi- 
mate relaxation function shown in the upper 
part of Fig. 7 (dotted line). 

It may also be of interest to see how closely 
this method allows an estimation of the distri- 
bution of Maxwellian relaxation times, K(log 7), 
from a creep curve. Using the same model as 
before, we can obtain three different distribution 
functions: 


1. The exact form of the Maxwell model. This is in- 
dicated by the infinitely sharp distribution curve (1) in 
Fig. 7 (lower part). 

2. K’ (log r) can be approximated from the exact relax- 
ation function by means of (V). 

3. The approximate relaxation function can be approxi- 
mated from the creep curve by (W), and then (V) can be 
applied as a second approximate step. The result is shown 
in Fig. 7 (lower part). 


IV. CRITICAL COMPARISON OF THE SEVEN 
METHODS OF REPRESENTATION 


Having discussed, in previous sections, seven 
methods of representing the mechanical proper- 
ties of viscoelastic materials, and the important 
mathematical relationships among these meth- 
ods, we now proceed to a critical evaluation of 
the methods. 

First of all—as indicated in the discussion— 
Class I methods would seem to be preferable to 
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Fic. 7. A relaxation curve and its distribution of 
Maxwellian relaxation times. 


the direct experimental functions of Class II 
from the standpoint of describing the properties 
of a material. This is because Class | methods 
are in the nature of differential equations which 
hold for any and all stress and strain sequences, 
whereas Class II methods are special solutions 
applying to only one type of experiment. A 
Maxwellian material, for example, is described 
in methods (A)-(C) by the differential equation : 


é= (1/G)8+(1/n)s. (48) 


In the special case of a constant strain, e, sud- 
denly imposed at the time ¢=0, the stress follows 
the specific equation 


s=Gee~"!", (49) 


The difference between Eq. (48) and the relaxa- 
tion function, y(t), is that the former holds for 
any sequence of stresses and strains, while the 
latter holds only for the one special experiment. 
For this reason (even though Eq. (49), together 
with the superposition principle, does implicitly 
define the properties of the material) we will 
limit the rest of this section to a relative evalu- 
ation of Class I methods as specifications of 
properties, followed by a relative evaluation of 
Class II methods as experimental tools. 
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A sharp distinction can now be made between 
methods (A) and (B) on the one hand, and (C) 
and (D) on the other. From the standpoint of a 
phenomenological study of the mechanics of 
continua, methods (C) and (D) are superior to 
(A) and (B). This is because in any phenomeno- 
logical study of mechanical behavior the only 
observable quantities are the components of the 
(total) stress and strain tensors; the ‘“‘strain 
contributions” of (A) and the “‘stress contribu- 
tions’’ of (B) are, strictly speaking, unob- 
servable quantities. This very fact, however, 
makes (A) and (B) more suitable than (C) and 
(D) for the formulation of molecular theories of 
deformation.'>—'§ In a molecular study, each con- 
tribution to the strain may often be identified 
with some specific molecular process; and hence 
the strain contributions as well as the total strain 
can be said to possess a physical significance. Like- 
wise, some authors!** have identified the various 
stress contributions of a Maxwell model with in- 
dividual ‘‘molecular mechanisms of supporting 
stress.” 

Our conclusion is that methods (A) and (B) 
are more suitable than (C) and (D) for a study 
of the molecular mechanisms of viscoelastic 


deformation, but that within the discipline of ° 


the mechanics of continua (C) and (D) are to be 
preferred. 

From the standpoint of the mechanics of con- 
tinua, there is a practical as well as an absolute 
criterion which must be considered. This is the 
ease with which a method can be used to predict 
behavior in special cases; i.e., the ease with 
which the fundamental differential equations can 
be solved subject to definite constraints. We may 
compare Class I methods on the basis of this 
practical criterion, as follows: 

Method (A), the Voigt specification, is ex- 
tremely well adapted to the solution of all 
problems in which the stress is specified as a 
function of time, and the resulting strain is 
desired. The reverse type of problem, in which 
a known strain is enforced upon the material and 
the stress is to be calculated, is very difficult to 
handle by method (A). 


% T. Alfrey, J. Chem. Phys. 12, 374 (1944). 
16 A. V. Tobolsky and R. D. Andrews, J. Chem. Phys. 
~ 13, 1 (1945). 

7 R. F. Tuckett, Trans. Faraday Soc. 38, 310 (1942). 

18 R. F. Tuckett, Trans. Faraday Soc. 40, 448 (1944). 

1” R. Simha, J. Phys. Chem. 47, 348 (1943). 

20 R. Simha, Ann. N. Y. Acad. Sci. 44, 297 (1943). 

* A. V. Tobolsky, R. E. Powell, and H. Eyring, The 
Chemistry of Large Molecules (Interscience Publishers, Inc., 
New York, 1943), 
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The Maxwell specification is the exact inverse. 
Problems in which the strain is specified can be 
solved very easily ; those in which s(t) is specified 
are very difficult. 

The operator equation (C) is intermediate in 
character. Problems involving a known s(t) or e(t) 
can both be solved readily, although perhaps not 
quite so simply as with the more applicable of 
methods (A) and (B). 

If the impedance function is known, problems 
of both types can be solved by the use of Fourier 
or Laplace transforms. 

We have observed that both the Voigt speci- 
fication (A) and the Maxwell specification (B) 
are amenable to molecular interpretation. The 
important question now arises: Can either of 
these two methods be said to be superior to the 
other from the standpoint of a molecular theory? 
This is a highly debatable question the answer 
to which at the present time would necessarily 
be based very largely upon the personal preference 
of the investigator. The authors’ personal pref- 
erence is for the Voigt specification. From this 
point of view, the molecular interpretation of 
viscoelastic behavior may be stated roughly as 
follows: In the presence of a stress a number of 
different molecular processes take place.” Each 
molecular process represents some distinct mode 
of molecular rearrangement. Each process has 
its own characteristic time dependence and 
thus results in a strain contribution which is a 
characteristic function of the time. The total 
strain in the sample is the sum of the contri- 
butions from all the different molecular mecha- 
nisms of response to stress. Using this formulation 
of the problem, a specific molecular theory of the 
viscoelastic response of one type of material— 
the amorphous linear high polymer—has been 
developed.'5 

In terms of the Maxwell specification of 
properties, the molecular description’ '® of 
viscoelastic behavior can be stated in the fol- 
lowing manner: The stress on a given cross section 
of the material is supported by interatomic or 
intermolecular forces exerted upon the atoms 
and molecules lying in the cross section. These 
forces may be given the generic name of ‘“‘bonds” 
connecting the material on the two sides of the 
cross section. Each of these bonds will have its 
own stiffness and its own characteristic rate of 
relaxation. A distribution of elastic relaxation 

2S. Bresler and J. Frenkel, Acta Physicochimica, 11, 
485 (1939). 


23 A.V. Tobolsky, I. B. Prettyman, and J. H. Dillon, 
J. App. Phys. 15, 380 (1944). 
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times is thus correlated on the molecular scale 
with a distribution in the relaxation times of the 
interparticle bonds. 

The experimental evaluation of the viscoelastic 
properties of a material can be conveniently 
divided into two parts—the relatively fast re- 
sponse (say faster than one second) and the rela- 
tively slow response (say slower than one sec- 
ond). In principle, either method (E), (F), or (G) 
can be used to evaluate both aspects of the visco- 
elastic behavior. In practice, this is not true. 

For the rapid part, determination of the 
steady-state response to sinusoidal stresses (as a 
function of frequency) is the most convenient 
method. Creep and relaxation studies of the 
response which occurs in a time scale of 0.01 or 
0.001 second are rather impracticable. 

On the other hand, creep and relaxation tests 


are much better suited to the evaluation of the 
slow response than are vibration tests. The 
dividing line of 1 second, above, is purely arbi- 
trary, and there is no reason why the range of 
applicability of the ‘‘slow’’ and “‘fast’’ methods 
cannot be made to overlap. 

The most promising experimental approach 
would seem to be the following : Map out the fast 
response of the material by means of sinusoidal 
stresses. Then map out the slow response by a 
creep curve or a relaxation curve. Finally, syn- 
thesize the two sets of data and express in terms 
of one of the fundamental methods of repre- 
sentation—e.g. a Voigt model. The resulting 
specification of the viscoelastic properties will 
cover the entire range of time scale—from the 
fastest frequency of the sinusoidal tests to the 
total duration of the creep curve. 





A Double-Yoke Plate Permeameter 


WILLIAM F. Brown, JR., AND SIDNEY M. RUBENS 
Naval Ordnance Laboratory, Washington, D. C. 


(Received June 25, 1945) 


A double-yoke permeameter, with which both normal and ideal permeabilities of moderately 
large steel plates can be measured, is described. The technique for making these measurements 


is also given. The permeameter is suitable for measuring plates several feet long, u 
wide, and from 3” to 14” thick; two approximately identical test plates are requir 


T is sometimes necessary to measure the mag- 
netic properties of iron or steel plates }’’ to 
13’ thick, several feet long, and up to 1’ wide. 
Although ring specimens or specimens suitable 
for use in standard permeameters may be 
machined from such plates, the procedure is 
laborious, destroys the plate, and in the case of 
soft materials affects the magnetic properties to 
an unknown degree. 

The double-yoke permeameter, illustrated in 
Fig. 1, was designed to make such measurements 
directly on plates. This device is a magnetic yoke 
in which a pair of test plates, both of the same 
kind of steel, form part of the magnetic circuit. 
The two test plates, which may be as long as 
desired and up to one foot wide, are mounted in 
the yoke so that over a considerable part of their 
length the plates are magnetized uniformly and 
in the same direction. This makes it possible to 
measure both B and H directly with suitable 
test coils. When the average permeability of 
large specimens is measured in this way, three 
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to a foot 


disadvantages usually present in measurements 
on ring specimens of the usual size (outside 
diameter <3’). are. obviated: inhomogeneity 
over short distances, anisotropy, and uncertainty 
of state of mechanical strain. 

The following terminology is convenient in 
describing the different states of magnetization 
of a specimen and the methods of producing 
them. To erase previous history, a specimen is 
subjected to a series of applied fields of alternat- 
ing sign and gradually diminishing magnitude; 
during this procedure, an additional constant 
field may be superposed. The state of mag- 
netization produced by this procedure is called 
the “‘ideal’’ state; and we shall call the procedure 
for producing it “‘idealization,’’ the field used to 
erase previous history the “‘idealizing” field, and 
the superposed constant field the “parameter” 
field. The idealizing field must begin with a high 
enough value to produce a result independent of 
the previous magnetic state, and it must decrease 
by small enough steps so that if they were to be 
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Fic. 1. Construction of 
the double-yoke plate per- 
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made still smaller, no change in the final mag- 
netization value would result. Idealization in 
zero-parameter field produees zero magnetization 
and is called “demagnetization.” ‘“‘ Normal’’ mag- 
netization is produced by first demagnetizing the 
specimen and then applying a field of fixed 
magnitude alternately in opposite directions 
until a cyclic state is attained. 

The double-yoke plate permeameter was 
designed primarily to measure normal perme- 
ability at low fields, but it may also be used to 
‘measure ideal permeability. 


DESIGN AND OPERATION 


The construction of the permeameter is shown 
‘in Fig. 1. The yoke consists of two U-shaped 
pieces of mild steel 3 feet long which have been 
sawed from a single plate 4 inch thick. The ends 
of these yoke pieces are milled smooth on each 
face so that they make good contact with the 
two test plates. The middle section A of each 
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show ‘“‘H”’ test coil. 





of the U-shaped pieces is wound uniformly with 
two electrically separate windings, which provide 
the primary flux in the yoke. Each winding con- 
sists of one layer of number 18 enameled wire. 
One of these windings may be used to provide a 
steady field; the other furnishes a field which is 
varied as may be required for demagnetizing, 
idealizing, or measuring. The yoke pieces are 
clamped together with the two test plates. The 
test plates are 3 feet long, are from 4 to 12 inches 
wide, and may be of any desired thickness. 
Several turns of fine wire are wound around 
each test plate. These serve as ‘B’ test coils. H is 
measured with a long flat rectangular coil wound 
on a Bakelite form. For plates 10 inches wide or 
more, an ‘H’ coil 12 inches long, 6 inches wide, 
and 3 inches thick is used to measure fields up 
to 1 oersted. Its effective NA (turns times cross- 
sectional area) is 4.52104 turns cm?. This coil 
lies in the space between the two test plates with 
its longitudinal axis parallel to the long dimen- 
sion of the plates. For narrower plates, an ‘H’ 
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coil 2 inches wide having about the same effective 
NA has been used. 

The yoke is assembled as shown in Fig. 1. The 
primary windings of the two sides of the yoke 
are connected so that the flux from the two sides 
passes through the test plates in the same direc- 
tion. The small windings on the short legs of the 
U-pieces each consist of 100 turns of number 18 
enameled wire wound on insulating forms. These 
can be moved along the short legs if necessary 
to compensate for flux leakage and improve the 
uniformity of the flux distribution over the 
central 12 inches of the test plate. The relative 
flux distribution obtained in a medium steel test 
plate ~ inch thick and 12 inches wide is shown 
in Fig. 2. This flux distribution was measured 
with a ‘B’ test coil whose distance from the 
center of the plate could be varied. 

In the measurements made on test plates of 
structural steels from 4 to 12 inches wide, no 
special precautions were required to obtain a 
sufficiently uniform flux distribution in the test 
plates. The plates were sawed to the desired size 
and filed clean where contact was to be made 
with the yoke pieces. They were clamped to the 
yoke as tightly as possible with brass C-clamps. 
For the 12-inch wide plates, the four small com- 
pensating coils just fitted between the long legs 
of the yoke and the test plate, so that their 
positions could not be adjusted. Nevertheless 
sufficiently uniform flux distributions in test 
plates of various thicknesses and permeabilities 
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Fic. 2. Test of uniformity of magnetization of test plates: 


flux distribution in }’’ medium steel test plate. 
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TABLE I. Uniformity of flux distribution in test plates. 








Distance 
from center 

of plate Medium steel Medium steel High tensile steel 

(inches) 4” K11.5” X36” a” x4” X36” a” X12” K36” 


6 0.5 —3. 
4 A 
2 
0 
—2 
—4 
—6 
Smallest 
detectable 
fractional 
difference 1/600 1/200 1/700 


% deviation from flux at center of the test plate 
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were obtained by ‘simply connecting the com- 
pensating coils in series with the main windings. 
Sample data on the flux distributions in 3 dif- 
ferent test plates are given in Table I. The 
smallest detectable difference in fluxmeter read- 
ings, expressed as a fraction of the reading, is 
given at the bottom of the table. If this uni- 
formity had not been so easily obtained with the 
wide plates, it would have been necessary to 
adjust the flux distribution by adjusting the 
currents in the compensating windings. 

In order to measure the flux density distribu- 
tion across a test plate at the center, small holes 
were drilled through a medium steel plate of 
dimensions }’”’ X 12” X36” so that three small ‘B’ 
coils could be wound through them. These coils 
were coplanar and were each 2 inches wide. One 
was at the center of the plate, and the centers of 
the other two were +3 and —5 inches from the 
center of the plate. B as measured by each of 


TABLE II. Variation in B across a test plate at the center. 











Position of ‘B’ coil B(gausses) 
At center of plate 141 
3 inches from center 132 
—5 inches from center 138 
Around entire plate 135 








these coils and by a test coil wound around the 
entire plate at the center under a given state of 
normal magnetization are listed in Table II. 
The average of the three small B coil measure- 
ments is 137, which differs from the value ob- 
tained by the coil surrounding the entire plate 
by less than 2 percent. 

The uniformity of H in the region between the 
plates was tested by measuring with the 12” 
long ‘H’ coil at several positions. The coil was 
moved both across the test plates and along 
them. The flux through the ‘H’ coil was constant 
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Fic. 3. Contributions of near poles to longitudinal mag- 
netizing force H at a point P in lower test plate, to longi- 
tudinal magnetizing force H’ at a point P’ above it, and 
to longitudinal magnetizing force H” at a point P’ below 
it. (a) Comparison of H’ and H: lower plate only present. 
(b) Comparison of H’ and H: both plates present. (c) Com- 
parison of H” and H. For identical pole distributions on 
equidistant surfaces, the contributions at P from points 
connected to it by solid lines are equal to the contributions 
at P’ or P” from points connected to it by solid lines; the 
contributions corresponding to dashed lines are unequal. 


to within 4 percent when the center of the coil 
was within 6 inches of the center of the plates 
longitudinally and 3 inches of the center trans- 
versely. The ‘H’ coil was also placed on top of 
the top plate, at the center, and the value of H 
was found to be a percent or two lower than in 
the space between the plates. 

The following procedure is used to measure 
normal permeability. The plates are demag- 
netized, the desired magnetizing field is applied, 
and the plates are brought into a cyclic state. 
Then the magnetizing current is reversed; B and 
H are measured by observing the deflections of 
suitable fluxmeters connected to the ‘B’ and ‘H’ 
coils. Directly B/H=u. 


VALIDITY OF THE H MEASUREMENT 


With any specimen which is at least approxi- 
mately uniformly magnetized over part of its 
length, the measurement of B presents no dif- 
ficulties; the more difficult problem is_ the 
measurement of H. Except in the case of ring 
specimens and ellipsoids, for which H can be 
calculated, it is necessary to rely on the con- 
tinuity of the longitudinal component of H across 
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the surface of the specimen. But the quantity 
actually measured is the average value, H’,y, of 
this component over a finite region adjacent to 
the specimen, and in many cases this average may 
differ by 100 percent or more from the average 
value, Hy, over the cross section of the specimen. 
In the present permeameter, suppose first 
that only the lower test plate is in place. The 
magnetizing force H at a point either inside or 
outside the plate, and near its center, may be 
analyzed on the basis of the sources into three 
parts: (1) that directly due to the windings; 
(2) that due to poles distant, say, more than 5” 
from the point in question; (3) that due to poles 
close to this point. Contributions (1) and (2) are 
approximately the same in the plate and in the 
region above it but within }” of its top surface. 
It is contribution (3) that is likely to have 
serious effects. As is indicated in Fig. 3(a), a 
pole at a point C on the upper surface of the 
plate contributes equally to the longitudinal mag- 
netizing force H’ at a point P’ above the top 
surface of the plate and to the longitudinal mag- 
netizing force H at a point P an equal distance 
below it; but a pole at a point D on the lower 
surface contributes less to H’ than to H. Hence 
if the contribution of local poles is appreciable, 
H’,, may differ significantly from J7,,. (The con- 
tribution of the volume pole density is neglected 
in this discussion, since it is usually small com- 
pared with the contribution of surface poles.) 
Addition of the upper test plate improves this 
situation considerably. An ideally simple case, 
in which the pole distributions on the four 
surfaces are identical and the distance between 
the plates is equal to their thickness, is shown in 
Fig. 3(b). (The former condition will be closely 
approximated, if contribution (3) is small com- 
pared with contributions (1) and (2).) The con- 
tribution at P’ from A, B, and C is equal to the 
contribution at P from B, C, and D; the dif- 
ference between H’ and H is due entirely to 
the inequality of the distances DP’ and AP. 
Since the contributions of such points as D to 
H’ and of such points as A to H are a small part 
of the whole contribution of local poles, which 
in turn is only a fraction of the magnetizing 
force, the error in equating H’,, to Hay should be 
very small under these ideal conditions. In the 
actual permeameter these conditions are not 
fully realized; but the error can be estimated by 
also measuring H”,,, the average magnetizing 
force over an area just below the lower plate or 
above the upper one. Fig. 3(c) shows that the 
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difference between H” and H is due largely to 
the difference between the contributions of A to 
H" and of D to H; these contributions differ by 
an appreciable fraction of the whole contribution 
of local poles. Therefore an agreement of H’ and 
H” to within 2 percent is strong evidence that 
the contribution of local poles is a small fraction 
of the magnetizing force, that the deviations 
from ideal conditions will therefore affect only 
a small part of the magnetizing force, and that 
H’ and H are equal to within considerably less 
than 2 percent. 

These arguments do not apply to double-yoke 
permeameters in which the specimens are bars 
with widths and thicknesses approximately 
equal. Therefore this plate permeameter appears 
to have advantages over such permeameters as 
far as reliability of the H value is concerned, 
especially since the H’’ measurement affords a 
method of estimating this reliability. For best 
results, the flux should be made as nearly uniform 
as possible by adjustment of the auxiliary coil 
positions and currents; the H’ measurement 
should be checked by a measurement of H”; 
and the pair of test plates used for the final 
measurements should be selected ones that give 
B and H” values very nearly equal for the two 
plates. 


COMPARISON WITH RING SPECIMEN DATA 


To check the reliability of the normal-perme- 
ability values obtained with the double-yoke 
permeameter, a ring specimen 4 inches in 
diameter was cut from a steel plate one-half inch 
thick, which had previously been measured in the 
double-yoke. The ring specimen was wound and 
measured. The results of both the ring and the 
double-yoke measurements appear in Fig. 4, in 
which y» is plotted against H. Although the 
points agree to within 1 percent, the reliability 
of the plate permeameter is not necessarily this 
good. The ring measures an average perme- 
ability for all directions of the magnetizing force 
in the plane of the sample, whereas the plate 
permeameter measures the longitudinal per- 
meability of the plates; therefore if both methods 
were perfect, they would give different results 
for anisotropic materials. The difference for this 
type of plate is estimated to be about 6 percent. 
In the results of Fig. 4 this difference is evidently 
compensated either by other differences between 
the quantities being measured, such as dif- 
ferences in the space-average permeabilities due 
to inhomogeneity, or by errors in one or both 
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Fic. 4. Comparison of ring specimen and double yoke 
measurements: variation of rmeability with mag- 
netizing force for a sample of medium steel. Crosses, 
double yoke; circles, ring specimen. 


methods. It is probably safe to conclude that the 
plate permeameter measurements are reliable to 
within 5 percent. 


OTHER APPLICATIONS 


The double-yoke permeameter may also be 
used, with somewhat less precision, to measure 
the ideal permeability of one of a pair of plates. 
A steady parameter field is applied by one of the 
primary yoke windings. The other winding 
applies the idealizing field. When the idealization 
is completed, the ‘H’ coil is rotated about an 
axis normal to the plates. (The ‘H’ coil was 
mounted on a circular cardboard disk which lay 
on the lower test plate, and the disk could be 
readily rotated by hand.) The resulting fluxmeter 
deflection gives H. The clamps are now removed. 
The top plate is gently lifted off the yoke pieces 
and the ‘B’ coil is slipped off the plate. The 
resulting fluxmeter deflection is proportional to 
B. In this way, also, states resulting from special 
procedures, such as incomplete idealization (in 
which the initial field is not high enough to erase 
previous history entirely, or the steps in reduc- 
tion of the field are too large to result in a truly 
ideal state), may be studied. A small ‘H’ coil 
placed on top of the upper plate is convenient 
for approximate measurement of the idealizing 
field. 

B and H values on a hysteresis loop may be 
measured by an obvious modification of the 
technique for obtaining normal magnetization 
data. 

The authors gratefully acknowledge the valu- 
able suggestions given by Mr. J. H. Sweer and 
the assistance rendered by Mr. J. R. Wright in 
constructing and testing this device. 
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An Electron Diffraction Camera for the Study of High Temperature 
Surface Reactions 


EARL A, GULBRANSEN 
Westinghouse Research Laboratories, Pittsburgh, Pennsylvania 


(Received June 28, 1945) 


An electron diffraction camera is described for the study of physical and chemical reactions 


occurring on the surfaces of metals and alloys in situ at elevated tem 


ratures. The specimen 


is mounted in a metal block and shielded from the surrounding cold surfaces by a series of 

metal shields. The metal block contains an internal electrical heater which allows the tem- 

perature of the specimen to be accurately controlled at temperatures of 800°C. Gas purification 

chains are included for the preparation of pure oxygen, hydrogen, etc. Several typical electron 
‘ diffraction reflection studies of the oxidation of iron and stainless steel are included. 


INTRODUCTION 


HE study of the structure of the reaction 

product resulting from chemical and 
physical reactions occurring on surfaces is a 
problem of increasing importance and one in- 
volving many difficulties. Electron diffraction 
methods have been applied recently with con- 
siderable success. 

During the past several years we have under- 
taken a program for the study of the kinetics 
of formation of oxide and other films on metals. 
In attempting to interpret the data thus obtained 
it was found necessary to study the structure of 
films by electron diffraction methods. 

For a complete interpretation it was felt that 
the specimens of metal should be studied in the 
electron diffraction camera itself and at given 
temperature and gas atmospheres. Thus, for 
example, the specimens could be oxidized in air 
at 600°C and the diffraction picture taken after 
evacuating at the same temperature. This 
would assure a minimum structural change and 
would allow a series of consecutive pictures to 
be taken without exposing the sample to the 
outside atmosphere or allowing the specimen to 
cool. The effect of heat and gas treatments as 
well as other processes could be studied in situ. 

An electron diffraction camera for this type 
of study requires (1) a furnace capable of being 
controlled at temperatures up to 800°C, (2) a 
series of attached gas chains for purifying the 
gases of interest, i.e., Oo, air, He, etc., (3) a rapid 
pumping system, a pressure measuring system 
and a valve system for isolating the camera from 
the pumps, and (4) a film or plate chamber in 
which 20 to 30 diffraction pictures could be 
obtained befo_e changing the film. 

Two cameras have been described which 
include specimen holders capable of being heated 
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to high temperatures. The first is the camera 
described by Yearian and Howe! which included 
a specimen holder for heating the specimens to 
600°C and which also allowed the specimens to 
be cooled to liquid air temperatures. The second 
is the camera described by Jackson and Quarrell.? 
A carbon arc was used for heating the end of a 
tube on which the metal specimen was mounted. 
The temperature of the surface was determined 
by a thermocouple. A small quantity of gas 
could be admitted and the reaction product 
studied in situ at the temperature desired. 


APPARATUS 


The essential features of our camera will be 
here briefly discussed. A schematic diagram of 









CATHODE 

—_ 

DISCHARGE 
TUBE 

ANODE 


THERMOCOUPLE GAGE 1 





LEAK VALVE 





MANIFOLD 
AND 








VALVE 











ee 
FOCUSSING 
COIL 
VALVE L DEFLECTING 
Tv 
DIFFUSION Y eins 
PUMP - = t- 
FURNACE 
SECTION 
“SAMPLE 
MANIPULATOR —}- 
| 


AND FURNACE 


FLUORESCENT SCREEN 
AND FILM HOLOER 








— 1 





GAS CHAINS 


Fic. 1. Schematic diagram, electron diffraction camera. 


1H. J. Yearian and J. D. Howe, Rev. Sci. Inst. 7, 26 (1936). 
2 Jackson and Quarrell, Proc. Phys. Soc. London, 51, 
237-43 (1939). 
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Fic. 2. Preparation chains, 


the camera is shown in Fig. 1. A cold cathode 
source is used to obtain the electrons. Hydrogen 
gas is used in the discharge tube, the pressure 
being regulated by varying the pressure of gas 
behind a porous porcelain leak. For a definite 
voltage, the current in the discharge tube is 
dependent only on the gas pressure. Using a 
voltage of 47.5 kv and a current of 50 micro- 
amperes in the discharge tube, exposures of 5 
to 10 seconds are necessary with a 0.004-inch 
anode aperture. The beam current in the camera 
section is of the order of 10-* ampere. An axial 
symmetrical magnetic field is used for focusing 
the electrons on the fluorescent screen and film. 

_Two sets of deflecting plates are used for 
moving the beam and, in combination witha 
movement of the specimen holder, the angle of 
incidence of the electrons relative to the specimen 
can readily be changed. 

The film chamber is designed for use with roll 
film. Twenty to thirty exposures can be taken 
without reloading the chamber. Above the film 
a fluorescent screen is placed which acts as a 
shutter to control the exposure time. This is 
manipulated by a lever and a sylphon bellows. 

A metal oil diffusion pump having a speed of 
250 1/sec. is used in combination with a Megavac 
fore pump to obtain the high vacuums. (See 
Fig. 2.) A three-inch valve is placed in the 
manifold between the camera and the diffusion 
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pumping, and leak system. 


pump. The system can also be pumped through 
an auxillary line direct to the fore pump. A 
vacuum of 10~ to 10-5 mm of Hg is maintained 
in the main section of the camera while a 
pressure of about 0.01 mm Hg is maintained in 
the gas discharge tube. 

The pressure in the system can be measured 
either by a manometer, by a thermocouple gage, 
or by a McLeod gage. The arrangement of the 
pressure measuring -system, gas preparation 
chain, and the vacuum system is shown in Fig. 3. 
Hydrogen gas is purified by (1) drying over 
CaCle, (2) removal of Oz by reaction with H, 
over platinized asbestos, and (3) removal of HzO 
and acid constituents by means of CaCl, and 
ascarite. Oxygen is purified by (1) drying over 
Ba(ClO,)e, (2) conversion of He, CO, and other 
combustibles to HO, COs, etc., by the use of 
platinized asbestos, and (3) removal of H,O, 
COs, etc., over Ba(ClO,4)2 and ascarite. 

To facilitate the removal of water and other 
volatile materials from the film a small tray of 
P.O; is placed below the film in the camera. The 
film is partially dried before using by pumping 
in a separate vessel. A special process 6-122 roll 
film with non-curling properties is used in the 
camera although Verichrome film has been used 
during the emergency. The high voltage is 
supplied by a conventional type of rectifier, the 
input of the high voltage transformer being con- 
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Fic. 3. Electron diffraction camera furnace and manipulator. 


trolled by a voltage regulator to +4 volt. The 
ripple at 50-microampere discharge tube current 
is less than 0.01 percent. The voltage is calibrated 
over the range by measuring the voltage drop 
across a 10,000-ohm resistor in series with a 1000- 
megohm resistor. The latter is sealed in a glass 
tube and kept in a nitrogen atmosphere. This 
external calibration is compared to a calibration 
based on the diffraction lines of evaporated 
aluminum or evaporated gold. Over a period of 
three years the high voltage has been found to 
be constant to better than 3 or 1 percent. 

The specimen holder, furnace, and manipulator 
are shown in Fig. 3. Two motions are possible. 
The first motion is a rotation so that several 
specimens can be examined successively. The 
second motion is an in and out motion perpen- 
dicular to the beam. This allows a’change of the 
incident angle for reflection studies. The rotation 
is accomplished by a‘rotating grease vacuum 
seal and can be set to an accuracy of 1/10 of a 
degree. The in and out motion operates with a 
micrometer screw thread on a sylphon bellows 
and can be set to one ten-thousandths of an inch. 

Details of the furnaces are shown in Fig. 4. 
Two furnaces will be considered. The first is 
shown in Fig. 4a. It consists of a pure silver 
block attached to the manipulator head by 
means of a thin-walled stainless steel section. 
The silver is furnace brazed to the stainléss steel. 
A silver cap with the-appropriate holes for the 
beam is fitted over the front face and part of the 
side. Thus the specimens are nearly surrounded 
by a uniform temperature and the surface tem- 
perature is close to the temperature of the 
furnace. The heater wire is wound on a Lavite 
core which is fitted into a recess in the silver 
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block. A chromel-alumel thermocouple is placed 
in a hole drilled centrally in the silver block to 
a point } inch from the front surface. The leads 
from the thermocouple and heater are taken out 
through stainless steel tubes. The furnace proper 
is attached to the head of the manipulator by a 
vacuum tight fuse wire gasket. The advantage 
of this arrangement is that the degassing of the 
wire and Lavite core does not effect the vacuum 
of the camera. The head of the manipulator is 
cooled by a stream of water. If it is necessary to 
operate the heater at extreme temperatures, the 
inside of the silver block can be evacuated. To 
reduce heat leakage, radiation shields are placed 
around the furnace and at other suitable points. 
The temperature of the surface of a steel sample 
has been checked by another thermocouple 
against that given by the block. The temperature 
of the surface has been found to be within 2 to 
14°C of the block temperature at a temperature 
of 500°C depending on the presence or absence 
of a gas atmosphere. 

A second furnace is shown in Fig. 4b. This is 
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Fic. 4. Furnaces for electron diffraction camera. 
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designed for one specimen and for use at tem- 
peratures up to 800°C. It is constructed of one 
piece of ‘Inconel’’ bar and is mounted on the 
manipulator head in the same manner as the 
silver furnace. Because of the poorer heat con- 
ductivity of the ‘‘Inconel”’ an additional thermo- 
couple can be mounted on the inside surface of 
the cap to check the cap surface temperature. A 
disappearing hot wire type of optical pyrometer 
is used for measuring the temperature of the 
furnace cap or specimen. The temperature of 
the surface is calibrated against the temperature 
of the block over the temperature range under 
which the furnace is used. 

For transmission type of specimens a special 
head can be fitted to the manipulator either for 
screen or for slit mounting of the specimens. In 
case heat treatment of the transmission speci- 
mens is desired, a steel block slightly larger than 
the usual 3-inch diameter and 3-inch long speci- 
men is fitted into the specimen chamber of the 
furnace. A hole has been drilled through the 
block parallel and near to the front face to 
accommodate the screen and to permit the 
passage of the incoming and diffracted electrons. 
This is useful for the heat treatment of evapo- 
rated metal films on silica or stripped oxide films. 


CALIBRATION 


Figure 5 shows schematically the variables 
involved in the analyses of reflection patterns. 
The poly-crystalline specimen at a distance L 
from the film diffracts electrons giving rings of 
radii Ray. The undiffracted beam is shown 
striking the film behind the shadow edge. 
Bragg’s law relating the diffraction angle to 
wave-length and lattice spacing can be applied 
to electron diffraction as well as to x-ray dif- 
fraction. For small angles Bragg’s law simplifies 
to Rari=AL/dnx. Here d is the wave-length and 
dnxi is the lattice spacing effective in producing 
the diffraction 4,2. 

The instrument can be calibrated by directly 
determining the voltage and thus \ or by measur- 
ing the diffraction pattern of known materials 
and calculating the values of \. Table I shows 
a comparison of the d;x: of evaporated aluminum, 
(taken by reflection) calculated from the voltage 
calibration and those listed in Hanawalt, Rinn, 
and Frevel’s* tables. The agreement is good to 
about 1 percent. This is what is expected because 


* Hanawalt, Rinn and Frevel, Ind. Eng. Chem. Anal. 
Ed. 10, 457-512 (1938). 
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Fic. 5. The reflection method. 


of the uncertainty of the plate to specimen 
distance due to the width of the specimen. The 
width of the specimen is 1 cm and the specimen 
to plate distance is 48 cm. 


PREPARATION OF SPECIMEN FOR 
REFLECTION STUDIES 


The technique of preparing the surface for 
study is extremely important since the diffraction 
of the electrons impinging on to the surface is 
markedly affected by (1) the flatness of the 
surface, (2) the area exposed to the beam, (3) the 
nature and size of the asperities present, and (4) 
the presence of impurities. According to Thom- 
son and Cochrane‘ the aim is to obtain a clean, 
flat, and slightly matted surface of the material. 
For many applications the specimen is ground 
on various grades of polishing paper using either 


TaBLE I. Diffraction pattern of evaporated aluminum. 











Electron diffraction X-ray tables 

Radius 

in cm dazi in A Intensity dari in A Intensity 
1.09 2.31 Strong 2.33 Strong 
1.25 2.015 Strong 2.02 Strong 
1.78 . 1.415 Strong 1.43 Strong 
2.07 1.215 Strong 1.219 Strong 
2.17 1.16 Weak 1.168 Weak 
2.73 0.922 Weak 0.928 Weak 
2.81 0.895 Weak 0.905 Weak 
3.08 0.820 Weak 0.826 Weak 
3.29 0.765 Weak 0.778 Weak 








4 Thomson and Cochrane, Theory and Practice of Electron 
Diffraction (The Macmillan Company, New York, 1939). 
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Fe FesxOu 





Fic. 6. Oxidation pure iron 300°C at 1.5 cm air. 
a. Vac. at 25°C. b. Oxidized 47 min. 


benzene, kerosene, or a solution of paraffin in 
kerosene. These solvents serve as a lubricant and 
aid in the prevention of oxidation of the metal. 
After removing the lubricant the specimens can 
be washed and dried in absolute alcohol. The 
specimen may be gently etched in order to 
remove surface oxide films and the surface layer 
of the metal. However, extreme care must be 
observed to remove traces of impurities left on 
the surface. 

Recently Vernon and co-workers have studied 
the effect of polishing on the surface film formed 
on chromium nickel 18-8 stainless steel by 
chemical analyses of the oxide film. They find 
that the higher the degree of polishing the thicker 
is the resulting oxide film and the greater the 
change in the chromium to iron ratio in the 
surface oxide. This work indicates that a high 
degree of polishing is to be avoided if the surface 
is to take part in subsequent chemical reaction. 

To prepare flat, clean, and abraded surfaces 
Mr. G. Rose of our laboratories has developed a 
precision abrader. To minimize oxide film forma- 
tion and changes in the composition of the surface 
layer it was felt that the heat generated in the 
abrasion process should be reduced to a minimum 
without the use of lubricants. This frictional heat 
can be partially removed if fresh surface is con- 
tinually exposed to the metal. This avoids metal 
to metal contact as a result of the polishing par- 

ticles picking up the metal from the surface. 


§ Vernon, Wormwell, and Nurse, Metallurgia 31, 19-22 
(1944). 
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The abrader consists of the following parts: 
(1) a glass plate over which a piece of polishing 
paper is tightly stretched, (2) a pair of connected 
precision ground parallel guide bars which can 
be shifted at will across the paper in definite 
strips, (3) a flat plate with suitable grooves for 
guiding on the parallel bars, (4) a brass chamber 
mounted in the center and perpendicular to the 
plate, and (5) a stainless steel piston in which 
the specimens are mounted. The abrading 
pressure can be regulated by adding weights to 
the stainless steel piston. After having been 
placed in the chamber the piston with specimen 
can take only two positions relative to the guide 
bars. This is to facilitate parallel and 90 degree 
abrading. The plate plus chamber and piston 
with specimen is pushed back and forth on the 
guide bars with the specimen in contact with the 
paper. 

The abrading treatment of a particular sample 
is determined by its hardness. Soft metals 
require different pressures and a more frequent 
change of polishing paper. The polishing papers 
used are 1, 0, 00, 000, and 0000. The last stage 
of the process is carried out exposing a new 
section of paper for each stroke and rotating 
90 degrees between each stroke. The abrading 
pressure is lowered during the last stage of the 
process. 

Using the above method we have successfully 
abraded pure aluminum and pure magnesium in 
the absence of lubricants. 

Following abrasion treatment the specimens 
can be given any pre-treatment desired before 
carrying out the experiment in the electron dif- 
fraction camera. The specimens are mounted in 
the instrument with the direction of abrasion at 
45 degrees to the electron beam. 


TYPICAL RESULTS 


Figure 6 shows the oxidation of pure electro- 
lytic iron at 300°C. The specimen was abraded 
with 1, 0, 00, 000, and 0000 polishing paper. 
Part (a) shows the diffraction picture of the 
metal taken by reflection at room temperature. 
Faint iron lines can be seen. Part (0) is the dif- 
fraction pattern obtained after 47 minutes of 
oxidation at 1.5 cm of air at 300°C. The pattern 
indicates Fe;O,. 

Figure 7 shows the oxidation of pure electro- 
lytic iron at 400°C and the effect of heat treat- 
ment on the resulting pattern. Part (a) is taken 
of the metal at room temperature while part (0d) 
is taken at 400°C in a vacuum before air is 
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Fic. 7. Oxidation pure iron 400°C 0.1 atmos. air. a. Polished No. 1 grit 25°C. b. Vac. at 400°C. c. Oxidized 60 min. 
air. d. Heated to 500° for 30 min. e. Cooled to 25°C overnight. 


added. The pattern is Fe;O,. Part (c) is the 
pattern resulting from 60 minutes of: oxidation 
at 400°C. and 0.1 atmos. of air. The pattern 
indicates a-Fe2sO;. The effect of temperature is 
shown in part (d) taken after heating to 500°C 
for 30 minutes. No change in the pattern is 
noticed. Part (e) is the diffraction pattern ob- 
tained after cooling to room temperatures. No 
change in the pattern is found. 

Figure 8 shows the oxidation of stainless steel 
(18-8) at 550°C. Part (a) shows the pattern ob- 
tained in a vacuum at 25°C. The lines can be fit- 
ted with either an iron or a chromium lattice since 
the lattice constants are similar. Part (b) shows 


»Fe or Cr a-Fe2O3 or Cr2O3 


a-Fe2Os or CroOs 


a pattern taken after oxidizing at 550°C in 7,5 
cm of air for 5 minutes. The structure is either 
a-Fe,O; or Cr2Qs3, since the lattice constants are 
very similar for the two oxides. Part (c) shows 
the same structure after continuing the oxidation 
for 50 minutes. To determine the equilibrium 
oxide present at 200°C the specimen was cooled 
to 25°C and then heated at 200°C for 16 hours. , 
Part (d) shows the resulting pattern which is 
the same as that found at 550°C. The specimen 
was now heated to 600°C. Part (e) shows the 
stability of the resulting oxide at 600°C. 
A complete study of the formation and sta- 

bility of the oxide films on the iron transition 


a-Fe2O; or Crz0s a-Fe20s or Cr203 





Fic. 8. Oxidation 18-8 stainless steel 550°C 0.1 atmos. air. a. Vac. 550°C. b. 5 min. air. c. 50 min. air. d. After cool- 
ing to 25°C and heating at 200°C for 16 hrs. (taken at 200°C). e. Same as (d) but taken at 600°C. 
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group metals and their alloys will be given in a 
later communication. 


SUMMARY 


An electron diffraction camera is described of 
special use in the study of physical and chemical 
reactions occurring on metals in special gas 
atmospheres and at elevated temperatures. 


Several furnaces are described for raising the 
temperature of the specimen. The chemical and 
physical reactions are carried out im situ and the 
specimen studied at the given temperature. Gas 
purifying chains are provided for the preparation 
of pure gases. 

The method of preparing samples is briefly 
discussed. Several typical results are given. 





Virtual Masses of Rectangular Plates and Parallelepipeds in Water 


YEE-TAak Yu 
University of Cincinnati, Cincinnati, Ohio* 
(Received June 26, 1945) 


The following is a report on experimental results obtained for the virtual masses of thin 


rectangular plates and rectangular parallelepipeds moving in water. 


Experimental data, 


together with an empirical equation, showing the manner in which the virtual mass of a rec- 
tangular plate varies with dimensions and direction of motion are given. Data and an empirical 
equation giving the virtual mass of a rectangular parallelepiped as a function of dimensions 
and direction of motion, for the special case in which the motion is parallel to one face, are also 


listed. 


INTRODUCTION AND SCOPE OF THE WORK 


T is a well-known fact that the force required 
to accelerate an object in a liquid is greater 
than the force required to give it the same 
* acceleration in space. This is an effect over and 
above any viscous or other dissipative force 
which may be acting. As a result, the effective 
mass of an object is greater in a liquid than in 
space. This increase in effective mass is referred 
to as “virtual mass.’’ The magnitude of the 
virtual mass of any particular object depends on 
the following factors: its size and shape, its direc- 
tion of motion, the density of the liquid, and the 
extent of the liquid in which it “is placed. The 
magnitude of the virtual mass of an object may 
be, and generally is, quite appreciable. Indeed, 
in certain cases, virtual mass may far exceed 
actual mass. 

To appreciate at once the importance of both 
experimental and theoretical work in this field 
it is only necessary to recall that the dynamical 
treatment of any mechanical system, which is in 
part or entirely submerged in a liquid, must 
involve a consideration of the virtual masses of 
- the submerged parts. 

Except for the consideration of a few simply 
shaped objects such as a sphere,! long cylinder, 


* Now at Research Laboratory, U. S. Naval Ordnance 
plant, Indianapolis, Indiana. 
! Horace Lamb, Hydrodynamics. 
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etc., very little theoretical work has been done. 
A search through the literature reveals that even 
less has been accomplished experimentally. 

The following is an account of the results of a 
series of experiments in which the virtual masses 
of thin rectangular plates and rectangular paral- 
lelepipeds in water were determined as a function 
of linear dimensions and as a function of direction 
of motion. The work herein reported is a con- 
tinuation of previous work? in this field. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


The apparatus consisted essentially of‘a simple 
torsion pendulum constructed as follows. A 
piano wire w, Fig. 1, having a diameter of 1 mm 
was stretched tightly between opposite legs of a 
strong metal frame F. A cross-arm AA consisting 
of a small brass rod fastened rigidly to the wire 
served as a support for a pair of identical rec- 
tangular plates or parallelepipeds BB. Large 
screws SS served to regulate the tension on the 
wire. A beam of light reflected from the small 
mirror m to a scale ten to fifteen feet away made 
it quite easy to obtain the period of oscillation 
of the arm for very small amplitude of motion. 
Since the entire frame was portable it was easy 
to obtain periods in air and in water. (Water was 

2 Yee-Tak Yu, “Virtual Masses and Me ya of Inertia 
of Disks and Cylinders in Various Liquids,” App: Phys. 


13, 66-69 (1942); ‘Double Torsion Sater le in a Liquid,” 
Am. J. Phys. 10, 152-153 (1942). 
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the only liquid used since theoretical considera- 
tion, as well as previous experimental work, had 
shown that all virtual masses are directly pro- 
portional to the density of the liquid.) The con- 
tainer for the water was a large cylindrical tank 
4 feet in diameter and 5 feet deep. The apparatus 
was submerged to a point where the mirror m 
was just above the surface of the water. 

The experimental procedure was simple. 
Identical rectangular plates or rectangular paral- 
lelepipeds BB were rigidly fastened to either end 
of the cross-arm of the torsion pendulum. 
Periods of oscillations, first in air and then in 
water, were carefully determined by timing a 
large number of. oscillations with an accurate 
stop watch. These two periods, together with 
certain other data, to be mentioned presently, 
were sufficient for computing the virtual mass of 
either object. 

The period of a simple torsion pendulum of 
the above type in air is given by 


Ta=2n(I./K)}, (1) 


and neglecting the damping due to the viscosity. 
(which on the basis of previous experience is quite 
negligible) the period of this pendulum in water 


‘ is given by 


T= 2(Iw/K). (2) 


Where J, and J, are the moments of inertia of 
the cross-arm with objects BB attached in air 
and in water, respectively, and K is the torsional 
constant of the wire w. Now let J,’ and I,’ 
represent the moments of inertia of the cross-arm 
only in air and in water, respectively. Then the 
increase in moment of inertia due to a single 
object B moving in water is given by 


T=43(Iw—Ia— (Tw! — Ia’) ). (3) 


This difference we shall refer to as the virtual 
moment of inertia of the object in the water. 
Letting Z represent the distance from the center 
of the wire w to the center of the object B and 
assuming that L is relatively long we can write: 


I=ML? (4) 


where M is now the virtual mass associated with 
B. Hence, assuming that the torsional constant 
K is known, Eqs. (1)—(4) give the final result 
we are seeking. It will be noticed that the same 
procedure can be used for finding virtual mass 
regardless of how the objects are oriented on 
the ends of the cross-arm. Thus this same 
apparatus and procedure was used as a means 
of determining the virtual masses of rectangular 
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Fic. 1. Diagram of apparatus. 


plates and rectangular parallelepipeds as func- 
tions not only of linear dimensions but of direc- 
tion of motion as well. 

The constant K was determined experi- 
mentally by replacing the cross-arm shown in 
Fig. 1 with a long uniform brass rod of known 
moment of inertia J;. J; is related to the period 
of oscillation T; in air by, 


K=4r(I;,/T}?). 


With this and previously given relations it is 
easy to derive the following convenient working 
equation, ; 


I 
M=——__{T,?—T,*-(T."*—T,')], (5) 
2T;7L? 





where T,,’ and T,’ represent the periods of the 
pendulum in water and in air respectively without 
objects BB. The values of T;, Ty’, and T,’ 
remain constant. However, it was found that 
relatively small changes in the tension would 
change these values slightly. Hence these 
periods were checked frequently. 

Rectangular plates for this work were cut 
from thin sheet lead. Rectangular parallelepipeds 
were cast from plaster of paris in a simple mold. 
After hardening and drying the plaster of paris 
casts were boiled in paraffin wax for some time 
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TaB_e I. Data for the parallelepipeds moving into water perpendicularly to their faces 2a X 2b 


(1k /Ti2 = 7.28 X 10°). 























2a X2bin em? 
T in sec. 
Lin cm 8x5 6x5 4x5 
M ing 
2¢ in cm. M, =133.2 M,; =90.4 M; =49.0 
Ta=3.72 T.=6.08 T.=3.32 Ty,=5.10 
0.15 L= 24.0 L=24.0 
M= 143.3 M =96.7 
M:= 10.1 M:2=6.3 
0.3 5.19 7.10 4.55 6.08 3.77 4.82 
23.9 23.9 23.9 
146.9 100.9 55.4 
13.7 10.5 6.4 
1.0 3.42 5.98 2.92 5.00 25.1 3.98 
23.3 23.3 23.3 
158.4 108.5 61.9 
25.2 18.0 13.0 
2.0 4.65 6.90 5.71 3.86 3.30 4.57 
23.3 23.3 23.3 
171.4 116.8 65.3 
38.2 26.4 16.3 
3.0 5.28 7.40 4.76 6.40 3.88 5.07 
23.3 23.4 23.3 
177.2 120.2 70.0 
44.0 29.8 21.0 
4.0 6.00 7.98 5.43 6.97 4.31 5.47 
23.3 23.3 23.3 
: 182.7 126.6 73.8 
49.5 36.2 24.8 
5.0 6.79 8.66 6.00 7.47 4.97 6.02 
23.3 23.4 23.3 
190.7 130.0 75.4 
57.5 39.6 26.4 
6.0 7.65 9.39 6.62 8.00 5.47 6.49 
23.4 23.3 23.3 
194.2 133.4 80.1 
61.0 43.7 31.1 
7.0 6.76 8.15 
23.3 
137.0 
46.6 
8.0 7.64 8.92 
23.3 
140.1 
49.7 








to render them impervious to water. A brass nut 
was cast in the face of each so that it could be 
screwed directly to the end of the cross-arm. 


EXPERIMENTAL RESULTS 


A. Virtual Masses of Rectangular Parallelepipeds 
With Motion Perpendicular to One Pair 
of Faces 


Following the procedure just outlined the 
virtual masses of three groups of rectangular 
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parallelepipeds were determined. Each group 
consisted of seven or more pairs of blocks of 
varying thickness. Each block of a particular 
group presented the same cross-sectional area 
(8X5, 6X5, and 4X5 sq. cm, respectively) per- 
pendicular to the direction of motion. Results 
are given in Table I. Letting 2a, 2b, and 2c 
represent the dimensions of a block (2c being 
the dimension parallel to the direction of motion) 
the empirical Eq. (6) gives a value of the virtual 
mass of any one of the blocks in close agreement 
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with experiment. A brief synopsis of the way in 
which this equation was arrived at will now be 
given. 

Graphs of Fig. 2 show the result of plotting 
virtual mass as a function of 2c for each of the 
groups. Consider any one of the curves. It is 
evident that it can be represented by M=M, 
+Mez where M, is the value obtained by extra- 
polating the curve to zero thickness. Mz is a 
function of thickness. (For any one curve 
2a X 26 is constant but 2c varies). A plot of log M2 
against log c for each family of blocks gives the 
straight lines of Fig. 3. The slope of each is 
approximately 3. Hence M, for any one of the 
curves, Fig. 2, is represented by M.=BC'. A 
comparison of the values of B, as obtained from 
each of the straight lines, shows that it is ap- 
proximately 3.5pab. The density of the liquid p 
is introduced here because it is known that in 
all cases virtual mass is directly proportional to 
this quantity. Thus we write: 


M=M,+3.5pabc'. (6) 


M,, the virtual mass for 2c approaching zero, 
undoubtedly depends on 2a and 2b. In order to 
find how M, depends on a and 6 the virtual mass 
of a number of lead plates, each having the same 
thickness (0.3 cm) but varying in area from 


M 


g 


g 


TOTAL VIRTUAL MASS /N GM 
S 
Ss 


/00 


80 


60 





40 
Q 





1 2 3 4 5 6 c 
THICKNESS (N CM 


Fic. 2. Virtual masses of rectangular parallelepipeds 
moving perpendicular to the faces of constant areas (8X5, 


— and 4X5 cm, respectively) as a function of thickness 
Ss 
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Fic. 3. Plot of log M2 against log C 
(Slope = approx. 4). 


4X4 to 4X12 sq. cm, was determined. Results 
are shown in Table II, column 1. (These same data 
are also an important part of Section B of this 
paragraph.) Since the plates are so thin, values 
of M are almost, but not quite, equal to M,. 
Assuming the correctness of Eq. (6) a value of 
M, was computed for each experimental value 
of M. From certain theoretical considerations to 
be mentioned later, it appeared likely that M, 
should be proportional to a’b?/(a?+5?)'. A plot 
of M, against this quantity gives the straight 
line shown in Fig. 4, the slope of which is 6.3. 
Hence we finally write as the total virtual mass 
of a rectangular parallelepiped having dimen- 
sions 2a, 2b, 2c and moving perpendicular to the 
face 2a X 2b, 


6.3a7b? 
M=o| ———+3.5abc\| (7) 
(a?+5?)} 


That this relation gives results in close agreement 
with experiment can be seen by comparing ex- 
perimental and computed values given in Table 
Ill. 


B. Variation of Virtual Mass of a Thin Rec- 
tangular Plate with Direction of Motion 


It is evident from elementary consideration 
that the virtual mass of a rectangular plate 
moving edgewise is less than when moving 
perpendicular to its face. Hence it is clear that 
virtual mass must in general depend on the 
direction of motion. In order to determine this 
dependence the following procedure was adopted. 
Exactly similar rectangular plates (the thin 
lead plates previously used) were rigidly fastened 
to either end of the cross-arm of the original 
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TaBLeE II. The variation of virtual mass of a rectangular plate with direction of motion. (In these data all plates have 
the same thickness, 0.3 cm. The effective arm L in each case was 23.55 cm. @ designates the angle between the 
normal to the face and the direction of motion.) 





























T in sec. C) 
M exp ing 0 w/i2 «/6 /4 x/3 54/12 x/2 
M comp in g 
2a X26 in cm? 
Re »  - 7. Ze To Te Te Te To Te To Te 
4x4 3.35 4.21 3.35 4.16 3.35 4.04 3.35 3.84 3.36 3.68 3.36 3.52 3.36 3.47 
pes Mexp Mexp Mexp am Mezp Mun 
40.3 37.0 31.0 20.7 12.1 4.6 1.8 
Meomp Meomp Meomp Meomp Mcomp M comp Memp 
41.0 38.3 31.2 21.4 11.8 4.4 
4x6 3.96 5.16 3.97 5.10 3.97 4.97 3.98 4.68 3.98 4.38 3.99 4.20 4.00 4.10 
70.2 64.7 54.8 37.6 19.6 7.4 2.7 
70.8 66.2 53.6 36.4 19.3 6.8 2.1 
4X7 4.28 5.64 4.28 5.58 4.29 5.36 4.29 5.07 4.30 4.74 4.30 4.48 4.30 4.40 
87.1 82.5 66.1 45.9 24.4 7.7 2.7 
87.1 81.2 66.0 44.6 23.4 YF a 
4x8 4.56 6.03 4.56 5.97 4.57 5.72 4.57 5.40 4.57 5.10 4.58 5.05 4.58 4.66 
102.5 96.2 76.1 52.6 28.2 9.1 2.9 
102.0 95.2 77.1 52.2 27.9 9.1 2.4 
4x12 5.46 7.38 5.47 7.27 5.48 7.00 5.49 6.55 5.50 6.08 §.51 5.74 5.52 5.62 
160.4 148.2 122.9 81.8 41.7 14.2 2.9 
159.4 147.2 120.7 81.2 42.2 15.0 3.0 


apparatus (Fig. 1) in such positions that the 
normal to the face of each was at a definitely 
known angle @ with respect to its direction of 
motion. This angle could be set to within an 


3 


MASS M, iN GM 
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(4) a, 7) 4§ 20 25 
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Fic. 4. Graph showing the relation between M, and 
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ator (where slope = 6.3). 
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accuracy of about +0.5 degree by means of a 
simple protractor arrangement which we shall 
not describe. A set of values of virtual masses 
for rectangules of various dimensions 2a and 28, 
but all with the same thickness (0.3 cm), for 
angles 0, 7/12, r/6, 2/4, 2/3, 54/12, and 2/2 
radians between the normal to the face and the 
direction of motion are given in Table II. These 
rectangles are so thin that they can almost be 
regarded as having only two dimensions. Plotting 
experimentally determined values of virtual 
mass against corresponding angles for each 
rectangle the curves of Fig. 5 were obtained. In 
plotting these graphs, values between 0 and 90 
degrees are repeated for corresponding angles on 
the other side of the normal. 

Any one of these curves is represented quite 
accurately, except for angle near +2/2, by the 
equation 


Me= Mo cos? 8. (8) 


That Eq. (8) does not give results in close 


TABLE III. Results obtained by computation from Eq. 
(7) and experimental results with rectangular paral- 


lelepipeds moving in water perpendicular to the face 
2a X 2b. 








2a X2bX 2c in cm 17x9xX4 8X4xX2 5X55 
M in g (Eq. 7) 1151 120 104 
M in g (exp.) 1148 118 104 
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agreement with experiment for angle near +2/2 
is owing largely to the fact that rectangles have 
a finite thickness. This fact will be taken account 
of in the work that follows. From the results just 
discussed it seems certain that, neglecting thick- 
ness, the virtual mass of a rectangle varies as 
the square of the cosine of the angle between the 
direction of motion and a normal to the face. 
Combining this with previous results we write 


6.3a7b? 
M= | cos? 6. (9) 
(a?-+b2)! 








Note agreement between Eq. (9) and experiment 
as shown in Table II. 


THEORETICAL CONSIDERATIONS 


The general theory of the virtual mass of an 
object moving in a perfect fluid has been com- 
pletely developed. Unfortunately, however, it is 
quite difficult to apply the theory to -most 
objects, even though their geometrical shapes 
may be rather simple. So far as is known, no 
theoretical work on rectangular plates or paral- 
lelepipeds has been reported. However, Horace 
Lamb does give M = pb’ as the virtual mass per 
unit length of an infinitely long strip of breadth 
2b when motion is perpendicular to its face. 
Hence, so far as agreement between theory and 
experiment is concerned, only two remarks can 
be made. The factor 6.3 occurring in Eq. (9) was 
obtained as the slope of the straight line Fig. 4. 
Now let a>} and @=0. Then neglecting } in the 
denominator of Eq. (9) and dividing through by 
2a we have M=3.15pb? which should represent 
the virtual mass per unit length of a very long 
strip. But this is exactly the relation given by 


Lamb. Secondly the cos? @ occurring in Eq. (9) . 


is to be expected from theoretical consideration. 
One of the boundary conditions over any surface 
is that the normal component of velocity of the 
fluid particles in contact with the surface must 


en Lamb, Hydrodynamics, Third edition, Chap. 
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Fic. 5. Virtual mass versus angle of motion. 


be equal to V cos @. Over a plane surface where 
@ remains constant V cos @ remains constant. 
Hence, the general expression for kinetic energy 
T will contain the square of this factor. 


CONCLUSIONS 


Insofar as comparison of experimental results 
with general theoretical considerations can be 
made, there is qualitative agreement. It is 
believed that the data herein given are suf- 
ficiently accurate to serve as a check on specific 
theoretical expressions which may later be 
reported. 

The empirical equations give results in close 
agreement with experiment and it is hoped, 
therefore, that they will be of immediate use in 
the field of applied science. 

The author wishes to express her thanks to 
Dr. D. A. Wells for many helpful suggestions 
and discussions. 
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Electron Microscopical Replica Techniques for the Study of 
Organic Surfaces 






R. BowLinG BARNES, CHARLES J. BURTON, AND RoBERT G. Scott 
Stamford Research Laboratories, American Cyanamid Company, Stamford, Connecticut 


(Received July 2, 1945) 


The problem of examining opaque surfaces in the electron microscope is discussed and it is 
shown that none of the methods heretofore described is completely satisfactory for the examina- 
tion of organic surfaces. Two new methods, a low pressure polystyrene-silica technique and a 
silver-silica double evaporation method, are described and it is demonstrated that they allow 
observation of organic as well as metallurgical surfaces. Electron micrographs of fiber surfaces, 
biological material, paint films, ores, and metals are included to illustrate the utility of these 


techniques. 





INTRODUCTION 


HE electron microscope has achieved con- 

siderable success in the observation of 
objects which are sufficiently thin or which are 
sufficiently dispersed to allow transmission of 
electrons of 30-100 kv energy. Bacteria, viruses, 
fibers, pigments, and smokes all make suitable 
subjects, and a large number of papers describing 
electron microscopical studies of these materials 
have been published. 

However, when efforts are made to observe 
“electronically opaque” specimens, as, for ex- 
ample, metallographic specimens, it is difficult 
to obtain microscopical data using the ordinary 
techniques. Typical metallurgical samples are 
far too thick to allow penetration of either light 
or electrons. The metallographic light microscope 
overcomes this difficulty by illuminating the 
specimen either from the side or through the 
objective and utilizes the light reflected by the 
surface to form the image. The same procedure 
is possible, at least in principle, with the electron 
microscope. Practical difficulties, however, have 
thus far prevented general utilization of this 
technique. 

Another method designed to accomplish the 
same purpose scans a small area of the surface 
with a very fine electron probe and reproduces 
the surface on the screen of a cathode-ray tube 
or facsimile recorder. This method was first sug- 
gested and developed into a practical instrument 
by Knoll' and von Ardenne.? More recently, 
Zworykin, Hillier, and Snyder* have constructed 
an instrument of this type which produces images 
having a resolving power of approximately 500A, 
considerably higher than that of the best optical 


1M. Knoll, Zeits. f. tech. Physik 16, 489 (1935). 

2M. v. Ardenne, Zeits. f. Physik 109, 553 (1938). 

*V. K. Zworykin, J. Hillier, and R. L. Snyder, A.S.T.M. 
Bull. No. 117, 15 (1942). 
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metallographs. However, the complexity of the 
equipment as well as the operating problems 
which were encountered have proved difficult to 
overcome. 

The limitations of direct methods of surface 
observation by electron optical methods led to 
the development of ingenious techniques for 
preparing thin film replicas of the surfaces of 
bulk materials. H. Mahl‘ was the first to suggest 
this possibility in connection with electron 
microscopical studies. He prepared thin oxide 
films on aluminum surfaces and floated them off 
after chemically etching away the base metal. 
He also suggested® the possibility of generalizing 
the method by using collodion, rather than 
natural oxide films, to form impressions of 
surfaces. 

A number of improved variations of this tech- 
nique have appeared during the past few years. 
Although similar in many respects, the methods 
vary chiefly in being either one-step or two-step 
processes. The direct, one-step replica technique 
consists in casting a thin film of resin on the 
prepared surface followed by mechanical strip- 
ping under water. This method, first described 
by Schaefer and Harker® of the General Electric 
Laboratories, -has been used by a number of 
workers because of the speed and simplicity of 
the technique. 

Two-step processes have been developed to 
eliminate the necessity for mechanical stripping 
of thin films, a procedure which is somewhat dif- 
ficult in any case and which is usually impossible 
when studying deeply etched or reentrant 
structures. The two-step method consists in first 
reproducing the prepared surface in a bulk 


4H. Mahl, Zeits. f. tech. Physik 21, 17 (1940). 

5H. Mahl, Jahrbuch AEG-Forschung 17, 67 (1940). 
6 ee Schaefer and D. Harker, J. App. Phys. 13, 427 
1942). 
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material which can be easily removed from the 
original and then making a thin replica of the 
first impression. The thin film is finally removed 
by some chemical method. 

Zworykin and Ramberg’ were the first to 
utilize this technique successfully in their silver- 
collodion method. A thick layer of metal, for 
example silver, is evaporated on the prepared 
surface by the usual techniques of vacuum 
evaporation and the metal film is then stripped 
mechanically from the original surface. A dilute 
solution of the final replica material, usually a 
collodion-amyl acetate solution, is flowed over 
the metal impression and allowed to dry. The 
metal film can then be removed by acid leaving 
a thin collodion replica which is examined in the 
usual way in a transmission electron microscope. 

The most completely satisfactory replica tech- 
nique to date is probably the polystyrene-silica 
method originated by Heidenreich and his col- 
leagues’ in the laboratories of the Dow Chemical 
Company. This is a two-step method in which 
the first replica is molded polystyrene; and the 
second thin film replica is evaporated silica. 
The polystyrene impression of the specimen is 
prepared in an ordinary molding press at a 
temperature of 160°C and with pressures in the 
range of 2000-5000 p.s.i. After cooling, the sides 
of the mold are sawed off and the original speci- 
men is either lifted off or dissolved in dilute acid 
or alkali. The surface is then exposed to silica 
evaporating in vacuum from a heated spiral tung- 
sten filament containing some quartz splinters. 
By dissolving the polystyrene base in some 
suitable solvent, this residual silica film may be 
utilized as a positive replica. 

This replica method is somewhat more difficult 
in technique than the others and may be a 
relatively time-consuming procedure. Good con- 
trast and resolution of the order of 100A are 
readily obtainable, however, and the replicas 
which result are easily duplicated and are rela- 
tively permanent. This method is perhaps the 
one most commonly employed at the present 
time and it has yielded outstanding results in 
many important metallurgical problems. 


LOW PRESSURE REPLICA FORMATION 


The restrictions attendant to all electren 
microscopical observations have been discussed 


7V. K. Zworykin and E. G. Ramberg, J. App. Phys. 12, 
692 (1942). 


8 R. D. Heidenreich and L. A. Matheson, J. App. Phys. 
15, 423 (1944). 


VOLUME 16, NOVEMBER, 1945 


at length in many places during the past few 
years. The action of an intense electron beam 
upon delicate specimens and the vacuum- 
induced dehydration effects are observational 
problems for which there as yet have been found 
no satisfactory solution. Consequently, the elec- 
tron microscopical studies of fibers and organic 
crystals which have been attempted have often 
been seriously circumscribed. 

In certain instances, partial solutions to the 
problem have been found. For example, in the 
study of cotton and wool fiber® it has been found 
that a mechanical disintegration of the material 
subdivides the parent fibers to such an extent 
that it is possible to observe the disintegrated 
specimen without any obvious artifacts arising 
from electron bombardment or vacuum desic- 
cation. 

However, in addition to obtaining information 
concerning the fibrillar structures, one is often 
desirous, for practical reasons, of observing the 
surface features of textile and paper fibers. 
Mechanical disintegration would in general 
destroy many of the features which might be of 
potential interest. This suggested the develop- 
ment of supplementary techniques which would 
allow surface examination of fibers, thereby 
augmenting the data obtained by the simple 
disintegration method. 

One fairly obvious approach to the problem 
is to consider the possibility of reproducing the 
surface of the specimen in a suitable plastic. 
Thin replicas suitable for electron microscopical 
studies could then be formed by one of the 
several available techniques. 

The chief modification necessary in an appli- 
cation of the polystyrene-silica replica technique 
to the study of organic, rather than metallic, 
surfaces is in a reduction of the pressure neces- 
sary to form the primary impression. Delicate 
organic materials would be crushed beyond all 
recognition, if one were to subject them to 
pressures of the order of 2000 p.s.i., as in the 
formation of metallurgical replicas. A reduction 
in the pressure, on the other hand, would neces- 
sarily entail a loss in resolution since the plastic 
would not be forced into all of the minute inter- 
stices of the surface. Nevertheless, it was decided 
to investigate low pressure methods in an effort 
to discover the possibilities and limitations of 
such a technique. Methods of this general nature 
have been successfully used in the past by 


®*R. B. Barnes and C. J. Burton, Ind. Eng. Chem. 35, 
120 (1943). 
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Fic. 1(a). Low pressure press used for forming initial 
polystyrene replica. 


EVAPORATING 
UNIT 














| SICA 


_—-# EPL /CA 








=<} 


. ~ i as 
SS 
S pre SS 
PUMP 








Fic. 1(b). Apparatus for vacuum evaporation of silica. 


optical microscopists'® and accordingly, the pos- 
sibility of producing impressions which would 
have adequate resolution to allow electron micro- 
scopical study was considered. After a number of 
preliminary experiments, a suitable method was 
devised and details of the technique are illus- 
trated in Figs. 1(a)—(d). 

The object under investigation is sandwiched 
with a }-inch square of 1-mil polystyrene film* 
between two ordinary glass microscope slides. 
The “sandwich” is then subjected to a pressure 
‘of 1 kilogram and is heated to approximately 
160°C. The heating was originally carried out in 
a small oven but it has recently been found that 
heating on a small hot-plate yields entirely satis- 
factory results. After cooling, the specimen is 
‘removed from the styrene impression either by 
mechanical stripping or by some suitable solvent. 


10 J. I. Hardy and T. M. Plitt, Wildlife Circular 7, U.S. 
Government Printing Office. 

*Obtained from Plax Corporation, Hartford, Con- 
necticut. 
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Fic. 1(c). Punch used for selecting desired portion of 
composite polystyrene-silica replica. 
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Fic. 1(d). Ethyl bromide boiler for removing poly- 
styrene from final silica replica. 


haatnesibstaell 


MOT PLATE 


Obviously, polystyrene is not the only resin 
suitable for formation of the original impression. 
In fact, other thermoplastic resins with lower 
softening points have in some cases been found 
preferable. One particularly useful material is 
isobutyl methacrylate, a resin having a softening 
point below 100°C, and consequently par- 
ticularly useful for forming impressions of the 
surfaces of water-containing substances. 

A film of silica a few hundred Angstroms in 
thickness is deposited on the polystyrene using 
the ordinary technique of vacuum evaporation." 
The desired portion of the polystyrene-silica 
composite may be selected and punched out with 
a tool which has an inside diameter slightly less 
than that of the ordinary specimen holder (ca. 
% inch), and the resulting disk placed in a 
specimen cap with a 200-mesh screen support. 
It has been found that the polystyrene, normally 


4 R. D. Heidenreich and V. G. Peck, J. App. Phys. 14, 
23 (1943). 
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Fic. 2(a). Polystyrene-silica replica of wool fiber showing 
scales and longitudinal striations. 9500. 





Fic. 2(b). Polystyrene-silica replica of human hair showing 
characteristic scale structure. 900 X< 





Fic. 2(c). Polystyrene-silica replica of human hair at 
higher magnification. Surface detail on scales is clearly 
visible. 9000 x. 
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Fic. 2(d). Polystyrene-silica replica of wool fiber showing 
gross scale features. 1200 


dissolved by immersing the composite replica in 
ethyl bromide, can be somewhat more effectively 
removed by placing the specimen holder in the 
condensing vapor of the solvent. 

-The cell used in this work consists of an ethyl- 
bromide boiler and a small condenser filled with 
ice water. The mount is placed in a region where 
the vapors condense, pick up some polystyrene, 
and drop back into the boiling flask. This brings 
the replica in contact with only pure ethyl 
bromide and minimizes fracture of the thin 
silica films by liquid pressure. Furthermore, once 
the polystyrene is dissolved, the silica film is 
properly oriented in the specimen holder and is 
ready without further handling for immediate 
insertion into the electron microscope. 

The first application of this new technique was 
made in a study of fiber surfaces. In order to 
show the entire width of the fibers, it was neces- 
sary to use relatively low electronic magnifica- 
tions. Even at these magnifications, however, the 
resolving power of the electron microscope is. 
adequate to allow considerable useful optical 
magnification. 

Figures 2(a)—(d) show electron micrographs of 
a group of wool fiber replicas. The scales and 
longitudinal striations are characteristic of fur 
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Fic. 3. Polystyrene-silica replica of nylon fiber showing 
“die-marks.”’ 1200X 


fibers and, judging from data obtained by means 
of the visual microscope, the reproduction is 
excellent. 

Surface replicas of nylon show no details of 
particular interest. This fiber is prepared by 
extruding the hot plastic through a small orifice, 
and it appears probable that the long grooves 
visible in the electron micrograph of Fig. 3 are 
die-marks resulting from imperfections in the 
spinnerette. 

Cotton is a fiber which has been exhaustively 
studied by many research workers. One of the 
more important microscopical features of cotton 
fibers is the cross-hatch structure formed by 
fibrils spiralling around the fiber axis. These 
spirals, one right-handed and one left-handed, 
form an angle of approximately 57° with one 
another, a value which has been substantiated 
by x-ray diffraction studies. A glance at Figs. 
4(a) and (b) clearly shows that the electron 
micrographs prepared from the silica replicas of 
cotton fibers fail to show this phenomenon, the 
structure shown by the replicas being more 
nearly reminiscent of a series of parallel single 
spirals. Further study, however, suggests two 
‘possible explanations. First, there may be a 
transparent membrane over the structure, or 
part of the structure, so that a plastic impression 
of the surface will not yield information directly 
correlative with the visual microscopical data. 
The second possibility is that a cotton fiber has 
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only a single series of parallel filamentary spirals 
wound around the fiber in “‘barber pole’’ fashion. 
If a visual microscope is focused at a point half 
way through the flat, ribbon-like cotton fiber, a 
slightly out-of-focus image will be obtained 
which will give the appearance of a double spiral. 
The fact that the angle formed by the filamentary 
bundles and the edge of the fiber is approximately 
30° lends considerable credence to this second 
theory. 

Thousands of electron micrographs of bacteria 
have been prepared during the past few years 
and it is known that a number of the electron 
microscopes in this country are being used 
exclusively on problems of this type. Neverthe- 
less, the limitations imposed by the high vacuum 
and electron bombardment to which the organ- 








Fic. 4(a). Polystyrene-silica replica of cotton fiber showing 
spiral nature of surface structure. 900. 








Fic. 4(b). Polystyrene-silica replica of cotton fiber showing 
structural features. 9000. 
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isms are subjected undoubtedly introduce arti- 
facts, although it is frequently difficult to 
estimate the extent to which these effects should 
influence interpretation of the photographs. 
Furthermore, the usual bacterial organisms are 
frequently too thick to allow adequate electron 
penetration so that internal morphology can be 
made visible. 

In view of this situation, the possibility of 
applying the replica technique to a study of the 
surfaces of cells and bacterial organisms was 
considered. The electron micrograph of red 
blood cells* reproduced in Fig. 5 shows one 
attempt in this direction. A blood smear was 
covered with a styrene film, was subjected to 
a force of 1 kilogram, and was then heated to 
150°C. After cooling, the plastic was pulled away 
from the smear leaving most of the cells behind. 
The micrograph shows the expected plate-like 





Fic. 5. Polystyrene-silica replica of red blood cells. 4000 x . 


appearance of the cells, together with a de- 
pression in the center. The irregular edges are 
probably caused by the film being forced under 
the cell while in the plastic state, then tearing 
when pulled away from the cell. Although the 
temperature and pressure to which the cells were 
subjected may seem excessively high, visual 
microscopical studies of preparations treated in 
this way fail to show any change in the cell 
morphology. A replica of B. Subtilis bacteria 
prepared in the same way is shown in Fig. 6. 


* A successful replica ef red blood cell surfaces was 
prepared by Mr. V. J. Schaefer of the General Electric 
Laboratories using the Formvar method and was described 
by him at the November, 1942 meeting of the Electron 
Microscope Society of America, Chicago, Illinois. 
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Fic. 6. Polystyrene-silica replica of B. Subtilis bacteria. 


The flagella appear as light lines throughout the field. 
9000 x. 





Efforts were also made to study heat-sensitive 
organic crystals by this new replica technique. 
In the past, numerous attempts have been made 
to examine vitamins and sulfa-drugs by the or- 
dinary methods of transmission electron micros- 
copy. These have proved generally unsuccessful 
for two reasons: (1) the crystals are usually so 
large that nothing more than silhouette images 
can be obtained and (2) the intense electron 
bombardment causes heating effects which often 
destroy the crystallographic identity of the 
materials. The replica method under discussion 
minimizes these problems, and the morphology 
of even large crystals is clearly reproduced. 
Figure 7 shows a replica of sulfadiazine crystals 
prepared following the method already discussed 
at length. There are also many metallurgical 





Fic. 7. Polystyrene-silica replica of sulfadiazine crystal 
showing crystallographic detail. 9000 X. 
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Fic. 8. Polystyrene-silica replica of pearlitic steel 
structure. 6000 X 


applications which may not warrant expenditure 
of the time and equipment required for suc- 
cessful utilization of the high pressure poly- 
styrene-silica technique."' In such instances, the 
low pressure polystyrene-silica method herein 
described yields results which are entirely satis- 
factory. Figure 8 shows an electron micrograph 
of a surface replica of pearlitic steel prepared by 
the low pressure method. This micrograph clearly 
indicates that many metallurgical problems will 
be amenable to study by this method. 

The use of diffraction grating replicas as a 
method of electron microscope magnification 
calibration, was first suggested several years ago 
in a publication” from these Laboratories. The 
simplicity and reliability of the method are well 
recognized by many workers in this field."* One 
rather serious drawback to more general use of 
this technique, however, is the fact that a 
relatively expensive original grating is required 
and this has in some instances proved a severe 
limitation. Sturdy nitrocellulose grating replicas, 
on the other hand, are inexpensive and widely 
used in spectroscopy and astronomical studies.* 
Accordingly, a method has been devised for using 
replica rather than original gratings. Silica is 
evaporated directly upon the ruled surface of 
such a grating replica. The nitrocellulose is then 
dissolved with acetone and the thin silica replica 
is ready for use in the electron microscope. 
Figure 9 shows a diffraction grating replica 


2C. J. Burton, R. B. Barnes and T. G. Rochow, Ind. 
Eng. Chem. 34, 1429 (1942). 

18 R. D. Heidenreich, J. Opt. Soc. Am. 35, 139 (1945). 

*Such replica gratings may be purchased, or may be 
prepared readily from an original grating, if available. 
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prepared in this manner. The permanence of 
these silica grating replicas makes them par- 
ticularly suitable for frequent recalibration of 
the electron microscope. 


SILVER-SILICA REPLICA TECHNIQUE 


Despite the success of the low pressure poly- 
styrene-silica replica method in allowing electron 
microscopical observation of a number of im- 
portant organic materials, certain difficulties and 
disadvantages of the method became apparent 
during the course of the work. In particular, only 
indifferent success resulted from efforts to obtain 
replicas of paper and texfile surfaces. Further- 
more, many problems may be encountered which 
require a study of surfaces of objects which in 
turn have a tendency to fuse with the thermo- 
plastic resin used for forming the impression and 
thus to destroy the very structures which it is 
desired to observe. 

Accordingly, alternative methods of replica 
formation were considered. One method in par- 
ticular seemed appealing: a combination of the 
Zworykin-Ramberg silver evaporation technique’ 
with the silica evaporation method of Heiden- 
reich.'* Such a double evaporation method would 
at once eliminate most of the objections which 
had been raised to the heat and pressure require- 
ments of the “low pressure polystyrene-silica 
method” described in the foregoing paragraphs. 
It would, of course, necessitate subjecting the 
original specimens to vacuum desiccation, but for 
many materials this would not be a serious re- 
striction. Moreover, the mechanics of the vacuum 

















Fic. 9. Silica replica of nitrocellulose diffraction grating 
replica. (11,000 lines/inch.) 7500. 


4 R. D. Heidenreich and V. G. Peck, Phys. Rev. 62, 292 
(1942). 
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evaporation technique should lead to higher 
ultimate replica resolution than with other 
methods where liquid surface tension effects 
play such an important part. 

In brief, then, the method can be described in 
the following way. A thick coat of silver is 
evaporated upon the surface of the specimen. 
This is then removed, either by mechanical 
stripping or chemical solution of the specimen. 
The desired portion of the primary silver replica 
is then selected and again placed in the evaporat- 
ing system where a thin coat of silica is deposited 
on the surface. The silver may then be removed 
by solution in dilute nitric acid and the silica 
replica examined in the usual way by means of 
the electron microscope. 

Figure 10 shows an electron micrograph of the 
surface of a sheet of alpha-pulp prepared by this 
double evaporation method. 





Fic. 10. Silver-silica replica of surface of alpha-pulp 
sheet. 900 X. 


Bacteriological specimens are also amenable 
to study by this method. A preparation of B. 
pseudomonas is shown in Fig. 11(a) and stained 
and unstained preparations of Escherichia coli 
are included as Figs. 11(b) and 1i(c). The 
increase in diameter of the stained flagella is 
obvious. 

Another problem of importance which would 
lend itself to this type of replica treatment is the 
study of paint surfaces. An electron micrograph 
of an enamel film is included as Fig. 12. 

As an illustration of the application of the 
silver-silica method to metallurgical studies, 
Fig. 13 shows a replica of the surface of a piece 
of stainless steel polished by the Chrysler Super- 
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Fic. 11(a). Silver-silica replica of B. pseudomonus. 9000 X. 





Fic. 11(b). Silver-silica replica of unstained Escherichia 
coli preparation. 7500. 





Fic. 11(c). Silver-silica replica of Escherichia coli stained 
with Liefson’s flagella stain. 7500X. 
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Fic. 12. Silver-silica replica of sprayed enamel 
film. 6000 . 


finish process.* The micrograph exhibits excellent 
resolution and appears to be entirely satisfactory. 

A polished surface of limonite ore has also been 
examined and the replica reproduced in Fig. 14 
obtained. 


INTERPRETATION OF REPLICAS 


Although much has been written regarding the 
mechanism of replica formation and the proper 











Fic. 13. Silver-silica replica of stainless steel surface 
polished by Chrysler Superfinish process. 6000 X. 


*Sample obtained from the Chrysler Corporation and 
photographed with their permission. The remarkable 
smoothness of this material is apparent from a considera- 
tion of the high magnification of the replica. The minute 
imperfections shown would be entirely invisible under the 
best visual microscope. 
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methods of replica interpretation, it is obvious 
that a considerable amount of research is still to 
be done in this direction. One optical illusion fre- 
quently occurs in replica photographs which to 
our knowledge has not been emphasized in the 
literature of electron microscopy and which has 
often been encountered in our own work. 

Thus, Fig. 15 shows two photographs of a 
diffraction grating replica, the only difference 
being that one of the positive prints has been 
rotated 180 degrees with respect to the other. 
Inspection of Fig. 15(a) clearly indicates a series 
of parallel grooves, whereas Fig. 15(b) certainly 
appears to comprise a set ef parallel projections. 
That this is not an isolated instance can be 








Fic. 14. Silver-silica replica of selectively polished 
limonite ore. 6000. 


proved by a careful examination of many of the 
electron micrographs appearing in this and other 
papers pertaining to replica studies. By rotating 
the photographs it is usually possible to cause 
the illusion of a complete reversal of topography, 
causing a depression to appear as a projection 
and vice versa. 

This effect may often be seen in current pic- 
tures of the war where aerial photographs of 
bomb craters are at times inadvertently reversed 
during the printing process and the craters 
appear to be projecting out of the earth. Rotating 
the picture will ordinarily produce the normal 
concave appearance. 

In order to illustrate the completeness of this 
illusion the photographs shown in Figs. 16(a) 
and (b) were taken. The object was a surface of 
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Fic. 15. Diffraction grating replica showing lines ap- 
pearing as (a) grooves and (6) projections. (30,000 lines/ 
inch.) 


sand in which a series of concave depressions had 
been made. This object was photographed ver- 
tically with a single source of illumination at 
approximately 45 degrees to the surface. The 
prints shown are identical save for orientation. 
A careful study of even the smallest details of 
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(a) 


(b) 


Fic. 16. Sand surface in which concave depressions have 
been formed. Photographic images appear concave in (a) 
and convex in (b), a complete topographical reversal having 
been effected by rotating the print through 180°. 


these photographsJ{clearly illustrates the phe- 
nomenon in question. 

Accordingly, unless considerable care is taken 
in preserving the proper orientation of electron 
micrographs of replicas, difficulties may be en- 
countered in interpretation. There are obviously 
solutions to this problem, perhaps the simplest 
method being the utilization of stereoscopy and 
contour finding.” 
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Stability of the Low Pressure Mercury Arc as a Function of Area at the Cathode Pool 


Carro_, W. Lurcy Anp Paut L. CoPpELAND 
Department of Physics, Illinois Institute of Technology, Chicago, Illinois 


(Received May 21, 1945) 


Using the method of Copeland and Sparing, the duration of an arc has been observed under 
fixed operating conditions a sufficient number of times to give a satisfactory estimate of the 
average life. The results show that the average life increases as the area of the cathode pool is 
reduced, and that the rate of change is particularly rapid for low arc currents. 





REVIOUS work! has indicated that an arc 

operating with a pool-type cathode has a 
finite chance of going out during each second of 
its operation. Conditions influencing the proba- 
bility of extinction are the temperature at which 
the cathode pool operates,’ the resistance in 
series with the arc,? and especially the arc 
current.! Work by Copeland and Sparing! has 
shown that conditions in the neighborhood of 
the cathode are primarily responsible for the 
extinction of the arc. The anchoring of the spot 
causes a very pronounced change in the’stability 
of the discharge. Although Copeland and Sparing 
had several anodes in their tube so that the 
geometry of the plasma and in particular the 
length of the discharge could be varied, they 
failed to detect a relationship between the 
stability and the length of the discharge. 

The method used in this investigation is 
basically that reported previously.'! The arc is 
repeatedly started and timed until it sponta- 
neously goes out. The average life, which is used 
as a quantitative measure of the stability, is 
obtained from a statistical analysis of the ob- 
served lives. Although the stability depends upon 
a great variety of operating conditions, the 
influence of each factor can be studied separately 
by varying it while maintaining all the others 
constant. 


EXPERIMENTAL ARRANGEMENT 


To obtain a convenient means for varying the 
cathode pool area, an experimental tube was 
constructed in which the mercury pool was con- 
tained within glass tubing which had been drawn 


down into the form of a cone as shown in Fig. 1. . 


Since the surface area depends upon the height of 
the mercury in the inverted cone, it is quite easy 
to change the cathode pool area by distilling 


1P. L. Copeland and W. H. Sparing, J. App. Phys. 16, 
302 (1945). sshitstitttied 
*W. Kaufmann, Ann. d. Physik [4] 2, 158 (1900). 
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mercury from a side tube into the main tube or 
vice versa. 

Previous work had shown that the stability of 
the discharge was influenced by the temperature 
of the bath in which the cathode pool is im- 
mersed, and the experimental arrangement con- 
tained a water bath which could be conveniently 
raised to surround the cathode or lowered to 
provide easy starting of the arc by means of the 
potential from a high frequency induction coil 
which was held against the glass wall of the tube 
just above the cathode pool. Throughout the 
experiments the bath was maintained at a tem- 
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Fic. 1. Experimental arc tube. 
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perature between 27°C and 28°C which corre- 
sponds to a vapor pressure of approximately 
0.0028 mm Hg. This temperature is high enough 
to provide easy starting, yet when the arc is in 
operation it is lower than that of other parts of 
the tube so that the mercury is retained in the 
cathode pool. The temperature is sufficiently 
close to the one at which Copeland and Sparing 
obtained maximum stability for a 3.0-ampere arc 
that rapid variation of stability with temperature 
is not to be expected in this range. 

Work by Ransom has shown the dependence 
of arc stability on circuit elements. A condenser 
placed across the arc reduces the stability, while 
an inductance in series with the arc increases it. 
The circuits used in these experiments may be 
understood by referring to Fig. 2. The objectives 
in the design of the circuit are simplicity and 
reproducibility. Circuits of two types were used 
in the course of the work. In the first of these an 
effort was made to reduce the influence of un- 
known conditions in the d.c. mains by placing a 
condenser of 88 mf capacitance across the arc 
and a series resistance of very low inductance. 
In this circuit it was possible to minimize the 
influence of any inductance in the meter meas- 
uring the arc current simply by placing it between 
the condenser and the mains in series with 
ordinary wire wound rheostats. The voltmeter 
across the arc had a resistance of 20,000 ohms/ 
volt, and it should not appreciably disturb the 
circuit. The resistors R; and R, directly in series 
with the arc across the condenser have negligible 
inductance. The wire coils have windings of the 
Ayrton-Perry type, and continuous adjustment 
is achieved by inserting a carbon compression 
rheostat R;. With the arc in operation the con- 
denser is maintained with a 75-volt potential 
difference. 


To achieve greater constancy in the supply, 
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Fic. 2. Circuit diagram. 


+R. A. Ransom, unpublished master’s thesis, Depart- 
ment of Physics, Illinois Institute of Technology (1941). 
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Fic. 3. Life data for 3.2-ampere arc. 


a second circuit was used for most of the work. 
A storage battery of 89-volt electromotive force 
and 90 ampere-hours capacity was connected as 
shown at 2. 


EXPERIMENTAL RESULTS 


Representative results in the life tests are 
shown by Fig. 3, in which the logarithm of N, 
the number of times that the arc survived for 
more than ¢ seconds, is plotted as a function of ¢. 
The arc current was 3.2 amperes. Previous work 
is confirmed in that the experimental points fall 
rather close to the straight lines drawn to repre- 
sent the data. The different lines represent data 
for different cathode pool areas indicated by the 
key in the upper right-hand corner of the figure. 
The negative of the reciprocal of the slope for 
each of these lines when divided by 2.303, the 
ratio of natural to common logarithms, gives the 
average life under the conditions of operation, 
and there is rather obviously a regular trend of 
the average lives as a function of cathode pool 
area. Similar results were obtained at 3.6 amperes 
of arc current. As would be expected from pre- 
vious work,' the average lives corresponding to 
the greater arc current are very much greater. 
A detailed examination of the data shows that 
in several cases the spread of the points is some- 
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Fic. 4. Life data with battery connected. 


what greater than would be expected on the 
basis of statistical fluctuations only. The natural 
interpretation would be that variation in the 
operating conditions was greater than it should 
have been. Statistical considerations‘ give the 
lower limit on the probable relative error to be 
expected in each set as 0.6745/N,', and since in 
most cases the probable error estimated from 
the spread of the points is quite close to that 
expected from fluctuation theory, the values of 
the probable error contained in Table I are 
based on statistical fluctuations only. 

When the 89-volt battery was placed across 
the condenser, the results shown in Fig. 4 were 
obtained with an arc current of 3.2 amperes. 
The trend of the average life with cathode pool 
area confirms that found in Fig. 3. Similar results 
were obtained with the arc current maintained 
at 2.8, 3.0, and 3.4 amperes. Table II sum- 
marizes all these data. 

The experimental tube described would not 
permit the use of very small cathode pool areas. 
For this reason a second tube, shown in Fig. 5, 
was constructed in which the cathode was con- 
tained in a capillary of area 0.03 cm’. There was 
such a marked increase in the stability of the 


*R. Peierls, Proc. Roy. Soc. A149, 467 (1935). 
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discharge for this tube that it was impractical 
to make average life determinations with the 
currents shown in Tables I and IJ. When this 
arc operates in a circuit of the first type without 
the storage battery connected, the average life 
is 8.7 sec. when the arc current is 1.7 amperes. 
There is an average life of 18.2 seconds for 1.8 
amperes and 27.4 seconds for an arc current of 
1.9 amperes. 


ANALYSIS OF THE RESULTS 


In view of the results presented there can be 
no doubt that tube geometry in the neighborhood 
of the cathode has a pronounced influence on the 
stability with which the arc operates. The natural 
interpretation is that when the cathode spot is 
free to move over the entire area of the pool the 
average life of the arc is a function of this area. 

The results obtained by Copeland and Sparing 
should correlate with those obtained in the 
present study. A cathode pool, having an area 
of 23.0 cm? which is quite a bit larger than any 
in the present study, was used. Notwithstanding 
the use of higher electromotive force and in- 
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Fic. 5. Experimental tube for small cathode areas. 
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ductive rheostats, both of which would tend to 
increase the stability, the observed average lives 
were a trifle lower than those reported here for 
corresponding arc currents and ambient tem- 
peratures. It would seem that anchoring the spot 
should be equivalent to a drastic reduction of 
the cathode pool area, and a spectacular increase 
in the average life was in fact observed when 
the cathode spot became anchored. These obser- 
vations suggest that, when the area is large, 
increasing it has a relatively small effect, but 
that very large average lives are obtained when 
the cathode spot is restricted to a small pool. 
The relevant data are summarized in Table III. 
In view of the low estimated accuracy of the 
results, the consistency of the trends shown is 
very satisfactory. 

A mathematical relationship in at least quali- 


TABLE I. Data on stability as a function of cathode pool 
area. (Potential difference across condenser, 75 volts, 
ambient temperature (27.5+0.5)°C.) 











Area of Probable 
Arc current cathode pool Average life relative error 
(Amperes) (cm?) (Seconds) (Percent) 
3.2 1.00 29.6 9.4 
3.2 1.42 25.9 5.0 
3.2 2.42 15.5 5.1 
ke 3.45 13.8 5.1 
3.2 4.5 10.1 5.0 
3.6 0.7 77.3 9.4 
3.6 1.1 64.0 9.4 
3.6 1.25 56.6 8.9 
3.6 2.00 47.8 6.7 
3.6 2.15 50.5 7.5 
3.6 2.72 47.6 4.8 
3.6 4.00 43.5 5.0 
3.6 4.5 saa 6.0 








tative agreement with this behavior is 
L=F(p, x, y, 2) exp (RI/A/), (1) 


where L is the average life of the arc the current 
in which is J, F(p, x, y, 2) is some function of 
the vapor pressure and of other undetermined 
factors represented by the parameters x, y, and z. 
A represents the area of the cathode pool, and 
through f(J) account is taken of the fact that 
the power of the pool area may depend upon 
the arc current. Taking the logarithm of both 
sides of Eq. (1) we have 


In [L/F(p, x, y, 3) ]=(RI/A™). (2) 
Repetition of the operation yields 


In In [L/F(p, x, y, 2)] 


=Ink+InI-—f(J)InA. (3) 
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Fic. 6. Variation of arc life with cathode pool area. 


In accordance with Eq. (3), it would be ex- 
pected that if In In L were plotted against In A 
in the relatively narrow range of A for which 
data were obtained, a linear relationship would 
be found. More precisely if In In L/F(p, x,y, 2) 
were plotted as a function of In A, the slope of 
the graph should give f(J) directly. Unfortu- 
nately, as is usually the case when a constant 
term multiplies an exponential, the value of the 


TABLE II. Data on stability as a function of cathode pool 
area. (Electromotive force of battery, 89 volts, ambient 
temperature (27.5+0.5)°C, probable relative error in aver- 
age life 6.7 percent.) 














Area of 

Arc current cathode pool Average life 

(Amperes) (cm?) (Seconds) 
2.8 0.65 16.5 
2.8 1.25 V2 
2.8 1.95 2.8 
2.8 4.9 2.25 
3.0 1.2 24.6 
3.0 1.86 15.4 
3.0 4.9 6.2 
32 0.8 28.2 
Sa 1.27 37.3 
32 1.83 17.7 
3.2 4.9 12.3 
3.4 0.9 46.3 
3.4 1.25 43.0 
3.4 1.85 37.0 
3.4 4.9 25.4 
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TABLE III. Variation of constants with area of 
cathode constriction. 





K F(p, x. y¥, 8) 
(amperes ~') (seconds) 


2.462 .06 0.011 
3.78 .40 — 
4.0 +1.0 
19.9 +0.5 


Reference 


C and S' : 
Cand L 4.9 

C and L 1.25 
C and S! Anchored 





0.00096 








constant F(p, x, y, z) cannot be determined with 
accuracy that is at all satisfactory. Perhaps the 
best that can be done is to assign to F(p, x, y, 2) 
the value 0.005 second, which is the order of 
magnitude expected from the work of Copeland 
and Sparing. (See Fig. 6.) Although the actual 
values obtained for f(J) depend upon F(p, x, y, 2), 


TABLE IV. Variation of arc stability with area 
at the cathode pool. 








Current 


Circuit type (amperes) 


1 
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the relative placement is not changed by large 
errors in this constant. (Calculations have been 
made setting F(p, x, y, z)=1, and F(p, x, y, 2) 
=10-*. The essential conclusions to be drawn 
from the calculations are the same in each case.) 
Using F(p, x, y, z) =0.005, the calculations yield 
for f(Z) the values presented in Table IV. 
Within the limits of experimental error inde- 
pendent sets of data taken in different circuits 
support the validity of Eq. (3), and they agree 
in supporting the qualitative behavior of f(J). 
The quantitative consistency happens to be 
better than the accuracy gives reason to expect. 


CONCLUSIONS 


The average life of a low pressure mercury arc 
with a free cathode spot is relatively large when 
the cathode pool area is small. Furthermore the 
life increases rapidly as the pool area is made 
very small. The variation of the average arc life 
with cathode pool area is rapid when the current 
is small and relatively slow when the current is 
large. This behavior may be represented within 
the accuracy of the present experiments by the 
relationship 


L=F(p, x, y, 2) exp (kI/A!). 
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